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Fast radio bursts (FRBs), millisecond-duration bursts prevailing in the radio sky, are 
the latest big puzzle in the universe and have been a subject of intense observational 
and theoretical investigations in recent years. The rapid accumulation of the observa- 
tional data has painted the following sketch about the physical origin of FRBs: They 
predominantly originate from cosmological distances so that their sources produce the 
most extreme coherent radio emission in the universe; at least some, probably most, 
FRBs are repeating sources that do not invoke cataclysmic events; and at least some 
FRBs are produced by magnetars, neutron stars with the strongest magnetic fields in 
the universe. Many open questions regarding the physical origin(s) and mechanism(s) of 
FRBs remain. This article reviews the phenomenology and possible underlying physics 
of FRBs. Topics include: a summary of the observational data, basic plasma physics, 
general constraints on FRB models from the data, radiation mechanisms, source and en- 
vironment models, propagation effects, as well as FRBs as cosmological probes. Current 
pressing problems and future prospects are also discussed. 
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|. INTRODUCTION 


Fast radio bursts (FRBs), milliseconds-duration radio 
bursts predominantly originate from cosmological dis- 
tances, are few remaining unsolved puzzles in contempo- 
rary astrophysics. The study of these mysterious events 
has a relatively short history. Some of the milestone ob- 
servations are highlighted as follows. 


e The first reported FRB was detected on July 24th, 
2001 (now called FRB 20010724), by the Parkes 
64-m telescope in Australia. It was not discovered 
until later by Duncan Lorimer and collaborators 
during an archival search for burst-like events. The 
discovery was published in 2007 in Science (Lorimer 
et al, 2007). This “Lorimer burst” was located 
3° from the Small Magellanic Cloud (SMC), had a 
peak flux density S, 2 30 Jy at ~1.4 GHz, a du- 
ration W ~ 5 ms, and a dispersion measure (DM) 
c 375 cm? pc. Its DM is in great excess of the 
value expected from Milky Way or SMC, suggest- 
ing that it likely originated from a cosmological 
distance. No repeated burst was discovered from 
the direction in 90 hours of follow-up observations 
(Lorimer et al., 2007). 


e Keane et al. (2012) reported another highly dis- 
persed burst-like event (later termed as FRB 
20010621A) with S, ~ 400 mJy at ~1.4 GHz, 
W ~ 7.8 ms, and DM ~ 746 cm"? pc. Since 
the burst was close to the Galactic plane, the ex- 
cess DM is not significant. The possibility that 
the burst was a giant pulse of an underlying pul- 
sar or from a Rotating RAdio Transient (RRAT) 
(McLaughlin e£ al., 2006) was not ruled out. 


e A strong support to the existence of extragalac- 
tic/cosmological FRBs was established (Thornton 


et al., 2013) when four more FRBs were discov- 
ered by the Parkes telescope. It was shown that all 
the events were from high Galactic latitudes, had 
large DM values in great excess of the MW values 
in those directions. Thornton et al. (2013) also es- 
timated that the event rate of FRBs is very high, 
about 10^ per day all sky above ~ 3 Jy ms fluence 
density threshold at 1.4 GHz. 


There have been FRB-like signals detected by the 
Parkes telescope that very likely have a terrestrial 
origin. These so-called “perytons” (Burke-Spolaor 
et al., 2011) were finally identified as artificial sig- 
nals caused during the magnetron shut-down phase 
of a microwave oven when a person impatiently 
opens the oven before heating is over (Petroff et al., 
2015c). Since those seemly-genuine bursts all hap- 
pened not during the dining time, this development 
suggested that FRBs are likely of an astronomical 
origin. 


Spitler et al. (2016) first reported in 2016 that 
one FRB source, named FRB 121102 (i.e. FRB 
20121102A, “R1”, or “Spitler burst"), emitted re- 
peated bursts with a similar DM as detected by 
the Arecibo 305-m radio telescope. The repeat- 
ing nature of the source allowed targeted observa- 
tions using the Karl G. Jansky Very Large Array 
(VLA) and the Arecibo telescope to detect addi- 
tional bursts and eventually localize the source us- 
ing the interferometric technique (Chatterjee et al., 
2017). This enabled the detection of a compact per- 
sistent radio source in association with the burst 
source (Chatterjee et al, 2017). Further very- 
long-baseline radio interferometric observations us- 
ing the European VLBI Network and the Arecibo 
telescope refined the persistent radio source to mil- 
liarcsecond scale, which corresponds to € 70 pc at 
the source (Marcote et al., 2017). It also led to di- 
rect identification of the source host galaxy in the 
optical band, which is a dwarf star-forming galaxy 
at redshift z = 0.19 (Tendulkar et al., 2017). This 
finally established the cosmological origin of FRBs. 


The polarization properties of FRBs have been 
studied closely over the years. Evidence of large 
rotation measure (RM ~ -186 rad m^?) in ex- 
cess of the Galactic value was reported for FRB 
20110523A, which suggested a dense magnetized 
plasma associated with the FRB (Masui et al., 
2015). The first repeater, rFRB 201211024, dis- 
played an even more extreme and evolving RM 
of the order 10°rad m^? (Michilli et al., 2018). 
It also showed an essentially non-varying polar- 
ization angle across each burst during individual 
bursts (Michilli et al., 2018). An opposite case was 
observed in another active repeating source FRB 


20180301A, which shows diverse polarization angle 
swings among different bursts (Luo et al., 2020b). 


rFRB 20121102A as the sole repeater did not 
last long before the Canadian Hydrogen Inten- 
sity Mapping Experiment (CHIME) discovered a 
few more repeating sources (CHIME/FRB Col- 
laboration et aL, 2019a,b). More repeaters were 
discovered through deep monitoring with the the 
Australian Square Kilometre Array Pathfinder 
(ASKAP) (Kumar et al., 2019) and the Five- 
hundred-meter Aperture Spherical radio Telescope 
(FAST) in China (Luo et al., 2020b; Niu et al., 
2021). 


Localizations of FRBs, both repeaters and non- 
repeaters, have been made via interferometry by 
the ASKAP collaboration, Deep Synotic Array 
(DSA) collaboration, and several other groups, 
which revealed a gallery of host galaxy types and 
positions of the FRBs within the hosts (Bannister 
et al., 2019; Bhandari et al., 2022; Macquart et al., 
2020; Marcote et al., 2020; Prochaska et al., 2019; 
Ravi et al., 2019; Xu et al., 2022) and the confir- 
mation of the theoretically expected DMjgm — z 
correlation (Macquart et al., 2020). 


An apparent ~ 16 day periodicity of a re- 
peating source, rFRB 20180916B (also called 
FRB 180916.701584-65), was reported from the 
CHIME observations (The CHIME/FRB Collabo- 
ration et al., 2020). Follow-up observations suggest 
that the active window is “chromatic”, with bursts 
detected in higher frequencies appearing at some- 
what earlier phases than those detected in lower 
frequencies (Pastor-Marazuela et al., 2021; Pleu- 
nis et al., 2021b). A tentative ~ 157-day period 
was also suggested for rFRB 20121102A (Rajwade 
et al., 2020). Bursting activities during the ac- 
tive windows are actually very sporadic. For rFRB 
20121102A, > 1600 bursts were detected by FAST 
in a total of 59.5 observing hours spanning 47 days 
during one active window (Li et al., 2021b), but 
there were no active bursts detected during some 
projected active windows later. 


On April 28th, 2020, an extremely high fluence, 
FRB-like event with two pulses was detected by 
CHIME (CHIME/FRB Collaboration et al., 2020) 
with one pulse also detected by the Survey for Tran- 
sient Astronomical Radio Emission 2 (STARE2) 
(Bochenek et al., 2020). The radio burst was as- 
sociated with a hard X-ray burst (HXRB) from 
a Galactic magnetar named Soft Gamma-ray Re- 
peater (SGR) J1935+2154 during one of its active 
phases (Li et al., 2021a; Mereghetti et al., 2020; 
Ridnaia et al, 2021; Tavani et al., 2021). This 


established a long-speculated connection between 
FRBs and magnetars. Deep monitoring of the mag- 
netar by FAST, on the other hand, suggests that 
the majority of X-ray bursts emitted by the mag- 
netar are actually not associated with FRBs (Lin 
et al., 2020), suggesting the rarity of the magne- 
tar FRB-HXRB associations. Deeper monitoring 
by FAST and European radio telescopes discovered 
fainter radio pulses from this source (Kirsten et al., 
2021; Zhang et al., 2020a). 


A repeating source rFRB 20200120E discovered by 
the CHIME FRB collaboration was found to be as- 
sociated with a nearby spiral galaxy M81 at a dis- 
tance of 3.6 Mpc (Bhardwaj et al., 2021). Follow-up 
observations surprisingly localized the source to a 
globular cluster in the host galaxy (Kirsten et al., 
2022). The bursts from the source have lower lu- 
minosities than typical cosmological FRBs. Some 
bursts have rapid temporal structures as short as 
60 nanoseconds (Nimmo et al., 2021). 


FAST detected another repeating FRB source, 
rFRB 20190520B (Niu et al., 2021), which resem- 
bles rFRB 20121102A in many aspects. Both 
sources have a point-source-like persistent radio 
source (PRS) associated with the burst engine and 
an extremely large and varying RM (Dai et al., 
2022; Michilli e£ al., 2018; Niu et al., 2021). rFRB 
20190520B, being located at z = 0.241 + 0.001, has 
the largest estimated host-contribution of DM ex- 
ceeding ~ 1000 pc cm" ?, which is likely due to the 
contribution of the PRS (Niu et al., 2021). 


Intense follow-up observations of the CHIME- 
discovered repeating source rFRB 20201124A by 
the FAST telescope (Xu et al., 2022) revealed pe- 
culiar short-term polarization property variations, 
including un-predictable RM evolution and non- 
evolution and oscillations of circular and linear po- 
larization degrees and linear polarization angles as 
a function of wavelength in a small fraction of 
bursts. Significant circular polarization was dis- 
covered from the source (Kumar et al., 2022b; Xu 
et al., 2022). Extreme RM variations, including 
a reversal of RM (Anna-Thomas et al., 2022; Dai 
et al., 2022), was observed from rFRB 20190529. 
All these suggest a dynamically evolving magne- 
tized environment around repeating FRB sources. 
Frequency-dependent polarization degree was no- 
ticed in a sample of repeating FRBs, which may be 
interpreted as a scatter of RM due to the multi-path 
propagation effect of radio emission (Feng et al., 
2022). 


One special source detected by CHIME, FRB 
20191221A, was identified to show a 216.8(1) ms 


periodicity with a significance of 6.50 (Chime/Frb 
Collaboration et al., 2022). It has a roughly 3 s long 
duration, making it an outlier in the FRB popula- 
tion. 


The physical understanding of FRBs also enjoyed a 
steady advancement in recent years, from knowing es- 
sentially nothing to painting a rough sketch of the FRB 
production mechanism. Similar to the field of gamma- 
ray bursts (Nemiroff, 1994), the early years of the FRB 
study also witnessed a large number of theoretical pa- 
pers dedicated to guessing the origin of FRBs based on 
very limited observational data (Platts et al., 2019). Not 
surprisingly, most of these ideas are quickly disfavored 
or completely rejected as data are accumulated. Rather 
than surveying all the proposed models (such a task has 
been carried out, see Platts ef al. (2019) and an online 
FRB theory Wiki page! ), this article focuses on a critical 
accessment of the leading ideas of interpreting FRBs that 
currently under active investigations. 


In the following, I will discuss the topics in the order 
of the decreasing certainty. I will first concisely summa- 
rize observational facts in 8II for the preparation of later 
discussion and refer the readers to more comprehensive 
observational reviews by (Bailes, 2022; Cordes and Chat- 
terjee, 2019; Petroff et al., 2019, 2022) and references 
therein?. After reviewing the basic plasma physics rele- 
vant to the FRB mechanisms (8III), I will discuss some 
generic theoretical arguments that pose constraints on 
any FRB models (8IV). The next section (§V) discusses 
possible mechanisms for generating the extremely coher- 
ent radiation of FRBs, with two general types of models 
(magnetospheric vs. relativistic shock models) discussed 
and compared in detail. This is followed by a survey of 
the source models (§VI) for repeating FRBs and some 
attractive ideas of generating genuinely non-repeating 
FRBs. The environmental models of FRBs are discussed 
un §VII and the propagation effects of FRBs are dis- 
cussed in SVIII. FRBs as various cosmological probes are 
summarized in 8IX. The review ends with a discussion 
of the problems and prospects in the field in 8X. Early 
theory reviews on the surveys of many theoretical mod- 
els can be found in Katz (2018b); Platts et al. (2019); 
Popov et al. (2018). Concise theory reviews on the phys- 
ical mechanisms of FRBs can be found in Zhang (2020c), 
Lyubarsky (2021) and Xiao et al. (2021). 


1 https: //frbtheorycat.org. 
? On the other hand, this review includes the most updated obser- 
vational progress not included in the previous reviews. 


Il. OBSERVATIONAL DATA 
A. Arrival times, coordinates, and naming convention 


A detected FRB is characterized by the time when it 
is detected on Earth (corrected to the barycentric time) 
and the spatial coordinate of the source. There have 
been different conventions to name FRBs. Since they 
are bursting events in nature, a widely adopted scheme 
is to name them based on the time when the burst was 
detected similar to gamma-ray bursts (GRBs), i.e. FRB 
YYMMDD. However, since some (probably most) FRB 
sources emit repeated bursts, one has to adopt the time 
when the first burst was detected to name the source. 
For example, the first repeater is widely named as FRB 
121102 or now officially FRB 20121102A. When CHIME 
came online, many detected FRBs flooded in. Since mul- 
tiple FRBs could be detected on the same day and some 
of them could be the repeating ones, the CHIME/FRB 
collaboration adopted a more informative/complicated 
name by combining time information and spatial infor- 
mation (right ascension [R.A.] and declination [dec]) of 
the source. For example, the second repeater detected 
by CHIME was named as FRB 180814.J04224-73 (now 
officially FRB 20180814A). There was also a suggestion 
to call repeaters as ‘R#’, where ‘#’ is an assigned num- 
ber based on the sequence of their discoveries. For ex- 
ample, FRB 121102 and FRB 180814.J0422+73 are also 
called ‘R1’ and ‘R2’, respectively. The 16-d periodic re- 
peater FRB 20180916B is ‘R3’. Another possibility was 
that one can add a prefix ‘r’ before the FRB name if a 
source is discovered to repeat. For example, ‘RI’, ‘R2’, 
and ‘R3’ may be also called ‘rFRB 20121102A’, ‘rFRB 
20180814A’ and ‘rFRB 20180916B’, respectively. There 
was an unofficial voting for the preferred naming con- 
vention among the attendees of the 2019-February FRB 
Workshop at Amsterdam, the Netherlands, but no con- 
sensus was reached. The commonly adopted naming con- 
vention in the literature now follows the Transient Name 
Server (TNS) convention ‘FRB YYYYMMDDabc’. Per- 
sonally, I think the information whether the source is 
a repeater is important. Throughout the review, I will 
follow the official TNS convention, but still add the ‘r’ 
prefix when the source is known to repeat. Other nick- 
names are also used occasionally. Note that the prefix ‘r’ 
is not a universally accepted convention but rather my 
personal preference. 


B. Temporal properties 


The typical observed duration (also known as width 
W) of an FRB is milliseconds. This duration is be- 
lieved to be the convolution of the intrinsic pulse du- 
ration at the source (W;), plasma scattering broadening 
(Tsc) during the propagation of the pulse, as well as in- 
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strumental broadening by the radio telescope (trei). As- 
sume uncorrelated Gaussian profiles of these components. 
The observed width may be written as (e.g. Cordes and 
McLaughlin, 2003; Lorimer and Kramer, 2012) 


W = [W2ü + 2)? 472 8], (1) 


where W; is the intrinsic duration of the FRB pulse in 
the source frame (the observed duration is longer by a 
factor of (1 + z) due to the cosmological time-dilation 
effect); 


1/2 


(2) 


is the scattering time, which includes the contributions 
from the Milky Way, intergalactic medium (IGM), and 


Tsc — [Tew + ToM + Tag (1 T z)"| 


the FRB host galaxy (see §VIII.A for a detailed discus- 
sion on scattering); and 


1/2 
tins = Ge samp F Atom F At?pw EE Ats,) (3) 


is the instrumental broadening (e.g. Cordes and 
McLaughlin, 2003; Petroff et aL, 2019), which in- 
cludes the data sampling interval, tsamp, the frequency- 
dependent smearing due to dispersion measure (DM), 


Atpm = (8.3 us) DMAvymzvai,. (4) 


the smearing due to the error of DM, Atspm, and the 
smearing due to the bandwidth, Ats, ~ (Av)! = 
1 HS (Avunz) | 

As shown in Figure 1, the lightcurves of FRBs show 
diverse behaviors. Many FRBs have one single pulse 
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FIG. 2 Examples of dynamic spectra of individual bursts from two repeaters that show down-drifting of pulses with frequency, 


also called the *sad trombone" effect. 


Upper: rFRB 20121102A (R1) (Hessels et al., 


2019). Red bars denote RFI excision; 


Lower: rFRB 20180814 (R2) (CHIME/FRB Collaboration et al., 2019a). 


(or indistinguishable multiple pulses). However, some 
FRBs (e.g. FRB 20121002) show an apparent temporal 
structure (Champion et al., 2016). The Galactic FRB 
20200428 had two pulses separated by roughly 30 ms, 
which may be also regarded as a repeating source that 
emitted two bursts. Some bursts clearly show an asym- 
metric pulse profile, with a longer decaying wing than 
the rising phase. This decaying wing is frequency depen- 
dent, with a longer tail at a lower frequency (e.g. FRB 
20111220, upper right panel, (Thornton et al., 2013)). 
The frequency-dependent scattering tail of these FRBs 
is consistent with Tse x v~* or Ts. x v~** as predicted 
by the plasma scattering effect (Cordes et al., 2016; Luan 
and Goldreich, 2014; Xu and Zhang, 2016). 


One interesting temporal feature of some FRBs is 


down-drifting of pulses with frequency (CHIME/FRB 
Collaboration et al., 2019a,b; Hessels et al., 2019), also 
called the “sad trombone” effect (Fig. 2). This is after 
correcting the standard dispersive delay due to propaga- 
tion and is likely related to the intrinsic radiation physics 
of FRBs. Such a behavior is often seen in repeating FRB 
bursts. The down-drifting is predominating. The oppo- 
site trend (up-drifting) has not been robustly claimed in 
the literature (CHIME/FRB Collaboration et al., 2021), 
even though the two apparently separated pulses in FRB 
20200428 indeed showed a higher peak frequency in the 
second pulse (CHIME/FRB Collaboration et al., 2020). 


The morphology of FRBs, especially for repeaters, has 
been studied extensively. Pleunis ef al. (2021a) stud- 
ied 536 bursts from 492 sources from the CHIME first 


catalog and identified four observed archetypes of burst 
morphology, namely “simple broadband,” “simple nar- 
rowband,” “temporally complex,” and “downward drift- 
ing”. Zhou et al. (2022a) studied more than 700 bursts 
from one repeating source rFRB 20201124A detected by 
FAST, and identified five morphological types based on 
the drifting patterns: downward drifting, upward drifting 
(a small fraction), complex, no drifting, and no evidence 
for drifting. Subtypes are introduced as needed based on 
the emission frequency range in the band (low, middle, 
high, and wide), and also the number of sub-pulses in the 
burst (1, 2, or multiple). Altogether, 18 morphological 
sub-types are identified. The longest burst includes 11 
pulses lasting 124 ms. There are no apparent correla- 
tions among duration, bandwidth, central frequency and 
flux. 


C. Spectral properties 


FRBs have been detected from 110 MHz (Pleunis et al., 
2021b) to at least 8 GHz (Gajjar et al., 2018). Non- 
detection at higher frequencies could be due to limited 
sensitivity (Law et al., 2017) or the difficulty to achieve 
strong coherence, but the lack of dispersion at high fre- 
quencies makes it difficult to differentiate RFIs from true 
signals, which might also contribute to the deficit. The 
non-detection at lower frequencies, especially with LO- 
FAR at 145 MHz, may suggest an intrinsic hardening of 
spectrum at low frequencies probably due to a certain 
absorption process (Karastergiou et al., 2015). 

'The spectral shape of some FRBs is not well measured. 
If one approximates the spectral shape as a power law 
function F, x v-*, the power law index a was observed 
to vary significantly from case to case. For example, 
the Lorimer burst had a = 4 + 1 (Lorimer et al., 2007), 
while FRB 20110523A had a = 7.8+ 0.4 (Masui et al., 
2015). Even for different bursts from the same repeating 
source, a can be very different. For example, the a values 
of rFRB 20121102A bursts ranged from —10.4 to +13.6 
(Spitler et al., 2016). Such a large variation may be the 
indication that the intrinsic spectrum of FRBs is nar- 
row. Multi-telescope studies of some repeater bursts of- 
ten show that the bursts detected in one band are not de- 
tected in another, e.g. for rFRB 20121102A (Law et al., 
2017) and rFRB 20180916B (Pastor-Marazuela et al., 
2021). This suggests that the spectra of these bursts 
are not simple power laws. Indeed, the dynamical spec- 
tra of FRBs (Figs.1 and 2) often show that the bursts are 
bright only in part of the whole observing bandpass. The 
Galactic magnetar burst FRB 20200428 had two pulses as 
detected by CHIME (CHIME/FRB Collaboration et al., 
2020), but only the second pulse that had a higher peak 
frequency was detected by STARE2 (Bochenek et al., 
2020), which has a higher bandpass than CHIME. This 
again suggests that the FRB spectra could be quite nar- 


row. A systematic study of the spectral properties of 
more than 700 bursts from rFRB 20201124A detected by 
FAST (Zhou et al., 2022a) suggests that the majority of 
repeating FRBs have narrow spectra, with the typical 
spectral band width of ~ 275 MHz. 


D. Repetition & Periodicity 


More than 20 FRBs have been reported to repeat 
(CHIME/FRB Collaboration et al, 2019a,b; Kumar 
et al., 2019; Luo et al., 2020b; Niu et al., 2021; Spitler 
et al., 2016). Since a repeating FRB is identified when- 
ever one more burst is detected from the same source, it 
is essentially impossible to claim that an FRB source is 
NOT a repeater. In fact, it is quite possible that all FRB 
sources repeat but with a wide range of repetition rate. 
Since the observed FRB rate density exceeds the rate 
density of supernovae, the most common catastrophic 
events, it is immediately inferred that the majority of 
the FRBs have to be from repeating sources (Luo et al., 
2020a; Ravi, 2019). The remaining question is whether 
all FRB sources repeat and whether there exists a mi- 
nority population of FRBs that do originate from catas- 
trophic events (Caleb et al., 2019a; Palaniswamy et al., 
2018). 

Some differences in the observational properties be- 
tween repeaters and apparent one-off FRBs have been 
noticed, but no conclusive results have been drawn. 


e The CHIME/FRB Collaboration (CHIME/FRB 
Collaboration et aL, 2021, 2019b; Pleunis et al., 
2021a) reported that repeaters tend to have wider 
widths than one-off FRBs. They also tend to 
have narrower spectra than one-off bursts. How- 
ever, the two populations have overlapping param- 
eter spaces, so that it is difficult to definitely tell 
whether an apparent one-off burst actually belongs 
to the repeater population. 


'The frequency down-drifting feature has been ob- 
served in several repeating sources (CHIME/FRB 
Collaboration et al, 2019a; Hessels et al., 2019). 
However, not all bursts from these sources and not 
all repeating sources show such a behavior. On the 
other hand, some apparently one-off FRBs show 
such a behavior, which may be regarded as candi- 
date repeating FRBs. 


Both supervised (Luo et al, 2023) and unsuper- 
vised (Zhu-Ge et al., 2022) machine learning algo- 
rithms applied on the first CHIME FRB catalog 
reached the consensus that repeaters and most non- 
repeaters seem to belong to different categories. 
Including both observed and derived parameters, 
both algorithms recognize brightness temperature 
and rest-frame spectral width as the two dominant 


traits to differentiate between the two categories. 
Some common candidate repeaters can be identi- 
fied from these two independent categories of ma- 
chine learning methods (Luo et al., 2023; Zhu-Ge 
et al., 2022). 


It is worth noting that some polarization properties, e.g. 
varying PA (Cho et al., 2020) or circular polarization 
(Dai et al., 2021), had once been proposed to be the 
unique properties of non-repeaters. However, later ob- 
servations showed that some repeaters also possess these 
properties (Luo et al., 2020b; Xu et al., 2022). It is now 
clear that polarization properties cannot be used to dif- 
ferentiate between the two categories. 

If all FRBs are repeaters, then at least some appar- 
ent one-off FRBs must have a very low repetition rate. 
Palaniswamy et al. (2018) and Caleb et al. (2019a) sug- 
gested that most FRBs cannot have a similar repeti- 
tion rate as rFRB 20121102A. Otherwise, many of them 
should have been observed to repeat. Indeed, extensive 
follow-up observations of some bright FRBs such as the 
“Lorimer burst” have so far failed to detect any repeated 
bursts (Lorimer et al., 2007; Petroff et al., 2015b), sug- 
gesting that they might have a different origin. Katz 
(2019) pointed out that the duty factor defined as D = 
(S)? / (S?) (S is flux density) may be used to differen- 
tiate repeaters from non-repeaters, with active repeaters 
such as rFRB 20121102A having D ~ 10^? while non 
repeaters having D ~ (1075 — 10719). 

With detailed simulations, Ai et al. (2021) suggested 
that tracking the evolution of observed repeater frac- 
tion Fp obs may shed light into the existence of genuinely 
non-repeating FRBs. This is because if genuinely non- 
repeating FRBs indeed exist, their numbers will linearly 
increase as a function of time. The number of repeaters, 
on the other hand, may approaching a limit with time. 
As a result, F; obs is expected to reach a peak and then 
decline. Therefore, detecting such a peak would strongly 
suggests the existence of genuinely non-repeating FRBs. 
In reality, however, depending on parameters and pos- 
sible evolution of source populations, the time to reach 
the peak could be long and the the duration at the peak 
could be also long. Long term monitoring of the sky us- 
ing CHIME-like wide-field survey telescopes will hold the 
key to place constraints on the existence of genuinely non- 
repeating FRBs. It is interesting to note that the recent 
CHIME observations suggested that F, obs stays constant 
for a few years already, which is consistent with the hy- 
pothesis that genuinely non-repeating FRBs do exist (Z. 
D. Pleunis, 2022, talk at the Cornell FRB workshop). 

Searches for periodicity of repeating FRB sources have 
been carried out extensively. The early targeted peri- 
ods in the searches were in the milliseconds to seconds 
range, similar to the periods of known pulsars and magne- 
tars. Deep searches of periodicity in this period range for 
rFRB 20121102A (Hewitt et al. (2022); Li et al. (2021b); 


Zhang et al. (2018b), see also an independent search by 
Katz (2022b)) and rFRB 20201124A (Niu et al., 2022; 
Xu et al, 2022) using thousands of bursts all led to 
null results, suggesting that FRB bursts are likely not 
giant pulses of rotating neutron stars. On the other 
hand, unexpected, very long periods (or active cycles) 
were found in some repeating sources. The most robust 
case is the CHIME-discovered rFRB 20180916B, which 
shows a ~ 16-d period with a ~ 5-d active window (The 
CHIME/FRB Collaboration et al., 2020). The duration 
and phase of the active window seems to be frequency- 
dependent, with the windows in higher frequencies ap- 
pearing earlier in phase and being narrower than the win- 
dows in lower frequencies (Pastor-Marazuela et al., 2021; 
Pleunis et aL, 2021b). Long-term monitoring of rFRB 
20121102A also revealed a possible long-term ~ 160-d 
periodicity (Cruces et al., 2021; Rajwade et al., 2020). 
Long-term monitoring of rFRB 20121102A with FAST 
suggests that bursts are often missing during the pre- 
dicted active window and the duty cycle of the periodic- 
ity becomes greater than 5096 (P. Wang et al. 2022, in 
preparation). This casts a shadow to the claimed period- 
icity. Finally, the *oddball" source FRB 20191221A was 
detected to have a 0.2168-s period with a significance of 
6.50 (Chime/Frb Collaboration et al., 2022). Since the 
total duration (~ 3 s) is much longer than other FRBs, 
this event likely has a different origin from the bulk of 
the FRB population. On the other hand, a deep pe- 
riodicity search of rFRB 20201124A bursts (Niu et al., 
2022) suggested that even though no global periodicity 
was found, fake local periodicity in adjacent burst clus- 
ters can be found with a significance up to 3.90. This 
cautions against claiming any periodicity from clustered 
bursts with significance < 4c. 


One interesting common feature of active repeaters is 
that the waiting time distributions of their bursts show 
two distinct peaks (Li et al., 2021b; Niu et al., 2022; Xu 
et al., 2022; Zhang et al., 2022; Zhou et al., 2022a). As 
shown in Figure 3, the first peak is around milliseconds 
and the exact value depends on how distinct bursts are 
defined. The second peak actually depends on the ac- 
tivity level of the source, ranging from 10s of seconds 
to 100s of seconds, even for the same source at different 
epochs. The bridge between the two peaks lie around 
10s of milliseconds. Since some FRB bursts show mul- 
tiple peaks, the short separations of bursts in the first 
component of waiting time distribution can be regarded 
as due to the similar origin as multi-peaks, which may 
be related to the continuous activity of the FRB source 
from one emission episode. Zhou et al. (2022a) defined 
“burst clusters" that include all the bursts whose rela- 
tive waiting times fall onto this first waiting time peak. 
'The second peak apparently scales with the global ac- 
tivity level of the source. More observations are needed 
to see whether the dip between the two components may 
carry information about the periodicity of the underlying 
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FIG. 3 The waiting time distributions of repeating FRBs. Left: The case of rFRB 20121102A during the 2019 active episode 
(Li et al., 2021b). The two peaks are at a few milliseconds and ~ 100 s, respectively. Right: The case of rFRB 20201124A 
during a four-day very active episode in September 2021. The second peak is at ~ 10 s, suggesting a very active episode (Zhang 
et al., 2022). In an earlier episode of the same source in April 2021, the second peak is ~ 100 s, suggesting that the same source 
can have very different activity levels, and hence, different waiting time distributions. 
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E. Dispersion measure and distance 


Radio waves in a plasma are dispersed, with waves 
with lower frequencies delayed with respect to waves with 
higher frequencies. The dispersion measure (DM) (see 
§III.B for details) describes the degree of such delay. The 
best-fit DM is obtained for each FRB when it is discov- 
ered?, and it carries the physical meaning of the column 
density of free electrons along the line of sight from the 
source to the observer (with the units of pc cm" ?). Since 
FRBs are from cosmological distances, the DM can be 
most generally written as 


D. mel) 


DM = A 
o 1l+2(l) 


dl, (5) 
where ne (a function of location denoted by 1) is the local 
electron number density, z is the redshift at that location, 


l is the comoving distance from the observer to a location 
along the path of propagation, and 


c [7 dz 


Da Jo BEY M" 


3 'The FRB search algorithm scans through a range of DM values 
to correct for such a delay. The DM of FRB is assigned either 
for the highest signal-to-noise ratio (S/N) or the finest burst 
temporal structure (Hessels et al., 2019). 


is the comoving distance from the observer to the source, 
where 


f(2) — exp E [ —M 208) 


142! 


Ho is Hubble constant, Qm, Qk and Qa are the energy 
density fraction of matter, curvature and dark energy, 
respectively, and w(z) = p(z)/p(z) is the dark energy 
equation of state parameter. For the concordance ACDM 
cosmological model, one has Qg = 0, Opg = Qa, w = —1, 
and f(z) =1. 

The observed DM is usually split into multiple terms 
(e.g. Thornton et al. (2013), Deng and Zhang (2014), 
Prochaska and Zheng (2019)) 


DMhost + DM; 
a a, 
(9) 
where DMww, DMnaro, DMicm, DMnost, and DMsre are 
the contributions from the Milky Way, its halo, the inter- 
galactic medium (IGM), the host galaxy, and the im- 
mediate environment of the source, respectively. No- 
tice that the observed contributions from the last two 
components are smaller by a factor of (1 + z), where 
z is the source redshift. The Milky Way term DMuw 
can be obtained using the MW electron density models 
derived from the radio pulsar data (Cordes and Lazio, 
2002; Yao et al., 2017) (with a > 50% uncertainty). The 
extended Milky Way halo contributes to an additional 


DM = DMyw + DMhnai + DMiam + 
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FIG. 4 The DM-z relation of FRBs with known redshifts. Left: The results of Macquart et al. (2020). The wide shaded region 
stands for the uncertainty of the Macquart relation due to density fluctuations of the large scale structure; Right: The latest 
DMz — z relation using 21 FRBs with redshift measurements. The NE2001 electron density model and DMhaio = 30pc cm? 
are adopted. The best-fit linear regression line is plotted. rFRB 20190520B has an abnormally large DMnost (Niu et al., 2021), 


which is marked separately. 


DMhalo ~ (30 — 80) pe cm? beyond DMyyw (e.g. Dolag 
et al., 2015; Prochaska and Zheng, 2019). 

The IGM component of DM is a function of redshift 
(Inoue, 2004; Ioka, 2003).The full expression reads (Deng 
and Zhang, 2014; Gao et al., 2014; Macquart et al., 2020; 
Zhou et al., 2014): 


. 8ceHoQfiam (^ xC + z')dz 
(DMiau(2)) = 81Gm; f E(z!) ; (10) 
where 
x(2) = Fxen(2)  gxeate(2) (11) 


noticing that the cosmological mass fractions of H and He 
are ~ 3/4 and ~ 1/4, respectively, Q is the energy den- 
sity fraction of baryons, frcm is the fraction of baryons 
in the IGM, and x,u(z) and Xe He(z) are the fractions 
of ionized electrons in hydrogen (H) and helium (He), 
respectively, as a function of redshift. The DM-z rela- 
tion is roughly linear at low redshifts (Inoue, 2004; Ioka, 
2003). With the standard cosmological parameters as 
measured by the Planck mission (Planck Collaboration 
et al., 2016), one can derive a rough linear relation at 
z « 3 (Zhang (2018a), see also Cordes et al. (2021); Pol 
et al. (2019)) 


Ho 
DMicm) ~ (855 a 
DN) = (Sip prta] (s km s-! gs) 


i (sans) (x) (3) i oo 


where friem is normalzied to ~ 0.83 (Fukugita et al., 
1998; Li et al., 2020b). In the literature, the DM-z rela- 
tion is also called the “Macquart-relation” to honor J-P 


Macquart’s leadership in the ASKAP collaboration to 
precisely localize a sample of FRBs and measure their 
redshifts to prove the theoretically motivated relation 
(10). Notice that Equations (10) and (12) apply to aver- 
age values. For individual FRBs, the measured DM can 
be either greater or smaller than the theoretical value due 
to the inhomogeneity of the IGM caused by large scale 
structures (Ioka, 2003; McQuinn, 2014). 

The redshifts of the localized FRBs (Table I) in- 
deed follow the theoretical expectations (Bannister et al., 
2017; Macquart et al., 2020; Marcote et al., 2020; 
Prochaska et al., 2019; Ravi et al., 2019; Tendulkar et al., 
2017). After deducting the Milky Way contribution, the 
external component of DM indeed shows a rough linear 
relation with z, with the best-fit line consistent with the 
prediction of the ACDM model (Macquart et al., 2020) 
(Fig. 4a). Using the Macquart et al. (2020) sample and 
systematically deducting an average DMpost value, the 
DM-z relation could give a constraint on fiam ~ 0.85 
(Li et al., 2020b), which is consistent with previous re- 
sults (Fukugita et al., 1998). 

Figure 4b gives the updated DMg — z relation with the 
21 redshift-known FRBs listed in Table I. The vertical 
axis is DMg = DM — DMyw — DMhnaio, where NE2001 
model and DMg4j, = 30 pc cm~? have been adopted. A 
simple linear regression best fit is presented. Using the 
YMW16 model or the average NE2001/YMW16 model 
lead to similar results, with slightly different regression 
results: 


DMg = 934.9z + 126.7, NE2001, (13) 
DMg = 979.7z + 103.1, YMWI106, (14) 
DMg = 957.3z + 114.9, average. (15) 
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TABLE I Published FRBs with measured redshifts, their observed DM values and the MW contributions. 


FRB z DM? DMmw(NE2001)? DMmw(YMWI16)° References 
rFRB 20121102A 0.19273 ~ 557 ~ 188 ~ 287 Tendulkar et al. (2017) 
FRB 20171020A 0.0087 ~ 114 ~ 37 ~ 25 Mahony et al. (2018) 
rFRB 20180301A 0.3304 ~ 517 ~ 152 ~ 254 Luo et al. (2020b) 
rFRB 20180916B 0.0337 ~ 349 ~ 199 ~ 325 Marcote et al. (2020) 
rFRB 20180924C 0.3214 ~ 362 ^ 41 ~ 28 Bannister et al. (2019) 
FRB 20181030A 0.0039 ~ 104 ^ 4l ~ 33 Bhandari et al. (2022) 
FRB 20181112A 0.4755 ~ 589 ~ 42 ~ 29 Prochaska et al. (2019) 
FRB 20190102C 0.2913 ~ 363 ~ 57 ~ 43 Macquart et al. (2020) 
rFRB 20190520B 0.241 ~ 1205 ~ 60 ~ 50 Niu et al. (2021) 
FRB 20190523A 0.6600 ~ 761 ~ 37 ~ 30 Ravi et al. (2019) 
FRB 20190608B 0.1178 ~ 339 ~ 37 ~ 27 Macquart et al. (2020) 
FRB 20190611B 0.3778 ~ 321 ~ 58 ~ 44 Macquart et al. (2020) 
FRB 20190614D 0.60 ~ 959 ~ 88 ~ 109 http://frbhosts.org/ 
rFRB 20190711A 0.5220 ~ 593 ~ 56 ~ 43 Macquart et al. (2020) 
FRB 20190714A 0.2365 ~ 504 ~ 39 ~ 31 Bhandari et al. (2022) 
FRB 20191001A 0.2340 ~ 508 ~ 44 ~ 31 Bhandari et al. (2022) 
FRB 20191228A 0.2432 ~ 298 ~ 32 ~ 20 Bhandari et al. (2022) 
rFRB 20200120E 0.0008 ~ 88 ~ Al ~ 33 Kirsten et al. (2022) 
FRB 20200430A 0.1608 ~ 380 ~ 27 ~ 26 Bhandari et al. (2022) 
FRB 20200906A 0.3688 ~ 578 ~ 36 ~ 38 Bhandari et al. (2022) 
rFRB 20201124A 0.0979 ~ 414 ~ 140 ~ 197 Ravi et al. (2022) 


a All DMs have the units of pc cm ?. 


b Calculated from the NE2001 model (Cordes and Lazio, 2002). Data provided by Ye Li who ran the script provided from 


https:/ /pypi.org/project/pyne2001/. 


© Calculated from the YMW17 model (Yao et al., 2017) using the website interface https://www.atnf.csiro.au/research/pulsar/ymw16/. 


Here the slope can be compared with the prediction in 
Eq.(12) and the y-intersection may be regarded as the av- 
erage (DMpost + DMgrc)/(1 + z). Comparing the fitting 
results to Eq.(12), one may tentatively draw the conclu- 
sion that friem > 0.9, which is greater than the estimate 
in the past (Fukugita et al., 1998). Considering the out- 
lier rFRB 20190520B with huge DMnost + DMsre (Niu 
et al., 2021) might have leveraged the y-intersection, an 
average value of DMpost + DM, ~ 100 pe cm~? would 
be reasonable. A systematically lower DMg than the lin- 
ear fit is noticeable at low redshifts, but this may be a 
result of large scale density fluctuations. More data are 
needed to judge whether there is a systematic deficit of 
DMg at low redshifts. 


F. Luminosity, energy and brightness temperature 


With measured redshifts, the isotropic-equivalent en- 
ergy and peak luminosity of FRBs can be measured pre- 
cisely. Because the DM — z relation has been confirmed 
from the data, for most FRBs without redshift measure- 
ments, the measured DM values can be used to estimate 
the redshift, and hence, the energetics of the FRBs. Lack- 
ing the geometric beaming information of FRBs, one can 
only estimate the isotropic-equivalent values of the peak 
luminosity and energy. The best estimates depend on 
the spectral shape of the FRB. If the FRB spectra are 
narrow-band with emission contained within the tele- 


scope observing band (which is the case for most bursts 
from repeaters, e.g. Zhou et al. (2022a)), it is more ap- 
propriate to multiply the bandwidth Av by the specific 
flux to obtain luminosity. On the other hand, if the FRB 
spectra are broad-band (which is relevant to some non- 
repeating FRBs, e.g. the Lorimer burst, Lorimer et al. 
(2007)) with emission extending beyond the telescope ob- 
serving band, it would be more appropriate to multiply 
the band central frequrncy ve by the specific flux to ob- 
tain luminosity (Zhang, 2018a). So, in general, one may 
write 


Av, narrow spectrum, 
Ve, broad spectrum, 
Di, j Sup (Av or ve) 
1028 cm / Jy GHz ’ 
(16) 


Lpiso & 4n D? Syp: { 


= (4n -10*? erg s7!) ( 


Eiso ~ 
SEC ds Ve, broad spectrum, 


|o 40-10% eg/ Dy V X 
7 1+2z 1078 cm / Jy-ms 


2 An Di. { Av, narrow spectrum, 
Es 


(Av or ve) 
GHz ' 
(17) 


where $,, is the specific peak flux density, F, is the 
specific fluence, and Dr, = (1 + z)D; is the luminos- 
ity distance. The isotropic peak luminosities of known 
FRBs vary from (Bochenek et al., 2020; Ravi et al., 2019) 
~ 10?5erg s~! to a few 1076 erg s~!. The corresponding 


isotropic energies vary from a few 10?? erg to a few 107? 
erg. The luminosity is extremely high by the radio pulsar 
standard, but is minuscule by the GRB standard. The 
true energetics of FRBs should be reduced by a beaming 
factor fẹ = max(AQ/4n,1/442) < 1, where AQ is the 
solid angle of the geometric beam, and y is the Lorentz 
factor of the FRB emitter (1/7 is the half kinetic beam- 
ing angle for an FRB emitter traveling close to speed 
of light). For an one-off FRB, a successful FRB engine 
should at least generate a luminosity and an energy of 
the order of f,L, and fy E, respectively. Observationally, 
the majority of hard X-ray bursts from SGR J1935+2154 
were not associated with FRBs (Lin et al., 2020). One 
possibility is that FRB emitters (at least those produced 
by magnetars) are narrowly beamed. If so, one would 
also expect to detect less-luminous but longer-duration 
radio bursts (“slow radio bursts”) with line of sight out- 
side the emission beam (Zhang, 2021). 

The combination of high luminosity and short variabil- 
ity timescale of an FRB defines an extremely high bright- 
ness temperature 7j. One may derive this by noticing 
that the observed specific intensity I, = S,/AQ, where 
S, is the observed specific flux, AQ = r(cAto)?/D% is 
the solid angle of the source viewed at the observer lo- 
cation (Ato = At/(1+ z) is the rest-frame duration of 
the burst, cAtg is adopted as the transverse scale, which 
is true for a non-relativistic, spherical, transparent emit- 
ter), and Da = D,/(1+ z) = Di /(1+ z)? is the angular 
diameter distance of the source. Considering an imagi- 
nary blackbody emitter with temperature Tp(vo) at the 
rest frame frequency vo = (1+z)v, and noticing I, (vo) ~ 
2kpTy(vo)(ve/c?) in the Rayleigh-Jeans regime (kp is the 
Boltzmann constant) and [,(vo) = I,(v)(1 + z)? (ie. 
I,/v? is constant), one finally obtains the brightness at 
the source frequency vo (Luo et al., 2023)4 


Tias) = S,DA(u-c-zp | S, D? 
bo) = Orkp(vAb?  Inkp(vAt)2(1 +z) 


S,( v AN? 
" 36 v 
Sey B ben) (=) 


x | (1 z) (Rg) > (18) 


ri Goss) 


The physical meaning of Tp is the imaginary temperature 
of the emitter if the photons and the electrons that emit 
the photons were in thermal equilibrium. This is appar- 
ently not the case for FRBs. The gigantic T, (~ 1099 K 


4 If the emitter is moving relativistically towards earth with a 
Lorentz factor I', the transverse size in the comoving frame would 
be l'cAt, so that T, is smaller by a factor of I? with respect to 
Eq.(18). The observer-frame Tj is boosted up by a factor of 
^ I, so the overall Tẹ is smaller by a factor of [ than Eq.(18) 
(see also Lyubarsky (2021)). Here we define Tẹ solely based on 
observables without assuming whether the source has relativistic 
motion. 
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for nominal FRB parameters) is much greater than any 
temperature allowed for incoherent radiation (see 8IV.E 
for details). This demands that the radiation mechanism 
for FRB emission must be “coherent”, i.e. the radiation 
by relativistic electrons is not only not absorbed, but is 
also greatly enhanced with respect to the total expected 
emission if electrons radiate independently (or incoher- 
ently). Before the discovery of FRBs, radio pulsars have 
been the only known sources of producing extremely high 
Ty's (typically ~ (1025 — 1099) K). FRBs further push the 
limit of the degree of coherent radiation in the universe. 


G. Polarization properties and rotation measure 


According to Petroff et al. (2019), early polarization 
measurements indicated a puzzling, heterogeneous pic- 
ture: the polarization properties can vary significantly 
among bursts. The high-quality polarization data ac- 
cumulated later suggested a more consistent picture: it 
seems that most FRBs have strongly polarized emission. 
The linear polarization degree is typically I, > 30%, 
sometimes nearly 100% (Cho et al., 2020; Day et al., 
2020; Luo et al., 2020b; Michilli et al., 2018). The ap- 
parent low polarization observed in some FRBs might 
be intrinsic, but could be also due to the large Fara- 
day rotation measure (RM, see Eq.(19) below) in these 
sources, as is the case of rFRB 20121102A (Michilli et al., 
2018). A frequency-dependent linear polarization degree 
has been observed in some FRBs, but it could be under- 
stood within a picture that the multi-path propagation 
effect introduces a scatter of RM so that the intrinsically 
strong polarization is smeared at low frequencies (Feng 
et al., 2022). Strong circular polarization has been ob- 
served in both apparently non-repeating FRBs (Caleb 
et al., 2018; Masui et al., 2015; Petroff et al., 2015a) and 
repeating FRBs (Kumar et al., 2022b; Xu et al., 2022). 
For linear polarization, the polarization angle (PA) re- 
mains constant across each burst for some FRBs (e.g. 
rFRB 20121102A, Michilli et al. (2018), see Fig.5 up- 
per panel). However, in some other FRBs, both appar- 
ent one-off ones (Cho et al., 2020) and repeating ones 
(Luo et al., 2020a), swings of PA across each burst are 
clearly observed, and the swing patterns are quite di- 
verse among bursts (Fig.5 lower panel). For the most de- 
tailedly studied repeater rFRB 20201124A, even though 
most of bursts are consistent with non-varying PAs, sig- 
nificant PA variations above 50 are observed in ~ 33% 
of bursts (Jiang et al., 2022). 

Linearly polarized radio waves propagating in a mag- 
netized medium would have the polarization angle under- 
going a frequency-dependent variation known as “Fara- 
day rotation". The degree of rotation is measured by the 
rotation measure defined by 


RM = (—0.81 rad m^?) [ [B (/nGIne(D (gj 
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FIG. 5 Examples of polarization angle variations across individual bursts from FRBs. Upper: constant PA in rFRB 20121102A 
bursts (Michilli et al., 2018); Lower: diverse PA swing patterns in rFRB 20180301A (Luo et al., 2020b). 


where Bj(I) is the [-dependent magnetic field strength 
along the line of sight (in units of micro-Gauss), n. is 
the number density of the medium along the line of sight 
in units of cm^?, and l is in units of pc. FRBs have 
a wide range of measured RM absolute values: whereas 
some of them have sizeable RMs ranging from a few hun- 
dreds to ~ 10? rad m^? in the case of FRB 20121102A 
(Michilli et al., 2018), some others have RMs consistent 
with being close to zero and could be used to place a con- 
straint on the magnetic field strength in the intergalactic 
medium (IGM) (Ravi et al., 2016). The distribution of 
RM/DM of FRBs, which gives a rough estimate of |B|, 
is slightly larger but not inconsistent with the distribu- 
tion of Galactic pulsars (Wang et al., 2020c). 


'The observed RM values of active repeaters show in- 
teresting variations. The first repeater rFRB 20121102A 
Michilli et al. (2018)) showed a secular decaying trend 
in RM. Short-term RM variation was observed in rFRB 
20180301A Luo et al. (2020a) and more clearly in rFRB 


20201124A (Xu et al., 2022). As shown in Figure 6 up- 
per panel, during an active episode of rFRB 20201124A, 
the RM of the source showed irregular variations during 
the first 36 days and turned to essentially invariant for 
another 18 days before the source quenched (Xu et al., 
2022). Another active repeater, rFRB 20190520B (Niu 
et al., 2021), showed an even weirder behavior. Its very 
large RM value of the order of 104 rad m^? underwent 
an unexpected reversal within 6 months (Anna-Thomas 
et al. (2022); Dai et al. (2022), see Fig.6 lower panel). 


H. Global properties 


The DM — z relation allows one to estimate the 
isotropic peak luminosity and energy of FRBs. For in- 
dividual sources, the estimated luminosity/energy can 
have a large error because of the uncertainty of the 
correlation. When a large sample of FRBs is consid- 
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FIG. 6 Examples of short-term rotation measure (RM) variations in active repeating FRBs. Upper: irregular RM variations 
observed during a 54-day compaign with FAST (Xu et al., 2022); Lower: surprising RM reversal from rFRB 201890520B 


(Anna-Thomas et al., 2022). 


ered, the uncertainties can be averaged out, so that the 
luminosity/energy function of FRBs can be reasonably 
studied. Independent groups (Hashimoto et al., 2022, 
2020; Lu et al., 2020; Lu and Piro, 2019; Luo et al., 
2018, 2020a; Zhang and Zhang, 2022; Zhang et al., 2021) 
reached the consistent conclusion that the bulk of the 
energy /luminosity function can be fit with a power law 
distribution: 


N(E)dE x E~®dE, N(L)dLe« L-"dL. (20) 


The index yg ~ Yr is not well constrained, e.g. 1.3- 
1.9 (Lu and Piro, 2019) or 1.5-2.1 (Luo et al., 2020a), 
but a central value 1.8 seems to be able to accommodate 
FRBs in at least 7 orders of magnitude, extending from 
~ 107° erg Hz-! for the Galactic FRB 20200428 to ~ 
1033 erg Hz-!, above which a possible exponential cutoff 
may exist (Lu et al., 2020; Luo et al., 2020a), see Fig.7 
upper panel. 

Besides global energy /luminosity distributions among 
FRB sources, for active repeaters one can derive detailed 
energy/luminosity distributions for individual sources. 
'The most comprehensive analysis has been done for a 
few active repeaters using FAST data. Li et al. (2021b) 
reported the detection of more than 1600 bursts de- 
tected from rFRB 20121102A in 47 days and found that 
there exist two components in the energy distribution. 
Whereas the high-energy part is consistent with a power 


law distribution, a distinct log-normal distribution com- 
ponent peaking at Eg ~ 4.8 x 10?" erg at 1.25 GHz is 
observed (Fig.7 lower panel). The energy distributions 
of rFRB 20201124A (Xu et al., 2022; Zhang et al., 2022) 
and rFRB 20190520B (Niu et al., 2021) show somewhat 
different shapes, but all require more complicated func- 
tions than the simple power law function. 


With the observed DM distribution, one can in princi- 
ple constrain the redshift distribution of FRBs. The ob- 
served DM distribution is the convolution of the intrinsic 
redshift distribution, FRB energy/luminosity function, 
and the instrumental fluence/flux sensitivity threshold, 
so inferring it is not straightforward. One needs to ap- 
ply a uniform sample (e.g. FRBs detected with the same 
telescope) to place the constraints. With the pre- CHIME 
data, Zhang et al. (2021) tested several astrophysically- 
motivated redshift distribution models, from a model as- 
suming FRBs tracking star-formation history to a model 
assuming FRBs tracking compact star merger events, 
which have a significant delay with respect to star forma- 
tion. They found that the available Parkes or ASKAP 
FRBs are not inconsistent with either model. James et al. 
(2022) showed that the simple non-evolution model is in- 
consistent with the data and found that the star forma- 
tion model is consistent with the ASKAP data. However, 
they did not test models invoking delays with respect to 
star formation. Hashimoto et al. (2020) suggested that 
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FIG. 7 Energy distribution of FRBs. Upper: The FRB 
isotropic energy distribution among different sources that 
shows a rough —1.8 power law distribution covering at least 
7 orders of magnitude. From Lu et al. (2020). Lower: 
The energy distribution of 1652 bursts detected from rFRB 
20121102A, which shows a bimodal distribution. From Li 
et al. (2021b). 


the limited data are consistent with no evolution with 
redshift. 

With the first CHIME-catalog, the FRB redshift dis- 
tribution can be further constrained. Zhang and Zhang 
(2022) pointed out that the DM distribution peaks at 
a value lower than predicted by the star formation his- 
tory model and suggested that the CHIME FRB data 
are consistent with a redshift model with a significant 
delay with respect to star formation. The conclusion was 
confirmed by Hashimoto et al. (2022) and Qiang et al. 
(2022), with the former group also claiming that the data 
are consistent with FRBs tracking the stellar mass rather 
than star formation rate. Using a reduced sample from 
the CHIME catalog, Shin et al. (2022) found that the 
CHIME bursts are still consistent with following the star 
formation history. However, this might be because Shin 
et al. (2022) have adopted criteria to remove low DM and 
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low S/N bursts, which have removed a significant number 
of nearby low-luminosity FRBs. However, those removed 
FRBs are the dominant population that demands a de- 
layed distribution from star formation. The existence of 
rFRB 20200120E in a globular cluster in M81 (Kirsten 
et al., 2022) suggests that such burst sources should be 
in abundance, which require significant delay from star 
formation. 


I. Host galaxies 


The first identified FRB host galaxy, that of rFRB 
20121102A, is a low-metallicity, dwarf star-forming 
galaxy, which is quite analogous to those of long- 
duration gamma-ray bursts (LGRBs) and superlumi- 
nous supernovae (SLSNe) (Nicholl et al., 2017; Tendulkar 
et al., 2017). On the other hand, the later identi- 
fied host galaxies, mostly for apparently non-repeating 
FRB sources, are typically Milky-Way-like massive spi- 
ral galaxies (Bannister et al., 2019; Bhandari et al., 2020; 
Heintz et al, 2020; Marcote et al, 2020; Ravi et al., 
2019). The positions of FRBs within the host galax- 
ies also carry clues for the origin of FRB sources. Even 
though rFRB 201211024 is located in an active star for- 
mation region of the host galaxy (Nicholl et al., 2017; 
Tendulkar et al., 2017), most other FRBs, especially ap- 
parently non-repeating ones, are not. Instead, many of 
them lie in the outskirt of the host galaxies with not 
particularly high star formation rate (Bhandari et al., 
2020; Heintz et al., 2020). The active repeater rFRB 
20201124A has a Milky Way-like massive host galaxy 
with high star formation rate (Fong et al., 2021; Piro 
et al., 2021; Ravi et al., 2022). Detailed observations 
with the Keck telescopes suggested that the host galaxy 
is a metal-rich, barred spiral galaxy, with the FRB source 
residing in a low stellar density, interarm region at an in- 
termediate galactocentric distance (Xu et al., 2022). This 
is inconsistent with the environment expected for long 
GRBs and superluminous supernovae. Cross comparing 
the host galaxy and FRB position properties with other 
astronomical transients, Li and Zhang (2020) showed 
that the global properties of FRBs are inconsistent with 
those of LGRBs and SLSNe, but are more consistent with 
Type II SNe and even compact object mergers. Overall, 
FRBs are not inconsistent with being all produced by 
magnetar engines, even though multiple formation chan- 
nels are also possible. Bochenek et al. (2021) compared 
the host properties of FRBs and core-collapse supernovae 
and reached the conclusion that the FRB environments 
are consistent with core collapse supernovae making mag- 
netars. 


J. Counterparts 


Most FRBs do not have counterparts detected in other 
bands or other messenger channels (e.g. gravitational 
waves and neutrinos). Searches have been conducted, 
and some putative counterparts were reported but not 
confirmed (e.g. DeLaunay et al., 2016; Keane et al., 2016; 
Sakamoto et al., 2021; Williams and Berger, 2016). So 
far, only two confirmed multi-wavelength counterparts 
have been observed for a few sources: 

First, both rFRB 20121102A (Chatterjee et al., 2017; 
Marcote et al., 2017) and rFRB 20190520B (Niu et al., 
2021) are found to be associated with a point-like per- 
sistent radio source (PRS). Incidently, these two sources 
are also active repeaters with relatively large RMs. It 
is suspected that all repeaters may have an associated 
synchrotron-emitting PRS (either a supernova remnant, 
a magnetar wind nebula, or a mini-AGN) but only the 
ones with a dense and highly magnetized environment 
(so a large RM) could be detectable (Yang et al., 2020a, 
2022a). 

Second, the Galactic FRB 20200428 (Bochenek et al., 
2020; CHIME/FRB Collaboration et al., 2020) detected 
from the magnetar SGR J1935+2154 was associated 
with a contemporary hard X-ray burst (Li et al., 2021a; 
Mereghetti et al., 2020; Ridnaia et al., 2021; Tavani et al., 
2021). Searches for X-ray/ y-ray emission in association 
with cosmological FRBs have been carried out for mul- 
tiple sources with null results (e.g. Cunningham et al., 
2019; Guidorzi et al., 2020; Laha et al., 2022a,b; Piro 
et al., 2021; Xu et al., 2022; Yang et al., 2019; Zhang and 
Zhang, 2017). The non-detection is expected since the 
predicted X-ray flux is below the sensitivity threshold of 
the detectors for cosmological FRBs even if the X-ray-to- 
radio luminosity ratio is the same as FRB 20200428. It 
is worth noting that there was a stringent optical upper 
limit (Z-equivalent 17.9 mag in a 60-s exposure) during 
the prompt epoch of FRB 20200428 (Lin et al., 2020). 
Since the prompt optical flux is very low even for the 
Galactic FRB, the chance of detecting a prompt optical 
counterpart for cosmological FRBs is slim. 

Searches for FRBs following some GRBs or superlu- 
minous supernovae have been carried out but with null 
results (Law et al., 2019; Men et al., 2019). Searches 
for progenitor explosions prior to some FRBs have been 
also carried out, with some candidates reported (Li et al., 
2022b; Wang et al., 2020d). 

Searches for gravitational waves (GWs) temporarily 
coincident with CHIME FRBs have been carried out, 
which led to tight upper limits on the GW fluxes (The 
LIGO Scientific Collaboration et al, 2022; Wang and 
Nitz, 2022). The null results imply at most O(0.01)% — 
O(1)96 of FRBs are associated with compact binary co- 
alescences (CBCs), which is consistent with the much 
higher rate density of FRBs than CBCs. Allowing a 
time difference between FRBs and GW events, a poten- 
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tial association pair between the NS-NS merger event 
GW190425 and a bright CHIME burst FRB 20190425A, 
with the FRB delayed by 2.5 hours with respect to the 
GW event, has been suggested (Moroianu et al., 2022). 
Its candidate host galaxy and the FRB environment are 
consistent with those expected for an NS-NS merger 
(Panther et al., 2022). 


IIl. BASIC PLASMA PHYSICS 


A plasma is a gas that contains a significant fraction 
of charged particles, usually with charge balance between 
negatively charged species (free electrons) and positively 
charged species (positive ions or positrons). An FRB is 
likely produced in a plasma and radio waves need to prop- 
agate through plasmas before reaching Earth. The dis- 
cussion of the physics of FRBs inevitably involves plasma 
physics, which we briefly review in this section. 


A. Plasma physics in the FRB context 


'The most important property of a plasma is the dou- 
ble reaction between particles and electromagnetic (EM) 
fields. While the EM fields would control the motion 
of the plasma, the motion of the plasma would generate 
currents and alter EM fields. The description of the phys- 
ical behavior of a plasma is therefore complicated (e.g. 
Kulsrud, 2005). In general, one needs to solve the evo- 
lution of the particle component (i.e. each species of the 
plasma) in six-dimensional (r, v) phase space in the form 
of the Fokker-Planck equation, and to solve the EM field 
component in three dimension in the form of Maxwell 
equations. For each particle and field component, one 
also needs to consider the physics in three scales: the 
large scale of smooth particle distribution and EM fields, 
the small scale of particle distribution and EM field varia- 
tions due to particle collisions, and the intermediate scale 
variation of particle distribution and EM fields dictated 
by various plasma waves. 

For the FRB problem, the most relevant scale is the 
intermediate one related to plasma waves. In many 
FRB radiation models, the observed FRB emission is 
related to certain types of plasma waves in the emis- 
sion region to begin with. The microscopic particle col- 
lisional/collisionless interaction processes are usually not 
important in interpreting FRB observations and we will 
not discuss them in the rest of the review. The largest 
macroscopic scale, on the other hand, could be impor- 
tant. This is particularly true if the emission region 
is from the magnetosphere of a rotating object (e.g. a 
magnetar), in which case the global magnetic field con- 
figuration and plasma density distribution play an im- 
portant role in defining FRB emission properties. For 
models invoking relativistic shocks, the globally ordered 


magnetic fields also play an important role in reproduc- 
ing some properties (e.g. high brightness temperature, 
high linear polarization degree) of FRB observations. 
More generally, radio waves associated with FRBs need 
to go through the plasmas between the source and the 
observer, undergoing dispersion, absorption, scattering, 
scintillation, and Faraday rotation and conversion for po- 
larized emission. In the rest of the section, we discuss the 
basics of dispersion and Faraday rotation and conversion, 
and leave more complicated multi-path effects (e.g. scat- 
tering, scintillation and plasma lensing effects) to Section 
VII. 


B. Radio wave propagation in a non-magnetized plasma 


Electromagnetic waves are oscillations of electromag- 
netic fields in both space and time in the form of 
exp i(k T — wt). When waves with a particular frequency 
go through a stationary plasma, even though their os- 
cillations in time (represented by angular frequency w) 
remain the same as in vacuum, their oscillations in space 
(represented by wave number k) would be modified in 
a frequency-dependent manner depending on the plasma 
properties. This leads to a varying wave propagation 
speed with frequency, known as dispersion. The rela- 
tionship w = w(k) is known as the dispersion relation. 

'The dispersion relation of EM waves propagating in a 
non-magnetized, globally neutral plasma can be straight- 
forwardly derived by introducing a space and time vari- 
ation of all quantities of the form of exp i(k - F — wt) in 
Maxwell's equations and Newton's second law equation 
involving the Lorentz force. The final dispersion relation 
reads (e.g. Rybicki and Lightman, 1979) 


ck? = ew, or ck nuu, or w? = ws -k?c (21) 
where 
4 2 
e=rni=l ud =1 (=) (22) 
iw w 


is the dielectric constant, n, is the index of refraction, c 
is conductivity defined by j = c E, and 


ATNee? id 
Up = ( : ) ~ (5.63 x 10* s 3)n1? (23) 
Me 
is the plasma frequency, where n, is the plasma den- 
sity, e and Me are the charge (absolute value) and mass 
of the electron, respectively. Noticing that e > 0 is 
required to have a real solution of the dispersion rela- 
tion w = w(k), one can see that w, defines a cutoff 
frequency, below which the EM waves cannot propa- 
gate. This is also the oscillation frequency of longitu- 
dinal waves (Langmuir waves) in a plasma?. If the FRB 


5 Note that the terms longitudinal (k || E) and transverse (k L E) 
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frequency (typically ~ GHz) is related to the plasma fre- 
quency, one requires n ~ (1.2x 101° cm?) vens 9, where 
VFRB,9 = vgnB/(10? Hz) and throughout the review the 
convention Qn = Q/10* is adopted in cgs units. 

The dispersion measure (DM) discussed in section 
(ILE) is defined through deriving the arrival time dif- 
ference of a pulse in two different spectral bands. One 
may start with the dispersion relation (21), which gives 
the group velocity of the dispersed wave 


2\ 1/2 
_ Ow | wr, 
vg(v) = aTe | — ] (24) 


This gives a frequeny-dependent arrival time of radio 
waves 
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where the approximation w >> wy has been adopted. The 
arrival time difference between two frequencies vo > nı 


can be expressed as 
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is defined. For a cosmological source, considering that 
the observed time tops = (1 + z)t and the observed fre- 
quency Yops = V/(1+ z), the final expression of DM is 
Eq.(5) when ¢ and v are expressed in terms of the ob- 
served values. Defining 


At ty, —t 
D= = 7 a 28 
Aap?) Ea” 28) 
1 2 
one can write 
DM = KD, (29) 


where (Kulkarni, 2020) 


K = 241.0331786(66) GHz ? cm? pes !, (30) 


indicate the direction of wave propagation with respect to the 
electric field E. EM waves are transverse waves. On the other 
hand, the terms parallel (k || B) and perpendicular (k 1 B) 
indicate the direction of wave propagation with respect to the 
magnetic field B. 


and DM is in units of cm^? pc. Notice that many as- 
sumptions have entered the above derivation (e.g. Kulka- 
rni, 2020): The motion of ions is neglected, the medium 
is cold, not moving with respect to the observer, and 
not magnetized. These factors are not important if the 
purpose is to give a rough estimate of electron column 
density along the line of sight but could be essential to 
perform precise measurements of arrival times and cross 
check the measurements of the same source by different 
detectors (e.g. the detection data of FRB 200428 between 
CHIME and STARE2). 


C. Radio wave propagation in a magnetized plasma 


1. General discussion 


When a plasma carries an ordered magnetic field B, 
the dispersion relation is much more complicated. Be- 
sides the plasma frequency w,, another characteristic fre- 
quency, the electron gyration frequency wg (also called 
Larmor frequency wr), is introduced. For non-relativistic 
motion, this frequency depends on B and fundamental 
constants, 1.e. 


TE ie RE (31) 


€ 


If the FRB frequency is related to wp, the required mag- 
netic field strength is B œ (360 G) vgng,s. Note that 
Qe = —eB/mec is defined as negative to contrast with 
the positive ion gyration frequency 


|. ZeB 
i mMc 


Q; = ZA., (32) 
Mii 
where m; is the mass of the positive ion and Z is the 
atomic number of the ion. For an electron-positron 
(ete!) pair plasma, one has Q; = |Q.| = wp. 
The existence of B introduces another special direction 

besides the wave propagation direction 

= Cs 

np = "i (33) 
'The dispersion relation becomes angle-dependent. Re- 
peating the exercise of wave expansion for the Maxwell's 
equations and Lorentz force equation for a global neutral 
plasma, one gets a dielectric tensor to replace the di- 
electric constant, which reads (e.g. Boyd and Sanderson, 
2003; Meszaros, 1992; Stix, 1992) 


» S —iD 0 
é=(iD sol. (34) 
0 OP 


'This is defined from 


ity, x (it x E) = e. Ë (35) 
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(which itself comes from the fourth Maxwell equation, 
j = dE with ë being the conductivity tensor, c is defined 
as @= [— An& [ ico, and I is the unit tensor), where E is 
the electric field vector of the waves, and the magnetic 
field direction is defined as the 2 direction. Here, 


S- (R+1)=1 ae oe (36) 
aan 00 
am coe a FAJ i 
S ud 
Bono 2j (40) 


where R, L, P denote parameters related to the “right”, 
“left”, and “plasma” modes, respectively, and S and D 
denote “sum” and “difference”, respectively. 

Very generally, 7, and B can have an angle 0. One can 
write 7, = (ny sin 0, 0, ny cos 0) without loss of generality, 
so that Equation (35) becomes (ii. E)ii, -n2 E+E E = 0, 
Or 


S—n2cos?0 -—iD  n2cos0sin0 Ez 
iD S-n 0 E, | =0. 
n2cos@sind 0 P—n?sin?6 E, 


(41) 
Taking the determinant of the coefficients, the general 
dispersion relation for waves propagating in a cold, mag- 
netized plasma becomes 


Ant — Bn? +C — 0, (42) 

where 
A= S sin? 0 + P cos? 6, (43) 
B = RLsin?0 + PS(1 + cos? 0), (44) 


C — PRL. (45) 


In the following, we consider the dispersion relations 
for a cold, magnetized plasma for different cases of the 
angle between k (or ñy) and B: 


2. k| B 
When the wave vector is along the magnetic field (e.g. 


for FRB waves propagating in the open field line region 
of a magnetosphere), Equation (41) becomes 


S-n? -iD 0 E, 
iD S-n? 0 E, | =0. (46) 
0 0 P E, 


Besides the P = 0 plasma mode (w? = w2), one has two 
transverse wave modes, i.e. the R and L modes’: 


n2 = R, R -— mode (47) 
n? = L, L — mode (48) 
with the dispersion relations 


w2 
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respectively. Note that hereafter for a pair plasma, the 
plasma frequency is defined as 


An. e? m 


in contrast to Eq.(23), where n+ = 2n, is the pair num- 
ber density, which is twice of ne for a neutral pair plasma. 
If one still uses the electron number density n, to define 
Wp, all the pair-related dispersion relations should have 
w2 replaced by 202. This is because in Eqs.(36)-(40), a 
small term w2; = 4ane(Ze)?/m; in parallel to w? has 
been ignored. This terms becomes comparable to w2 , in 
the case of pairs. 


Setting R = 0 and L = 0 and looking for positive 
solutions", one can define two cutoff frequencies 
Qi — Qe : We 9 +e 
vm [ag Rez A era 
MICE EE A 
B (we ub us. pair 
and 
NEPNCED INCL 
E Ed 4 AE. 
[dre hf —w/2, jon qus 
"d ua y/2, pair 
respectively. Here Q; < |Q-| and €; = |Q,| (Z = 1) 


have been adopted for an ion plasma and a pair plasma, 


6 Notice that opposite conventions of R-model and L-mode defini- 
tions have been used in different textbooks. For example, Boyd 
and Sanderson (2003) defines right(left)-handed with respect to 
the photon propagation direction while Rybicki and Lightman 
(1979) defines right(left)-handed with respect to the line of sight 
direction towards the source. We adopt the Boyd and Sanderson 
(2003) convention in the following discussion. 

7 Negative frequencies simply mean waves propagating in the op- 
posite direction. So solving positive solutions is complete in solv- 
ing the propagation problem. 
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respectively. The propagation condition for the R-mode 
and L-mode waves depends on the sign of the denomina- 
tors in Equations (38) and (39), respectively. 

Setting R — oo and L > oo, one can define two prin- 
ciple resonances at Wres.R = |Qe| = wp and wes, = Qi. 
'The frequency range that radio waves can propagate is 
defined by n? > 0, which is 


R — mode, (54) 
L — mode. (55) 


W Ug, Oro «Up, 


Up, orw < Qi, 
It is interesting to consider two asymptotic regimes. 


e In the regions far from the magnetosphere of a neu- 
tron star (e.g. in the ISM or IGM), one has w > wp 
Wp > wp and |Q.| > Q;. In this case, one has 
WR C wL ^ wp. The wave propagation condition 
is w > wy for both R- and L-modes, which is the 
same as a non-magnetized medium. 


e In the regions within a neutron star magnetosphere 
and for a pair plasma, one has w « wp, Wp X wp 
and |Q,| = iw = wp. In this case, one has wg ~ 
wy, ~ (w2+w3)"? ~ wp and the resonances are also 
WB- The R-mode and the L-mode become the same 
and are essentially transparent in all frequencies. 


3.k LB 


In another extreme case when the wave vector is per- 
pendicular to the magnetic field (e.g. for FRB waves 
propagating in the closed field line region of a magneto- 
sphere), Equation (41) becomes 


S —iD 0 E; 
iD S-n? 0 E, | =0. (56) 
0 0 P-n? E. 


One can also define two modes: the ordinary (O-) and 
the extraordinary (X- or E-) modes, i.e. 


2—P, O-— mode (57) 
L 
n= 2 X — mode, (58) 


with the O-mode dispersion relation 


2,2 w2 
E Z pP=1- 2 (59) 


w w2’ 


and the X-mode dispersion relation 


(w? —w2)?— ww, 
ek? RE Wa un-gzy: ion (60) 
BP CN paquete. air 
wu? —w2 3 P 
respectively. The O-mode corresponds to the case that 


the wave electric field vector is parallel to the background 


magnetic field vector, i.e. Ey || B, so that electrons mov- 
ing in response of E, oscillations do not feel the existence 
of the B field. As a result, the dispersion relation is the 
same as the non-magnetized medium case, and hence, 
the mode is called “ordinary”. The X-mode corresponds 
to the case of Ë, L B. The electrons in response of B. 
oscillations would also undergo gyration motion around 
the background B field, the hence, the mode is called 
“extraordinary”. The X-mode has cutoffs (k — 0) at 
wr (R = 0) and wr (L = 0), and principle resonances 
(k — co) at S — 0, which defines two (upper and lower) 
hybrid resonance frequencies 


a _ (opt +o 
Wres,H m D 


AM Me(w2 — 0,0.) V ^ 
t {1+ (we TOP 02)2 : (61) 
p i e 


For an ion plasma, since Q2 > Q2, it is interesting to 
note that the second term in the square root is always « 
1. One therefore has 


x 
= 
Q 


Whos, UH > we +F Q2, (62) 
0,0, w2 — Q0, 
Wes LH — 2 \ Z 2 2 (63) 
i we + OF + 02 
For an e* plasma with 0? = 0? = w2, one has 
Weos, UH a we + WR, (64) 
Wes, LH = WE. (65) 


The frequency range that radio waves can propagate 
(n2 > 0) is 


w > wp, O — mode, (66) 


W > UR, 
OT WL <W < Wres,UH; 
Or W < Wres,LH- 


X — mode (67) 


One can again consider two asymptotic regimes. 


e In regions far from the magnetosphere of a neutron 
star (e.g. in the ISM or IGM), one has w > wp, 
Wp > wp and |Qe| > 0;. In this case, one has wp ~ 
Wres, UH, WL ^ Wp and WrestH ~ wp. The wave 
propagation condition is w > wp for both O- and 
X-modes, which is the same as a non-magnetized 
medium. 


e In regions within a neutron star magnetosphere and 
for a pair plasma, one has w < wp, wp « wp 
and |Q.| = Qi = wg. In this case, one has 
WR ~ Wres, UH S5 wL = (w2 + wA)? ~ wg and 
Wres,LH = wp. So the X-mode is essentially trans- 
parent in all frequencies. The O-mode, however, 
can only propagate when w > wp. Because of 
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this, when radio waves propagate across the closed 
field line regions of a neutron star, the X-mode Ew 
vector would adiabatically rotate to maintain per- 
pendicular to the local B until reaching the radius 
where w > wp is satisfied, at which the polarization 
vector is frozen out (Lu et al., 2019). 


4. Oblique propagation 


When (k, B 
persion relation should take Equation (41), which is 
more complicated (not discussed below due the limited 
space, but see Boyd and Sanderson (2003); Stix (1992)). 
Nonetheless, the treatments in the two extreme cases are 
helpful to discuss the general behavior of the dispersion 
relations when 0 is small or close to 7/2: 


) — has an arbitrary angle, the dis- 


e When 0 < 1, one has the quasi-parallel case. The 
dispersion relations can be modified from the R- 
and L-mode relations (Eqs.(49) and (50)) by re- 
placing wg by wg cos 0. 


e When 6 — 7/2, one has the quasi-perpendicular 
case. The O-mode dispersion relation is revised to 


ek? wu? — we 
2 ~ 3 2 noa A? (68) 
w w? — We COS 0 


which can be reduced to Eq.(59) when 0 = 7/2. 
The X-mode dispersion relation can be modified 
directly from Eq.(60) by replacing wg by wg sin 0. 


In the literature, for the oblique cases, the X-mode and 
O-mode are usually defined as the cases when E, is per- 
pendicular and parallel to the (k, B) plane, respectively. 
Note that the O-mode defined this way is not completely 
“ordinary”, since there is still a E, component that is 
perpendicular to B. One should be cautious to extend 
the properties of the O-mode in the k L B case to the 
more general O-mode. For example, the statement that 
O-mode cannot propagate in a neutron star magnetop- 
shere is only valid in the quasi-perpendicular regime. In 
the quasi-parallel regime, even the *O-mode" is essen- 
tially extraordinary, i.e. a significant E, component is 
perpendicular to B. The waves can therefore also prop- 
agate. 


D. Faraday rotation 


Let us take a closer look at the propagation of radio 
waves in the case of k | B in an ion plasma. Dropping 
out 9;, the R(L)-mode dispersion relations (Equations 
(49) and (50)) can be generally written as 


2 2.2 wp 
TT EN M 
1+ (wp/w) 


ekeruüc "m (69) 


where the wp < w approximation has been adopted in 
the second equation, which is usually valid for the ISM 
and the IGM. 

Following the same procedure in SIILB and replacing 
ws by w;(1 +wg/w) (again valid for wg < w), one gets 


5 1/2 
vv) = h-3 082) l (70) 


w2 w 


Further requires wp < w, one can derive 


is [2 ji (48). (71) 
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D 
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(72) 


One can see that the effect of B field in the arrival time 
has a v~? dependence, which is much smaller than the 
DM term. It depends on i neBydl (a proxy of the ro- 
tation measure discussed below), but this term is practi- 
cally not measurable. 

A measurement of [d neDjdl is achievable by mea- 
suring the rotation of the polarization angle (PA) of lin- 
early polarized waves as a function of frequency known as 
“Faraday rotation". Since linearly polarized waves can be 
decomposed as the superposition of a right-handed and 
a left-handed circularly polarized components and since 
the two modes (R- and L-modes) have different propaga- 
tion speeds, the PA of the observed waves would display 
a frequency-dependent variation. Mathematically, this 
can be denoted as the variation of the phase difference of 
the circularly polarized waves as a function of frequency. 
Noticing kàc? = Rw?, k2c? = Lw?, and the phases of the 
R/L mode waves RL = i in x dl, the rotation angle 
can be written as 
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where 
e? " 
RM = —-——— e By dl 
E ou Pe 
y. ne By dl 
c (—0.81 a 4 
(—0.81 rad m^? mr ie (74) 


For cosmological sources, the observed wavelength is A = 
(1+z)Asr, so a more general expression is Equation (19). 


E. Faraday conversion 


More generally, Faraday rotation is a special case of 
^Faraday conversion". In general, a polarized electro- 
magnetic wave can be characterized by four Stokes pa- 
rameters (e.g. Rybicki and Lightman, 1979) 


E 

Q = ef cos 2 cos 2x, 
US Ea cos 21 sin 2x, 
V = ee sin 2y, 


where I = 4/Q? + U? + V? is the total intensity, L = 
\/Q? + U? is the intensity of the linear polarization, V 
is the intensity of the circular polarization, eg = VI is 
the amplitude of the elliptically polarized EM waves, v» = 
(1/2) arcsin(V/T) is a proxy of the circular polarization 
degree II, = V/I which is intrinsic to the waves, and x = 
(1/2) arctan(U/Q) is the angle between the semimajor 
axis of the ellipse and the z-axis defined by the telescope, 
which is extrinsic to the waves. Notice that (I, 2v and 
2x) are spherical coordinates in a imaginary Poincare 
sphere, and (Q, U, V) defines a polarization vector P 
from the center to a point on the sphere in the Cartesian 
coordinate system, which defines the polarization state of 
the wave. Faraday rotation is simply the rotation of the 
P vector around the V axis. When P rotates around axes 
other than the V axis, there would be conversion between 
linear polarization L and circular polarization V. The 
waves would undergo Faraday conversion (Melrose et al., 
1995; Zheleznyakov and Zlotnik, 1964). 

The physics of Faraday conversion can be understood 
as follows. Any polarization state can be decomposed 
into superposition of two fundamental modes, either two 
circular polarization modes (e.g. R- and L-modes) for 
the quasi-parallel case or two linearly polarization modes 
(e.g. O- and X-modes) for the quasi-perpendicular case. 
'The different phase velocities of the two eigen modes 
would make the two modes out of phase and introduce 
modified polarization behaviors after superposition. For 
the quasi-parallel case, the different velocities of R- and 
L-modes introduce rotation of the superposed linear po- 
larization angle, and hence, Faraday rotation. For the 
quasi-perpendicular case, on the other hand, the differ- 
ence in the propagation velocities in the O- and X-modes 


would make the two modes out of phase, making the su- 
perposed polarization elliptical. Effectively, part of linear 
polarization is converted to circular polarization. The 
amplitude of Faraday conversion is smaller than that of 
Faraday rotation by a factor of wp/w, which is < 1 for 
waves propagating in a medium far outside of the neutron 
star magnetosphere. 

Mathematically, one may consider that the vector P 
undergoes rotation around an imaginary vector axis in 
the direction of 


Q = (g,h, f) (79) 


on the Poincare sphere. The variation of the circular 
polarization degree can be described by dP /dz = Q x P, 
where the z-axis is the direction of the V component 
(Gruzinov and Levin, 2019). The three components of Q 
are 


f=-- - Bz, (80) 
C W 
2,42 . » 
h+ig=-—-— "3 (Be + iByy, (81) 


where (Bz, By, B=) is the unit vector Ê = B/B, f de- 
notes the traditional Faraday rotation rate discussed in 
Equation (73), and (h 4- ig) describes the Faraday con- 
version rate. To order of magnitude, one can see that 
h/f ~ g/f ~ wg/w, which is « 1 for waves propagating 
far outside a neutron star magnetosphere. This means 
that © is essentially parallel to the V direction and that 
Faraday conversion is a small-order effect compared with 
Faraday rotation. 

If one measures oscillations of Stokes parameter V , one 
may define a conversion measure (CM) as (Gruzinov and 
Levin, 2019) 


(Iv) = CM)?, (82) 


where IIy = |V|/J and (IIy) is the rms value of Iy. The 
CM can be related to RM through 


CM ~ “8 RM"? ~ (107?m-?) RMM/?(B/G), (83) 
Ww 


where B is in units of Gauss and RM is in units of 
rad m" ?. This is strictly valid for a small conversion 
angle 0; (the final angle by which the linear-polarization 
Q — U plane rotates). For a large 0j, a more precise 


expression is (Gruzinov and Levin, 2019) 


(Iiv) = "E [e-(€M »2)?/2 — e(CM XY], — (84) 


When both CM and RM are measured, one can directly 
measure B using Equation (83). 

Physically, for a cold plasma Faraday conversion hap- 
pens when the B field is quasi-perpendicular. Astrophys- 
ically, this may be (but is not necessarily) related to the 


22 


reversal of Bj along the line of sight (Gruzinov and Levin, 
2019; Melrose, 2010). Another possibility of having Fara- 
day conversion is when electrons are no longer “cold” but 
are mildly relativistic with a mean Lorentz factor ye > 3. 
This is because when considering the response tensor or 
electrons with a general energy distribution, the expres- 
sions of the h, g, and f parameters depend on the ye in 
the medium (Huang and Shcherbakov, 2011). As ye in- 
creases, h and g increase and f decreases so that conver- 
sion becomes progressively more important and rotation 
becomes less important. The non-detection of Faraday 
conversion in rFRB 20121102A has been used by Vedan- 
tham and Ravi (2019) to place an upper limit on the 
Lorentz factor of the electrons in the medium that gen- 
erate the conversion, i.e. Ye < 5. 

Faraday conversion can be more generally described 
using the transport equation (Huang and Shcherbakov, 
2011; Li et al., 2022a) 


B €r n nr O0 mv 

dS EL "n n pv 0 |g 

DE n = S 85 

ds 0 0 —pv ù pr sf) 
€v "y 0 -pr m 


for the Stokes vector 


Š= (86) 


where in the second equation we have replaced Q with L 
by adopting a coordinate system with U = 0 without loss 
of generality. Here e's are the emission coefficients, 7’s are 
absorption coefficients, py (the same as the f parameter 
in Equation (80)) is the coefficient for Faraday rotation 
and pr, (essentially the amplitude of h + ig in Equation 
(81)) is the coefficient for Faraday conversion. 

Apparent oscillations of L and V have been discovered 
in some bursts from rFRB 20201124A (Xu et al., 2022). 
These features may be interpreted as Faraday conversion 
or polarization-dependent absorption, which in any case 
demands a complex magnetized environment around the 
source (Li et al., 2022a; Xu et al., 2022). 


F. Plasma radiation mechanisms 


In classical electrodynamics, charged particles radiate 
when undergoing acceleration. Below we briefly discuss 
three well discussed radiation mechanisms involving elec- 
tron acceleration in electric fields, magnetic fields, and 
electromagnetic waves, respectively. 


1. Bremsstrahlung 


An electron in the Coulomb electric field of an ion 
would radiate through bound-bound (line emission), free- 


bound (recombination) and free-free (bremsstrahlung) 
processes. The opposite processes give respective absorp- 
tion processes of the photons. 

For a plasma in thermal equilibrium with temperature 
T, the plasma thermal bremsstrahlung (free-free) emis- 
sivity reads (Rybicki and Lightman, 1979) 


Ee — = 
Y” dVdtdv 
25566 2m 1/2 
= Z?Neni —hv/kpT = 
3m,c3 (xin) TG e 


= (6.8 x 107? erg cm ? s ! Hz!) 


x Znen; T Ve hvIkaT gg, (87) 


where c, kp, e, and m, are standard fundamental con- 
stants, T' is the gas temperature, n; is the number density 
of ions, Z is the atomic number of the ions, and gg is the 
Gaunt factor. The reason for the factor nen; is that the 
emissivity of each electron depends on the number den- 
sity of ions and the total emissivity is proportional to the 
number density of electrons. Since Zn; = n, is needed 
to maintain charge neutrality, n2 enters the problem, so 
a convenient emission measure 


D 
em = f n?dl 
0 


can be defined for a radio source, which may be related 
to the DM of the source through EML ~ DM?, where L 
is the characteristic size of the source. 

'The opposite process of bremsstrahlung, i.e. free-free 
absorption, is relevant to constrain the physical condition 
to allow the FRBs with the extremely high brightness 
temperatures to be observed. The absorption coefficient 
can be expressed as (Rybicki and Lightman, 1979) 


"2 4e6 2T 
" — 3mehe \3kgme 
= (3.7 x 108 em!) Zznan;T- 12 y73(1 — e^ PI kBT gg, 
(89) 


(88) 


1/2 
á 7) Z?neniv™’ (1 — e P" FBT) gg 


or, in the Rayleigh-Jeans regime 
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Integrating over distance, one gets the optical depth 
(Cordes and Lazio, 2002) 


D 
dim] a di 
0 


= (5.47 x 10 TZ 15 gaz (91) 


c em-6* 


An FRB is transparent only if 7# < 1 is satisfied in 
the emission region and also in the local environment 


surrounding the FRB source. 
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For a relativistically hot plasma, the emissivity and 
absorption coefficient should be multiplied by a correc- 
tion factor. The frequency-integrated correction factor 
is (1+ AT), where A = 44 x 101? K^! (Rybicki and 
Lightman, 1979). 


2. Cyclotron, synchrotron and curvature radiation mechanisms 


Electrons gyrate in magnetic fields and radiate. For 
non-relativistic electrons, the emitted cyclotron radiation 
spectrum is line-like, with the main power at the Larmor 
frequency wp and progressively lower powers at its higher 
harmonics. 

A relativistic electron with Lorentz factor ye radiates 
synchrotron radiation with a characteristic radiation fre- 
quency (Rybicki and Lightman, 1979) 


3 2 eB, 
WSR =~ alem a 
e 


(92) 


The power 2 for ye is due to the following three factors: 
(1) the relativistic mass is larger by a factor of ye; (2) 
the fraction of the orbital time with radiation beamed to- 
wards an observer is smaller by a factor of 2/4, due to the 
relativistic beaming effect; and (3) the observed timescale 
is shorter than the emission timescale by roughly a fac- 
tor of (1 — 8) ~ 1/(2y2), where 8 is the dimensionless 
speed of the electron. If synchrotron radiation is respon- 
sible for the FRB emission (e.g. within the framework of 
the synchrotron maser model), the required condition is 
12 By RE (360 G) VFRB,9- 

The relativistic beaming effect for synchrotron radi- 
ation is valid under the vacuum approximation. In 
a plasma, with the refraction index n, = ye < 1, 
the beaming angle 0, becomes 4/1 — n282 rather than 
V 1-— 8?. If n, deviates from unity much more than f, 
one has 6, ~ J/1— n2 = w,/w and synchrotron radia- 
tion is suppressed (Rybicki and Lightman, 1979). One 
may define the Razin frequency by equating 0, and 1/7, 
which gives 

WRazin = YVeWp- (93) 
Synchrotron radiation is suppressed when w < Wrazin- 
Matching the Razin frequency with GHz, the condition 


is y2ne ~ (1.2 x 101? cm~?) vag s. 


In a strong magnetic field environment such 
as the magnetosphere of a pulsar or a magne- 
tar, the synchrotron cooling timescale, tesp ~ 


Jemec? /[(4/3)7282cor(B?/8n)] ~ (8x 107? s) Brea}, 
is extremely short. As a result, charged particles stay at 
the lowest Landau level and essentially slide along mag- 
netic field lines in the local inertial (co-rotating) frame. 
Since the field lines are usually curved, particles will ra- 
diate when they accelerate in the curved trajectory. The 
characteristic frequency of such curvature radiation can 


be calculated by replacing the electron gyration radius 
in the synchrotron radiation formula by the curvature 
radius p of the field lines so that 


3 4c 
WCR = rey (94) 


The origin of 42 is similar to synchrotron radiation, ex- 
cept that there is no yYeMe suppression in gyration fre- 
quency for synchrotron radiation (the mass does not en- 
ter the problem since the curvature radius of the field line 
does not depend on mass). To match the GHz emission, 
the parameters should satisfy Y oie c 1.4 vrRB.. 


3. Compton and inverse Compton scattering 


An initially at-rest electron oscillates in electromag- 
netic waves and emits essentially isotropically at the 
same incident frequency if iw; < mec? with cross sec- 
tion equals the Thomson scattering cross section o = 
(87 /3)(e?/mec)? ~ 6.65 x 10775 cm?. When the elec- 
tromagnetic waves have an extremely large amplitude so 
that the electron reaches a relativistic speed (relevant 
to FRBs near the FRB generation site), the electron 
motion trajectory becomes complicated and cross sec- 
tion much enhanced (e.g. Yang and Zhang (2020), see 
Section VIILC for details). The existence of a strong 
background magnetic field further complicates the pic- 
ture (Beloborodov, 2021a; Qu et al., 2022). 

When an electron moves relativistically and interacts 
with electromagnetic waves with angular frequency wi, 
it would inverse Compton scatter the waves to a higher 
frequency 


ws ~ 72(1 — B cos6;)uj. (95) 


Such a process could be relevant to FRB radiation (Sec- 
tion V.B.3). 


IV. GENERAL CONSTRAINTS ON THE MODELS 


In order to interpret FRBs, a competent model needs 
to invoke a radiation mechanism model to address in- 
dividual burst properties (brightness temperature, po- 
larization properties, spectral down-drifting, radio effi- 
ciency, high-energy emission, etc) and a source model 
that accounts for the global properties of the bursts (en- 
ergetics, burst rate, luminosity/energy function, redshift 
distribution, host galaxy properties, etc.). Before dis- 
cussing these in detail in Sections V) and VI), one may 
place some generic, essentially model-independent con- 
straints on the models based on some basic observational 
facts and physical principles. 
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A. Burst duration (width) and engine size 


After correcting for the convolution effects from scat- 
tering and instrumental effects (§II.B), the intrinsic du- 
ration W; of an FRB defines a length scale 


Ri = cW; = (3 x 10’ cm) W3, (96) 


where W_s3 is the intrinsic duration in units of millisec- 
onds. The size of the FRB central engine Ro should 
satisfy Ro SZ R;. This is straightforward if the FRB 
emitter does not move with a relativistic speed. The rea- 
son is that if Ro > Ri, even if the emission region is lit 
up simultaneously everywhere, the duration of the event 
should be Ro/c > W; due to the light propagation delay 
between the front end and the rear end of the emission 
region with respect to the observer. 

If the FRB emitter is moving towards the observer with 
a relativistic speed (which is likely the case as discussed 
in 8IV.B and 8IV.F.1 below), the situation is more com- 
plicated but the conclusion of Ro < R; remains valid. 
Most generally, let us assume that the emitter travels 
with a bulk Lorentz factor I in a direction with an an- 
gle 0 with respect to the line of sight. In the lab frame, 
let us consider that the central engine sends off two light 
signals towards the relativistic emitter (the fastest causal 
connection is through propagation of photons), and the 
emitter promptly reacts to the two signals and send off 
two signals to the observer immediately after receiving 
the two central engine signals. Approximating the emit- 
ter as a point source and ignoring cosmic expansion, one 
can write the following relation between the three inter- 
vals (Zhang, 2018c): 


1— B cos0 
I-B 

where Ateng is the time interval for the engine to emit 
two signals; At, is the time interval for the relativistic 
emitter to receive the two signals from the engine and 
also the time interval for the emitter to send off two sig- 
nals; and At,ps is the time interval for the observer to 
detect the two signals. Here ( is the dimensionless speed 
of the emitter, 0 is the angle between the direction of mo- 
tion and line of sight, and the factor (1 — 8 cos0) (which 
~ 1/21? for 0 = 0) is a factor accounting for the propa- 
gation effect, and T = (1 — 82)-1/? is the Lorentz factor 
of the emitter. One can see that even though the emit- 
ter timescale is stretched due to its motion, the observed 
timescale Atop; still track the central engine timescale 
Ateng (tobs = teng for 0 = 0)9. As a result, W; can be 


Ateng = (1 — 8 cos 0)Ate = Atobs, (97) 


8 If the line of sight is outside the emission beam, (1 — B cos 0) ~ 
(1— cos0) which becomes > (1— 8), so Atobs becomes > Ateng- 
One can see a longer burst with a lower flux. Such “slow ra- 
dio bursts” may be detectable from Galactic magnetars or other 
FRB-emitting sources (Zhang, 2021). 


used to constrain the size of the central engine in any 
case. 

Equation (96) immediately suggests that the most 
compact, stellar-mass objects in the universe, i.e. a neu- 
tron star or a stellar-mass black hole, are the most likely 
candidates for FRB engine. Larger objects (e.g. white 
dwarfs, stars, and even supermassive black holes) have 
been invoked to interpret FRBs in some models, but these 
models must invoke contrived conditions to allow only a 
small enough region to power an FRB. 


B. Variability timescale and emission radius 


The rapid variability timescale, in particular the ~ 60 
ns timescale observed in rFRB 20200120E from the M81 
globular cluster, can be used to further constrain the 
emission radius of FRBs (Beniamini and Kumar, 2020; 
Lu et al., 2022). For an on-beam FRB (i.e. 0 ~ 0 for a 
point source, or 0 « 0; for a conical jet with an opening 
angle 6;), a natural variability timescale? is 


5t ~ PERB 
E CN 


(98) 


This timescale defines both the observed time for the 
emitter to travel to the emission radius R in the rising 
phase, and also the angular spreading timescale due to 
the propagation delay of a spherical jet front in the de- 
caying phase. In principle, if one is allowed to arbitrarily 
increase the Lorentz factor of the emitter, any small ót 
can be reproduced for any R. So Eq.(98) alone is not con- 
straining. Interesting constraints can be posed when the 
duration of the burst W is considered. For certain mod- 
els, for example, the synchrotron maser model invoking 
the external shock (Metzger et al., 2019), the emission 
radius can be estimated as Rpggp ~ I?cW. This immedi- 
ately suggests that ôt cannot be significantly shorter than 
W. The 60-ns variability from the ms-duration bursts 
of the M81 globular cluster FRB (Nimmo et al., 2022) 
therefore disfavors the external shock model of FRBs 
(Lu et al., 2022). The synchotron maser internal shock 
model (Beloborodov, 2020) is still allowed. However, it 
suffers from other drawbacks. For example, the frequency 
down-drifting feature, which the external shock model in- 
terprets as the shock propagating to progressively larger 
radii (Metzger et al., 2019), is no longer straightforwardly 
interpreted within the internal shock models. On the 
other hand, rapid variability of FRBs is not a challenge 
to the magnetospheric models, as a 0.4-nanosecond pulse 


9? Scintillation (see 8VIILA for more discussion) can introduce 
modulations in shorter timescales but with small amplitudes. 
Distinct pulses in an FRB lightcurve should be intrinsically re- 
lated to the size of the source or emission region. 
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has been observed from the magnetosphere of the 33- 
ms Crab pulsar (Hankins and Eilek, 2007) (even though 
FRBs are more energetic than nano-shots). 


C. Periodicity 


The special source, FRB 20191221A, was detected to 
have a 0.2168 s periodic separation during a 3 s duration 
(Chime/Frb Collaboration et al, 2022). Since known 
sources of a sub-second period are all rotating neutron 
stars (pulsars), this source offers a definite clue that at 
least some FRBs originate from pulsar-like objects. Fur- 
ther arguments can be made that the FRB radiation re- 
gion (at least for this source) is the magnetosphere of 
an underlying pulsar or magnetar (Beniamini and Ku- 
mar, 2022; Chime/Frb Collaboration et al., 2022). This 
is because models invoking emission regions outside the 
magnetosphere do not have well-defined geometric win- 
dows to maintain a strict periodic window. 

The lack of periodicity from active repeaters such as 
rFRB 20121102A (Hewitt et aL, 2022; Li et al., 2021b; 
Zhang et al., 2018a) and rFRB 20201124A (Niu et al., 
2022; Xu et al., 2022), on the other hand, places less 
constraints on the models. Katz (2020a) argued that this 
suggests a black hole rather than a neutron star origin of 
repeating FRBs. This argument is not strong because un- 
like pulsar emission, FRB radiation pressure is so strong 
that the magnetospheric structure is likely significantly 
distorted so that a well-defined magnetospheric window 
(the conventional open field line region) likely does not 
exist and it is entirely possible that an FRB emitting 
neutron star emit bursts at random phases. The radio 
bursts from the magnetar SGR J19354-2154 seem to be 
emitted from a much wider phase window than the nar- 
row window for pulsed emission (Zhu et al., 2022). With 
the burst data alone, it appears that the source does not 
have a strict periodicity, even though the magnetar has 
a strict 3.24-s period. 

So far, only rFRB 20180916B has been confirmed to 
possess a long-term 16-day periodicity (Pastor-Marazuela 
et al, 2021; Pleunis et al, 2021b; The CHIME/FRB 
Collaboration et al., 2020). Its origin is subject to de- 
bate. The most natural interpretation would be to at- 
tribute this to the orbital period of a binary system, with 
the emission from the FRB emitter only reaches the ob- 
server in a particular phase window (Ioka and Zhang, 
2020; Lyutikov et aL, 2020; Wada et al., 2021). Other 
interpretations to the 16-d period of rFRB 20180916B 
include magnetar precession (Levin et al., 2020; Yang 
and Zou, 2020), slowly rotating magnetars (Beniamini 
et al., 2020), and even precession of a black hole ac- 
cretion disk (Katz, 2022a). None of these models were 
theoretically predicted before the discovery of the rFRB 
20180916B periodicity. So, it would be uncomfortable, at 
least to me, if such long-term periodicty is a common fea- 


ture of active repeaters because that would require such 
periodicity being at the heart of FRB generation mech- 
anisms (Zhang, 2020d). It is now clearer that such a 
long-term periodicity is not commonly observed among 
active repeaters (the case of rFRB 20121102A is to be 
confirmed, see §II.D). Whatever mechanism that is oper- 
ating in rFRB 20180916B likely applies in rare cases and 
is probably due to a chance coincidence. 


D. Energetics, radio emission efficiency, and beaming 


The derived isotropic energies of individual bursts and 
the energy-dependent burst rates for repeaters can be 
used to place interesting constraints on the average lu- 
minosity and total energy budget of the underlying FRB 
source, which may be used to constrain FRB source mod- 
els. 

For one-off FRBs, the true peak luminosity and energy 
of the burst are 


Ly — Lsssafotir 3 
E- Eisen. s (99) 


where Ly iso (Eq.(16)) and Eiso (Eq.(17)) are the mea- 
sured isotropic radio peak luminosity and energy directly 
from observations, 7, is the radio emission efficiency, and 
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fo (100) 
is the beaming factor of an individual burst, with dQ be- 
ing the solid angle of the burst. Note that fy reduces 
and 5.! increases the energy budget of the source so 
that the effects of the two factors tend to cancel out each 
other. Both factors are not well constrained from obser- 
vations. The X-ray burst associated with FRB 200428 
was more than 10+ more energetic than the radio burst 
itself (Li et al., 2021a; Mereghetti et al., 2020; Ridnaia 
et al., 2021), so for this particular event, the upper limit 
of n, is ~ (107^ — 107°). Various X-ray flux upper limits 
for extragalactic FRBs places a lower limit on 7,, which 
is of this order or even smaller (Laha et al., 2022a,b; Piro 
et al., 2021). 

For repeating sources, one should consider the aver- 
age energy-dependent bursting rate dN/dtdE;,. during 
the active phase and the observational duty cycle of the 
active phase Ç (e.g. for the rFRB 20121102A observing 
campaign with FAST (Li et al., 2021b), the observational 
duty cycle is about 60 hours out of 47 days). One should 
also introduce a global beaming factor 


jum. (101) 


which can be larger than f, of the individual bursts. This 
is because the global emission beam can have a larger 
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FIG. 8 A cartoon picture of the beaming factor of individual 
bursts (with a solid angle 62) and global beaming (with a 
solid angle of AQ. A fan-beam from a magnetospheric rotator 
is illustrated as an example for global beaming, but a more 
general geometry is possible. 


solid angle AQ, inside which each burst could have a nar- 
rower beam (See Fig.8). The average source luminosity 
that is used to make FRBs would be 


L i EN) moo (Fynz t i)dEso, (102) 
src — Bao dtd Eiso iso br iso; 


where Fisom and Ei; are the minimum and maxi- 
mum isotropic FRB energy from the source. The rea- 
son that P, rather than f is adopted is the following. 
Even though each burst has a beaming factor fj, there 
are altogether AQ/dQ such bursts on average most of 
which are not detected but are added to the energy bud- 
get of the source. The final beaming factor is therefore 
(AQ/6Q) fy = Fy. 

Assume that the repeater source has a lifetime 7 and 
the activity level remains unchanged during lifetime. The 
total energy budget of the source may be estimated as 


T pEiso.M dN 
Esc = Fore Eiso F, e dE iso t 
I J, (a) (Fin ! C )dEssodt, 


(103) 

Interesting constraints on the energy budget of re- 
peaters have been made. For rFRB 20121102A, Li 
et al. (2021b) reported that the 1652 bursts detected 
in ~ 60 hours during a 47-day observational campaign. 
The total emitted radio energy (corresponding to the 
integral in Eq.(103) without the (Fyn-'¢~! factor) is 
~ 3.4 x 10! erg. Considering ¢ = 60/(47- 24) = 
0.053, nr = 10-4, 4, and F, = 0.1F,,-1, one can de- 
rive that the total source energy used to make FRBs 
is Esc = (6.4 x 1099 erg) Fs, 10,.1,(6/0.053)-!. This 
is already ~ 496 of the total dipolar magnetic energy 
(Ep ~ (1/6) B2 R? ~ (1.7 x 10** erg) B2, Rẹ) of a magne- 
tar. Since rFRB 20121102A already existed more than 
a decade, this observation already poses a significant en- 
ergy budget issue for the magnetar model, unless 7, is 
larger or Fp is smaller. The magnetospheric models sat- 


isfy these constraints!?, but the synchrotron maser shock 
model is already severely constrained by the data (see 
§VI.A for details). Note that before the FAST observa- 
tions, Margalit et al. (2020b) already posted a tight en- 
ergy budget constraint on rFRB 20121102A based on the 
previous data within the framework of the magnetar syn- 
chrotron maser model. The many more bursts detected 
with FAST (Li et al., 2021b) only tightened the con- 
straints for the source. Even more stringent constraints 
on the magnetar synchrotron maser model have been es- 
tablished for another active repeater rFRB 20201124A 
(Xu et al., 2022; Zhang et al., 2022). 


E. Brightness temperature and coherent radiation 


Astronomical objects emit four levels of electromag- 
netic radiation with increasing complexity (Table II): 
blackbody radiation, thermal radiation, incoherent non- 
thermal radiation, and coherent non-thermal radiation. 
Blackbody and thermal radiations require that the emit- 
ting particles are in thermal equilibrium (defined by the 
gas temperature T), with the blackbody radiation hav- 
ing an additional requirement that photons have a large 
enough optical depth to reach thermal equilibrium as well 
(Rybicki and Lightman, 1979). Thermal radiation in- 
cludes blackbody radiation, but also allows photons not 
to achieve thermal equilibrium. One example is thermal 
bremsstrahlung which has a different spectral shape from 
blackbody with the cutoff energy defined by the gas tem- 
perature T'. If particles are accelerated to deviate from 
thermal equilibrium, say, in shocks or magnetic reconnec- 
tion regions, the radiation becomes non-thermal. As non- 
thermal particles typically have a (one segment or multi- 
segment) power-law distribution, non-thermal radiation 
spectra are typically (broken) power-laws. In particu- 
lar, in the low-frequency regime, non-thermal radiation 
is usually subject to self-absorption from the opposite 
process of the emission mechanism if particle radiation is 
incoherent. 

Coherent non-thermal radiation may be defined in dif- 
ferent ways, but the most straightforward way is through 
its ability of overcoming the self-absorption limit. The 
self-absorption-defined specific luminosity limit at a par- 
ticular frequency is the blackbody specific luminosity at 
that frequency for a gas with the maximum temperature 
and source size allowed by the emitter. For electrons 
with a characteristic Lorentz factor ye in the comoving 
frame, the comoving-frame effective temperature would 


be kT’ ~ yemec?. For a synchrotron radiation source, 


10 Because of unidentified coherent mechanism of FRBs, the ra- 
dio efficiency in the magnetospheric models cannot be predicted. 
However, the radio emission efficiency of radio pulsars can range 
from 1078 to close to unity (Szary et al., 2014). 
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Ye = Max(Ym, Ya) is usually the maximum of the follow- 
ing two: the minimum injection Lorentz factor 7, and 
the corresponding Lorentz factor for self-absorption Ya 
(e.g. Kumar and Zhang, 2015). As a result, 


kT’ = y,mec? = max(Ym, Ya)Mec?, (104) 


or 


T' = (5.9 x 10! K)ye2 (105) 


defines the maximum incoherent brightness tempera- 
ture in the comoving frame, where ye ~ 100 has been 
adopted. For radio galaxies, Kellermann and Pauliny- 
Toth (1969) showed that the observations had a max- 
imum brightness temperature of about Tmax ~ 107? 
K, which corresponds to a typical electron Lorentz factor 
ye ~ (102—103). They argued that this limit is physically 
related to the requirement that the second-order Comp- 
ton scattering power does not exceed that of the first- 
order (synchrotron self-Compton) through self-regulated 
Compton cooling. For systems like GRBs or blazar jets, 
Je is related to the bulk Lorentz factor of the jet, which 
directly defines the internal energy density (and hence, 
the effective temperature) in the emission region. 

For a relativistic emitter beaming towards earth, the 
allowed maximum radio specific flux at a frequency is 
larger than the comoving value by a factor of D for an 
extended source or D? for a point source (Zhang, 2018c), 
where 


1 
Dz——————er 


T(1 — 8cos0) Gi 


is the Doppler factor, IT and 0 carry the same meaning 
defined earlier, and the final approximation applies to 
the regime 6 < 1/T. Compared with the observation- 
ally defined brightness temperature (Eq.(18)), one may 
perform either of the following two approaches. One is 
to derive brightness temperature in the comoving frame 
(T7) and compare with the maximum 7"; the other is to 
derive the maximum allowed T in the observer frame and 
compare with observationally-defined T, (Eq.(18)). We 
adopt the more straightforward latter approach and de- 
rive the condition that coherence is required by the data 
if 


Ty > Dyemec?/k ~ (5.9 x 1015 K)T 25, 5. (107) 


This result is consistent with Lyubarsky (2021) who 
adopted the opposite approach. Since FRBs have the 
observed T, much greater than this value, their radiation 
mechanisms must be coherent. We will discuss various 
coherent mechanisms in SV. 


F. Attenuation processes 


In order to have high-T, radio pulses detectable from 
Earth, the radio waves must overcome various absorption 


TABLE II Astrophysical radiation mechanisms 
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Mechanisms Particles Photons Examples 
blackbody thermal equilibrium thermal equilibrium CMB, stars 
thermal thermal equilibrium may/may not in thermal equilibrium disks, intracluster medium 


non-thermal 
non-thermal 


incoherent non-thermal 
coherent non-thermal 


subject to self-absorption limit 
not subject to self-absorption limit 


SNRs, GRBs, blazars 
radio pulsars, FRBs 


or scattering processes along the propagation paths. The 
three important processes to attenuate the radio emis- 
sion flux are induced Compton scattering, free-free ab- 
sorption, and synchrotron absorption, which we discuss 
below in turn. 


1. Induced Compton scattering and Lorentz factor lower limit 


With the existence of free electrons, photons with a 
particular frequency can be scattered out of the state 
and other photons with different frequencies can be scat- 
tered into the state. The Thomson scattering optical 
depth can be estimated as vp ~ n,oqR, where n, is 
the electron number density, or = (87/3)(e?/mec?)? ~ 
(6.65 x 107?5 cm?) is the Thomson cross section, and 
Ris the size of the emission region. This optical depth 
is relevant for scattering of high-frequency photons (for 
very high-energy photons, the Klein-Nishina correction is 
needed), but in the low-frequency regime, scattering can 
be enhanced significantly by induced Compton scatter- 
ing if Tj is high enough (Kompaneets, 1957; Lyubarsky, 
2008; Thompson et al., 1994; Wilson and Rees, 1978). 
'The essence of this mechanism can be summarized as 
follows (Wilson and Rees, 1978). Consider two photon 
states (not electron states) a and b (defined by both the 
energies and directions of the photons) with photon oc- 
cupation numbers na and nj, respectively. The sponta- 
neous change in n4 because of scattering from a to b is 
dn, /dt x —n4. However, since photons are bosons that 
satisfy the Bose-Einstein statistics, the existence of pho- 
tons at b actually boost the scattering rate from a to b by 
a factor of (ny +1), i.e. dna /dt = —na(1 +n). Similarly, 
the scattering rate from b to a is dna /dt x (na+1)ny. The 
net change at level a is dn, /dt « [(na 4- 1)ny — n4 (1 +na)] 
which essentially cancels out but leaves a small term re- 
lated to the recoil frequency shift due to Compton scat- 
tering, ie. Av/v = (hv/mec?)(1 — cos0), where 0 is the 
angle between the directions of a and b. It is found that 
for cold electrons without bulk motion, induced Compton 
scattering becomes important when (kgTy/m,c?)? > 1 
(where Q is the solid angle of the uniform beam) (Wilson 
and Rees, 1978)!!. As a result, the optical depth due to 


11 The factor (kgTp/mec?)Q? is the product of the photon oc- 
cupation number kpTy/hv and the fractional change of energy 


induced Compton scattering is enhanced with respect to 
Thomson scattering by the same factor, i.e.!? 


= 


mec? 


TO o) Tr œ (6.4 x 1074)0?T; 367r. (108) 
The detailed expression depends on the explicit prob- 
lems one is addressing. For example, if one considers 
the induced Compton scattering constraint in an emit- 
ting source, the expression can be written as (Lyubarsky, 
2008) 


(109) 
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where S95 is the observed specific flux of the FRB, v is 
the FRB frequency, Dy, is the luminosity distance of the 
source, rg is the radius of the launching point, and Z is 
an integral that has the dimension of ro/c and carries the 
information of €. For another example, if one considers 
the induced Compton scattering by a medium as an FRB 
from a separate source passes through it, the expression 
becomes (Ioka and Zhang, 2020) 


" 3ov LynecAt 
~ 822? rmen? ' 


TC (110) 
where r is the distance between the FRB source and the 
scatterer. Note that the above discussion applies to an 
unmagnetized plasma with wg « wp. In a highly magne- 
tized environment such as the magnetosphere of a neu- 
tron star, charged particles are confined in strong mag- 
netic fields, so that the particles required to have the right 
directions and energies for induced Compton scattering 
are not available. As a result, there is no need to consider 
the induced Compton scattering constraint in the emis- 
sion region if FRBs are emitted from the magnetosphere 
of a central engine. 

When the emitter is moving relativistically with a bulk 
Lorentz factor I, the induced Compton scattering opti- 
cal depth drops significantly (Lyubarsky, 2008). From 
Eq.(108), noticing Q c L^? and T} = T/D ~ Ty/T, one 
gets 


kpTy Tr 


— = 111 
Mec? TS’ (111) 


TQ ^ 


~ hv /mec?. 
12 A coefficient 3/87? is added with precise calculations (e.g. W. 
Lu, 2021, unpublished notes). 


which is significantly smaller than the case without bulk 
motion. 

The induced Compton scattering optical depth also 
drops if the electron gas is relativistically hot. Lu and 
Kumar (2018) suggested that for a narrow Gaussian-like 
spectrum with a characteristic electron energy ye, the 
approximated optical depth is 


kT, TT 
mee 7 


TQ ^ (112) 
For à power law photon spectrum, the results depend 
on the spectral index p (convention I, œ v?) but the 
suppression factor is shallower than »y; ?. 

Some FRB emission models invoke relativistic shocks 
as the emission site (e.g. Beloborodov, 2017, 2020; 
Lyubarsky, 2014; Metzger et al., 2019; Plotnikov and 
Sironi, 2019). The emission region would be also rela- 
tivistically hot in the comoving frame. Combining Equa- 
tions (108), (111) and (112) and noticing 9? is already 
included in the 1/T? suppression factor, one may derive 
that the condition of ro < 10 is? 


Tye > (5.8 x 104) Ty er > 


(113) 
Since ye x I is generally expected!^, one may place a 
lower limit of I as 
1/10, 1/10 
D 0 uS (114) 
where ye = el has been assumed. Similar constraints 
have been derived by Lyubarsky (2008); Murase et al. 
(2016). Note that within the relativistic shock models, 
a Lorentz factor of this order is also required by the du- 
ration and variability constraint (Eq.(98)), so induced 
Compton scattering constraint is usually satisfied in the 
shock model without introducing an additional condition 
(e.g. Beloborodov, 2020; Metzger et al., 2019). 


2. Free-free absorption 


Radio emission can be also attenuated via free-free 
absorption, the inverse process of free-free emission or 
bremsstrahlung. The importance of free-free absorption 
for FRBs has been discussed by various authors (Kumar 
et al., 2017; Kundu and Zhang, 2021; Luan and Goldre- 
ich, 2014; Metzger et al., 2017; Murase et al., 2016; Yang 
and Zhang, 2017). 


13 Notice that induced Compton scattering mainly modifies the 
shape of the spectrum rather than exponentially attenuate pho- 
ton flux. As a result, a larger optical depth than unity, e.g. 
Tc = 10, is adopted as the transition point where the effect be- 
comes important. 

14 'This is straightforwardly expected for external shocks. For in- 
ternal shocks, ye is more related to the relative Lorentz factor 
between the shocks, which may also scale with r. 
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The free-free absorption coefficient (Eq.(90) together 
with the relativistic correction factor (1+ AT) (8IILF.1) 
can be used to estimate the optical depth against free- 
free absorption. An FRB is transparent if the optical 
depth is below unity. 

Free-free absorption is important when the density of 
the emitter or environment medium is high. Therefore, 
the free-free absorption constraint was adopted (Luan 
and Goldreich, 2014) to disfavor an early FRB model 
invoking flaring stars (Loeb et al., 2014). For repeating 
FRB models invoking a young magnetar born from a su- 
pernova explosion, free-free absorption was used to place 
a lower limit on the age of the supernova remnant be- 
fore which the remnant shell is too dense to allow FRBs 
to escape freely (Metzger et al., 2017; Yang and Zhang, 
2017) (see §VII for details). For FRB systems invoking 
relativisic shocks, either as the site of FRB emission or 
as a screen in front of FRB produced at an inner radius, 
free-free absorption in the hot shocked plasma could be 
important if the total kinetic energy exceeds ~ 10“4 erg, 
which may account for the frequency down-drifting fea- 
ture observed in some FRBs (Kundu and Zhang, 2021). 


3. External synchrotron absorption 


For active repeaters surrounded by a persistent radio 
source (PRS) (Chatterjee et al., 2017; Niu et al., 2021), 
coherent FRB emission needs to pass through the PRS, 
which is likely powered by synchrotron radiation. Under 
certain conditions, FRBs could be absorbed by the PRS 
via synchrotron absorption and the PRS source could 
be subsequently heated up by the absorbed FRBs (Yang 
et al., 2016). 

Assuming that the nebula electrons have an initial dif- 
ferential number density spectrum N (ye, 0) = Kyz?, one 
can estimate the synchrotron optical depth as 


?KR 1 CE 
€ V 
Tu SR = p (2) f(p), 


4mec vg Nvp 
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where vg = eB/(2xm.,c), R is the radius of electron 
acceleration region, and f(p) is a function of order unity. 
Solving T, sr = 1, one can derive synchrotron absorption 
frequency (Yang et al., 2016) 


(116) 
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The spectrum of the nebula needs to be solved numer- 
ically by including electron injection, synchrotron cool- 
ing, as well as heating by FRBs, which would give rise to 
complicated spectra for both electrons and photons. The 
predicted spectra (Yang et al., 2016) turn out to share 
the general shape of the later observed PRS spectrum of 
rFRB 20121102A (Chatterjee et al., 2017; Marcote et al., 
2017), as shown in Li et al. (2020a). The small nebula 


size and not too high a synchrotron self-absorption fre- 
quency constrain the parameter space for such models in 
general (Metzger et al., 2017). 


G. Ordered magnetic fields and strengths 


The fact that FRB emission is linearly polarized with 
a high polarization degree poses a generic constraint, 
namely, there must exist ordered magnetic fields in the 
FRB emission region. Indeed, current leading models to 
interpret FRBs invoke either magnetospheres of magne- 
tized central engines or relativistic shocks with ordered 
magnetic fields. 

Further constraints on the strength of magnetic fields 
have been discussed in the literature (Kumar et al., 2017; 
Lyutikov, 2017). The argument is that the electromag- 
netic wave energy density in the emission region should 
not exceed the magnetic energy density of the emitter 
in the same region before the FRB is emitted. Such a 
constraint can be placed if the FRB emission originates 
from dissipation of magnetic fields, or the magnetic field 
in the emission region confines the generated FRB emis- 
sion. Note that such a condition in general is not al- 
ways necessary for producing intense electromagnetic ra- 
diation. For example, the fireball model for GRBs does 
not require to abide by such a condition, with the elec- 
tromagnetic energy of radiation generated from the ther- 
mal energy or the dissipated kinetic energy in the fireball 
(Zhang, 2018c). In the case of coherent radiation, on the 
other hand, many models require that ordered B fields 
should remain ordered during the emission processes. As 
a result, such a condition is quite relevant. 

The electromagnetic wave energy density, indepen- 
dent of the emission frequency, may be estimated as 
Liso/ (4r Rĉgpgc), where Rpgnp is the radius where FRB 
emission is radiated. The condition 


Liso B? 
——5 — X = 117 
AnRiagc 8r (aday 
gives 
2Lio 1 = 
B> ~ (8.2 x 101? G)Liss 45 pig. (118) 
HrnB 


The key is how to estimate Rprp. If one assumes 
Rrra = cW; = (3 x 10’ cm) (Wms), one obtains 
B > (2.7 x 108 G)Liso.42(Wms)~', which leads to the 
conclusion that the emission region has to be within 
the magnetosphere of a neutron star (Lyutikov, 2017). 
This argument, however, is flawed, since Rprp cannot 
be always simply estimated as cW;. If the emitter is 
moving relativistically with a bulk Lorentz factor T, as 
is envisaged in the synchrotron maser models, one has 
Rrrp = I?cW; = (3 x 1018 em)T2 W_3. The B-field 
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constraint becomes 
QLiso 1 -25p-— 
B>4| c Paw, = (2.7 x 10° G) Lis, 2 D3?W 1. 

(119) 

Note that the magnetic field strength at the light cylin- 

der of a magnetar is By, ~ B,(cP/2rR,) ? = (9.2 x 

10? G)B, iP ? R2. So this estimate allows the emis- 

sion region to be outside of a neutron star magnetosphere. 


H. Afterglow 


A genetic constraint can be placed on the brightness 
of the multi-wavelength afterglows of FRBs. Afterglow 
observations for GRBs have been essential in identifying 
their multi-wavelength counterparts and host galaxies as 
well as measuring their redshifts. In the case of FRBs, 
the isotropic energy is typically more than 10 orders of 
magnitude smaller than GRBs (Exo,rnp ^ 10?? erg vs. 
Eio,anB ~ 10?? erg). The expected FRB afterglow emis- 
sion is expected to be much fainter (Yi et al., 2014). One 
possible way of enhancing afterglow emission is to as- 
sume that the FRB radiative efficiency 7, is very low 
so that the afterglow kinetic energy can be boosted by 
a factor of 5. !. According to the standard GRB after- 
glow model (Mészáros and Rees, 1997; Sari et al., 1998; 
Zhang, 2018c), the characteristic synchrotron frequency 
of injected minimum-energy electrons and the peak syn- 
chrotron specific flux for a relativistic jet being deceler- 
ated by a constant-density medium read 


Vm = (33 x 105 Hz 1 z) Pt; Peg, 
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where the blastwave kinetic energy in normalized to 1044 
erg (which assumes 7, = 10-9 for Eprg = 10°8 erg), 
€e and eg are shock equipartition parameters for elec- 
trons and magnetic fields, respectively, p is the power 
law index of the injected electrons, n is the medium den- 
sity, tq is the observing time in units of day, and Dz 2 
is the luminosity distance of the source in units of 1028 
cm. One can see that the afterglow emission peaks in 
the radio band and is extremely faint. Detailed calcula- 
tions (Yi et al., 2014) suggest that a detection is possible 
only if the source is extremely nearby and the FRB is 
extremely energetic (i.e. the radio efficiency is very low), 
e.g. E = Epgp/n- = 10°’ Egnp 40r, 7 erg. For a rela- 
tivistic, mildly magnetized jet, the reverse shock emission 
could be brighter than the forward shock emission, which 
would ease the detection of the afterglow (Yi et al., 2014). 

No confirmed FRB afterglow has been detected so 
far (even for the Galactic FRB 200428). 'This is con- 
sistent with the theory and suggests that 7, is not ex- 


tremely low. It is worth noting that in the synchrotron 
maser model invoking external shocks (Metzger et al., 
2019), the multi-wavelength counterpart associated with 
the FRB could be regarded as its own “afterglow”, even 
though the electron energy distribution is assumed to be 
thermal rather than a power law. No Fermi acceleration 
of particles is envisaged, which could be a problem the- 
oretically. The two hard spikes observed in the X-ray 
counterpart (Li et al., 2021a; Mereghetti et al., 2020) of 
FRB 200428 (Bochenek et al., 2021; CHIME/FRB Col- 
laboration et al., 2020) can be interpreted within this 
model as the external shock emission (Margalit et al., 
2020a), even though it is more naturally interpreted as 
emission within the magnetar magnetosphere (Lu et al. 
(2020), Yang and Zhang (2021), Ioka (2020)). 


V. COHERENT RADIATION MECHANISMS 


Coherent radiation mechanisms invoke fundamental 
plasma physics, which could be shared among different 
source models. For example, coherent curvature radi- 
ation by bunches has been discussed in many different 
contexts involving magnetospheres, such as radio emis- 
sion from the inner magnetospheres of pulsars and mag- 
netars (Katz, 2014; Kumar et al., 2017; Ruderman and 
Sutherland, 1975; Yang and Zhang, 2018), from ejected 
magnetospheres from “blitzars” (Falcke and Rezzolla, 
2014; Zhang, 2014), from kinetic-energy “combed” mag- 
netospheres (Zhang, 2017), from magnetopsheres during 
asteroid-NS collisions (Dai, 2020; Dai et al., 2016; Geng 
and Huang, 2015), as well as from the global magne- 
tospheres formed by merging charged objects (Zhang, 
2016). The synchrotron maser mechanism in relativsitic 
shocks, on the other hand, has been invoked in the 
magnetar internal (Beloborodov, 2017, 2020) or exter- 
nal (Lyubarsky, 2014; Metzger et al., 2019) shock mod- 
els, shocks from low-B compact objects (Long and Pe’er, 
2018; Waxman, 2017), and even black hole accreting sys- 
tems (Sridhar et al., 2021a). Therefore, it is reasonable to 
detach radiation models from source models and discuss 
the general physics behind each radiation model. This is 
the task of this section. 


A. Coherent radio emission overview 


Following the discussion in SIV.E, we can summarize 
two fundamental properties of a coherent radiation mech- 
anism: (1) the observed luminosity, Lops, exceeds the 
sum of the emitted power P, for individual particles, i.e. 
Lobs > NeP-, where Ne is the total number of electrons; 
and (2) the observed luminosity is not subject to self- 
absorption, so that Equation (107) is satisfied. 

There are several ways to classify coherent radiation 
mechanisms. Based on differences in general physics, one 
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may classify the mechanisms in the following three types 
(e.g. Melrose, 1978). Each mechanism has its emission 
properties and back-reaction mechanisms. 


e Coherent emission by bunches (or the “antenna” 
mechanism): In this mechanism, emitting particles 
are physically clustered in six-dimensional phase 
space, i.e. in both 3-D position space and 3-D mo- 
mentum space. This is how coherent emission is 
emitted from antennae in radio stations. Within 
this mechanism, microscopic particles (e.g. elec- 
trons) are physically bunched together to radiation 
as a global particle with a total charge Ne pe, where 
Nep» is the number of charges in each bunch, typ- 
ically distributed within a unit volume defined by 
the wavelength of the radio waves (Nep ~ ney? X3, 
where ne is the charge number density, A is the 
wavelength, and y is the bulk Lorentz factor of 
the bunch). The emission power of the bunch, de- 
pending on the degree of coherence, can reach a 
maximum of N2,P. (e.g. Yang and Zhang, 2018). 
The total luminosity of the system would be ~ 
N2,NyP., where Ny ~ Ne/Ne,» is the number of 

bunches in the emission region. The back-reaction 

effects of such bunched emission are two folds: due 
to internal Coulomb repulsion, bunches tend to dis- 
perse in space. Radiation reaction may also make 
the particles disperse in the momentum space (Mel- 
rose, 1978). 


e Hydrodynamic instabilities (or “plasma masers"): 
In this mechanism, some oscillation modes in a 
plasma exponentially grow with time, with macro- 
scopic particles clustering in the momentum space. 
The MHD waves eventually escape in the form of 
electromagnetic waves in the radio band. The back- 
reaction effect is that as the mode grows, dispersion 
in the momentum space occurs and the instability 
would suppress itself. 


e Kinetic instabilities (or ^vacuum masers"): In this 
mechanism, electromagnetic waves detached from 
the plasma fluid would undergo negative absorp- 
tion in an energy-population-inverted medium so 
that the amplitude of emission grows with distance, 
reaching a high brightness temperature. The effect 
of back-reaction is that masers tend to reduce pop- 
ulation inversion so that the instability also sup- 
presses itself. 


Only a few types of objects are observed to emit coher- 
ent radio emission, e.g. Sun, Jupiter, astronomical maser 
sources, pulsars, and FRBs. The mechanisms operating 
in different types of objects can achieve different degrees 
of coherence (i.e. different values of Ty). Melrose (2017) 
reviewed the mechanisms of coherent emission in different 
types of objects and suggested that they have different 


origins: (1) Plasma emission at the plasma frequency 
wp, Which invokes Langmuir plasma waves (longitudi- 
nal oscillations) through a streaming instability as the 
trigger mechanism, likely applies to solar radio bursts; 
(2) Electron cyclotron maser emission at the cyclotron 
frequency wp, which invokes a cyclotron plasma insta- 
bility, likely applies to Jupiter and Earth Aurora; (3) 
pulsar coherent emission must have a different mech- 
anism, which has at least four possibilities: curvature 
emission by bunches, linear acceleration emission, rela- 
tivistic plasma emission, and anomalous Doppler emis- 
sion. However, all four mechanisms encounter difficulties 
and the pulsar coherent mechanism remains an enigma 
after more than half a century of study. 

The prospect of understanding FRB coherent emis- 
sion is not bright, either, since it involves more extreme 
processes to produce coherent emission. In any case, 
many mechanisms have been discussed in the literature 
and some have been briefly reviewed in Zhang (2020b), 
Lyubarsky (2021) and Xiao et al. (2021). In the fol- 
lowing, we will present a critical review on various FRB 
coherent radiation models, which are generally grouped 
into two types based on the emission region: those in- 
volving magnetospheres (also called “close-in” or “pulsar- 
like” models) and those invoking relativistic shocks far 
outside of the magnetopsheres (also called “far-away” or 
“GRB-like” models). 


B. Magnetospheric models 


All the pulsar-like mechanisms proposed for FRBs, as 
expected, have been proposed to interpret pulsar radio 
emission. In the following, we will discuss these mecha- 
nisms in turn, each with a brief introduction within the 
pulsar context, and then with a critical evaluation on 
its motivations and issues to account for FRB emission. 
Some pulsar mechanisms that have not been reinvented 
for FRBs are discussed in the end. 


1. Pulsar magnetosphere basics 


Before going over detailed pulsar-like models, it is in- 
formative to review the basic physics of pulsar magneto- 
spheres. 

Consider a pulsar that carries a plasma-loaded mag- 
netic field and rotates with an angular velocity Q. Let us 
make two idealized assumptions: 1. The plasma has in- 
finite conductivity so that the net force received by each 
particle is zero, Le. e(E + (1/c)((( x F) x B)) = 0 (the 
ideal MHD condition, which is also the force-free condi- 
tion as explained below); 2. The rotating magnetosphere 
is in a steady state so that the 0/Ot terms in Maxwell 
equations are zero (strictly applies to a uniformly rotat- 
ing, ax fip = 0 rotator, where fig is the direction of the 
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magnetic axis, which is either parallel or anti-parallel to 
the direction of the spin axis Q). From Maxwell equa- 
tions and with some basic vector calculus, one can derive 
that everywhere in the magnetosphere within the light 
cylinder radius 


P 
Rio = = = — = (4.8 x 10° cm)(P/1 s), 


122 
Q Qn (122) 


the net charge density as observed in the inertial frame of 
an observer who watches the star rotates is the Goldreich- 
Julian density (Goldreich and Julian, 1969) 


| ÖB 1 GB 
Pe = PGJ = ie. (Rey re? 


Cc 


where B is the local magnetic field at a location in the 
magnetosphere, and for a dipolar field, its strength falls 
with radius r as B ~ B,(r/R)~3, where R is neutron star 
radius and B, is the surface magnetic field strength. The 
last approximation applies to the region well within the 
light cylinder. This corresponds to a net charge number 
density 


nas = pa3/e ~ (6.9 x 101? cm?) Bj9P^!. — (124) 
By definition, with such a density there is no E com- 
ponent parallel to the local B vector (ie. Ej — 0), and 
(E x B) drift velocity is just the velocity v to allow parti- 
cles to be frozen in the magnetic fields and co-rotate with 
the star, Le. (Ex B)/B? = i/c. The local current density 
can be simply denoted as j = pet’, so the ideal MHD con- 
dition E + (1/c)(% x B) = 0 condition can be also trans- 
lated to the “force-free” condition peE+(1/c)(jx B) = 0. 
For an oblique rotator (Q x fig # 0), the 0/dt = 0 as- 
sumption is no longer satisfied, but particle-in-cell (PIC) 
simulations show that the GJ density is still an excel- 
lent description of the local charge density in a force- 
free magnetosphere (Spitkovsky, 2006). Note that the 
Goldreich-Julian density does not depend on the specific 
assumption regarding the magnetic field configuration. 
A force-free magnetosphere is boring, with no particle 
acceleration and emission. In reality, however, maintain- 
ing a force-free magnetosphere is not easy. One needs 
to have abundant electron-positron pairs with a number 
density n+ = ngj and a multiplication factor € > 1, 
in order to maintain a net charge density matching the 
GJ density everywhere in the magnetosphere. Without 
copious pair production, deviation from the GJ density 
would be quickly built up even if initially a GJ mag- 
netosphere is realized. This is because the centrifugal 
force drives particles away due to the rapid spin of the 
star. Asa result, various charge deficit regions, or “gaps”, 
where |p| < |pa;| is satisfied, would form in the magneto- 
sphere (Arons and Scharlemann, 1979; Cheng et al., 1986; 
Muslimov and Tsygan, 1992; Ruderman and Sutherland, 


1975). In these gaps, Ej no longer vanishes. Charged 
particles are accelerated and radiate curvature radiation 
or inverse Compton scattering, producing e* pairs via ei- 
ther the yB or yy QED processes (Daugherty and Hard- 
ing, 1996; Hibschman and Arons, 2001; Zhang and Hard- 
ing, 2000). The pairs subsequently redistribute in the Ej, 
forming an opposite Æ field and eventually “screen” the 
original Ej. The magnetosphere then again approaches 
the GJ force-free configuration. Such processes are likely 
unsteady, driving refreshed generation of pairs. Produc- 
tion of pairs has long been regarded as the necessary 
condition to power pulsar radio emission, with the ra- 
dio pulsar “deathline” defined such that pair production 
conditions fail (Ruderman and Sutherland, 1975; Zhang 
et al., 2000). 

Another way of modifying the GJ magnetosphere is to 
introduce a global current J in the magnetosphere (Be- 
loborodov, 2009; T'hompson et al., 2002). In this case, the 
net charge density as observed by a lab-frame observer 
becomes (Thompson et al., 2002) 


Pe = PGI + Ptwist > (125) 


where 


= 


T = - 
Ptwist = —— hQ- [rx (V x B) 2 SO: (rx J) 
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describes a new charge density component to induce a 
twisted magnetic field component around the current 
(Ampere’s law). A twisted magnetosphere can be still 
force-free, but is not in a steady state and would grad- 
ually untwist via dissipation within the twist-supported 
current with a non-zero potential (Beloborodov, 2009). 
Chen and Beloborodov (2017) showed from PIC simu- 
lations that an electric *gap" with unscreened parallel 
electric field can form in a twisted magnetar magneto- 
sphere, which continuously accelerate particles and main 
pair production. Twisted magnetospheres are usually 
discussed within the context of the magnetars after X-ray 
flares, which undergo secular untwisting in an extended 
period of time. 

Recent PIC simulations revealed that besides charge- 
depleted gaps for pair starved magnetospheres, an- 
other promising energy dissipation and particle accel- 
eration site for a pair rich magnetosphere is the equa- 
torial current sheet region outside the light cylinder 
(Kalapotharakos et al., 2018; Philippov and Spitkovsky, 
2018). This region is regarded as a possible new site for 
high-energy emission from pulsars. 

Phenomenological studies and geometric modeling of 
pulsar radio emission suggest that there are potentially 
three types of pulsar radio emission: 


e Inner magnetospheric radio emission: Radio emis- 
sion from old, slowly rotating pulsars is consistent 
with emission from the inner magnetosphere in the 
open field line regions. The double-peak pulse pro- 
file and its "radius-to-frequency mapping" (wider 
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separations at a lower frequencies) as observed in 
a large sample of pulsars strongly support this ge- 
ometric configuration. Modeling suggests that the 
radius of the emission is about 10s of stellar radii 
(Rankin, 1993). 


e Outer magnetospheric radio emission: Young pul- 
sars such as the Crab pulsar have a pair of pulses 
that clearly align with the high-energy (y-ray, and 
X-ray) pulses (Hankins and Eilek, 2007). Since 
the latter has to be emitted from the outer mag- 
netosphere (the predicted high-energy cutoff due 
to yB pair production from inner magnetosphere 
models for y-ray emission was not detected), this 
radio component must be generated from the outer 
magnetosphere or even in the current sheet region 
outside the magnetosphere. 


e Magnetar radio emission: Magnetars are poor ra- 
dio emitters and usually to not emit radio pulses 
during the quiescent state. However, they can be- 
come transient radio pulsars after bursting activi- 
ties. When they emit, the radio pulses sometimes 
show a broader pulse profile and a flat or even ris- 
ing spectrum, in apparent contrast to the pulses 
from normal pulsars (Camilo et al., 2007). SGR 
J1935+2154 was detected by FAST to show a pul- 
sar phase five months after FRB 200418, with 795 
pulses detected in 16.5 hours over 13 days (Zhu 
et al., 2022). Unlike the radio pulses of radio pul- 
sars, these pulses have an opposite phase with re- 
spect to the X-ray pulses from the magnetar. It 
is unclear whether magnetar radio emission shares 
the same origin as one of the two mechanisms op- 
erating in normal pulsars or has its distinct origin. 


FRB emission has a typical luminosity ~ 10 orders 
of magnitude higher than pulsar radio emission. It is 
unclear whether any of the three above mentioned mech- 
anisms can apply to FRBs. 


2. Coherent curvature radiation by bunches 


This mechanism has been widely discussed in both 
the pulsar and FRB fields. Within the pulsar con- 
text, Ruderman and Sutherland (1975) suggested that 
unsteady vaccum gap discharges release “sparks” com- 
posed of secondary electron-positron pairs, which collide 
at a distance of 10s of neutron star radius. Two stream 
instabilities drive the formation of bunches (Melikidze 
et al., 2000; Usov, 1987), which radiate coherently in 
curved magnetic field lines to produce pulsar radio emis- 
sion from the inner magnetosphere. The mechanism was 
found “user-friendly” to account for the phenomenology 
of pulsar radio emission, including the characteristic fre- 
quency, radius-to-frequency mapping, polarization prop- 


erties, etc. (Ruderman and Sutherland, 1975). The for- 
mation and maintenance of the bunches were regarded 
as the main drawbacks for such a mechanism (Melrose, 
1978), but various suggestions to overcome these criti- 
cisms have been discussed in the literature (e.g. Melikidze 
et al., 2000). 

The application of this mechanism to FRBs has been 
discussed by several authors (Cooper and Wijers, 2021; 
Katz, 2014, 2018a, 2020b; Kumar et al., 2017; Lu and 
Kumar, 2018; Lu et aL, 2020; Wang and Lai, 2020; 
Wang et al., 2022b,c; Yang and Zhang, 2018; Yang et al., 
2020b). Because of the extremely high T, of FRB emis- 
sion, some novel aspects of the mechanism have been 
noticed. The key ingredients of such a mechanism can 
be summarized as follows: 


e Characteristic frequency: According to Eq.(94), 
the frequency of curvature radiation is vog ~ 
0.72 GHz 425p;!. For 1 GHz radiation, the re- 
quired electron Lorentz factor is 

Ye c 1101/3 p2/8, (127) 
which is in the range of 10? — 10? for a wide range 
of curvature radius p, from ~ 107 cm (10Ryg) to 
~ 10? cm (around the light cylinder radius). 


e Emission power of a bunch: The emission power of 
an individual electron is 


Pe 2 y*e?c 


E 2 


6 x 1071? erg $71) 42, pz? 


4/3 —2/3 


e (4 
e (7.2 x 107" erg a") 14 "p; (128) 


A bunch of Ne» electrons would emit with a power 
~ N2,P, (Strictly, this is the maximum value, 


Yang and Zhang (2018).). The number Ne» in a 
bunch can be estimated as 


Neb = AgÀne & AgAGnaj & 3 x 107" Ap gus  Gmaj40. 
(129) 

where Ç is the net-charge multiplicity with re- 
spect to the Goldreich-Julian density, and Ay is 
the cross section of the bunch, whose radial size 
is fixed roughly as the wavelength A of the emis- 
sion (Fig.9). The most conservative estimate gives 
Ap min ~ 7(YeA)*, which requires that the trans- 
verse coherence region covers the wavelength in 
the electron comoving frame (Kumar et al., 2017; 
Wang and Lai, 2020). The bunch cross section 
can be in principle much larger, up to the radius 
whose projection in the direction of line-of-sight is 
A, ie. ri4 ~ VroAd (the Fresnel length, Fig.9); 
but is limited to the casually connected region size 
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FIG. 9 An illustration of the shape of bunch. The radial size 
is approximately limited by the wavelength, ie. ~ A. The 
maximum transverse size is at least ~ yeA, but can be as large 
as the Fresnel length ~ rgo. Note that in order to show the 
geometry clearly the bunch size is greatly exaggerated. In 
reality, the distances from the two edges of the bunch to the 
NS as well as ro are similar to each other. 


ri2~ p/y. So one may write 


Ap max c: min [roA, T (o/y)*] 
~ min(9.4 x 108 cm rovg 1,3.1 x 101? em pog’). 
(130) 


Here rg is the distance between the FRB emission 
region and the effective origin of the field line tan- 
gents!?. This is especially the case when field lines 
are nearly parallel in the outer magnetospheres. 


e Observed luminosity: Because of the light propa- 
gation effect as discussed in 8IV.A, the observed 
power of an individual emitting electron is greater 
than its emitted power by a factor of ~ (1 — 
B.cos0)-! ~ 42 (when 0 € 1/7). Considering 
the possible existence of M, independent bunches 
that contribute to the observed luminosity at an 
epoch, one may write the total true luminosity (not 
isotropic equivalent) as 


Le N,N? Pei 
~ (7.8 x 1075 erg s 1) Noe Af orna 10; (131) 


where Mp, A and naj are normalized to their re- 
spective typical values (Ny may be estimated as 


15 The introduction of rg is purely for a geometric purpose to es- 
timate the maximally allowed the size of the bunch. Physically, 
particles are ejected from the inner magnetisphere of the neu- 
tron star and whether there is transverse coherence up to 7r9A 
depends on the detailed particle injection and bunch formation 
processes. 


te 


Ar/X = (3.3x10° cm)Arzvo, where Ar is the depth 
of the field line that contributes to instantaneous 
radiation). Interestingly, vg and p7 are apparently 
cancelled out in Equation (131). This "true" lumi- 
nosity from the model may be compared with the 
beaming corrected luminosity derived from the ob- 
served isotropic luminosity as discussed in §IV.D, 
i.e. 


L ~ Liggmax(0; /4, (452) 1), (132) 


where the solid angle of an individual FRB 9€) is 
written as 710; where 0; is the half opening angle of 
the min-jet. Note that when 6; < y7}, our treat- 
ment is consistent with (Kumar et al., 2017), who 
used a ^2 parameter to make a connection between 
the emitted power and the observed isotropic lumi- 
nosity. Our treatment is more general that includes 
the 6; > 4; ! regime. One can see that for plausible 
parameters, the observed FRB isotropic-equivalent 
luminosity can be reproduced. 


Cooling time and the required Ej: Kumar et al. 
(2017) first pointed out that for models invoking 
curvature radiation by bunches, a steady Ej is 
needed in the magnetosphere to continuously in- 
ject energy to the bunches to maintain the ob- 
served luminosity for the typical FRB duration. We 
reproduce this conclusion in the following with a 
slightly different derivation. Consider a bunch of 
Ne,» electrons radiating coherently. The total en- 
ergy is Ey = Ne»yemec? and the total emission 
power is N2,P,. So the cooling timescale is 


VeMeC* 


= ~ (45x10 7 s)Apgordy naio (133) 


Ne P. 


'This is much shorter than the typical FRB dura- 
tion. In order to maintain FRB emission power, the 
electrons need to continuously gain energy from an 
electric field Ej such that (Ne pe)Eje = N2,F.. 
'This gives 


Ne VP. 1/3 


= —*" * ~ (1.4 x 109 esu)us p7? Apa nanan: 


ec 


(134) 
The existence of such a field is required to apply 
the coherent curvature radiation by bunches to ex- 
plain FRB emission. Kumar and Bošnjak (2020) 
proposed that such an Ej may be provided by the 
propagation of Alfven waves to a charge starved re- 
gion at an altitude of 10s of neutron star radii. Lu 
et al. (2020) argued that such a mechanism can ac- 
count for FRBs with a wide range of luminosities. 
Cooper and Wijers (2021) investigated the maxi- 
mum luminosity of this mechanism by considering 
the effect of the induced current of the emitting 
bunches. They confirmed that the mechanism can 
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generate emission with FRB luminosities. Qu et al. 
(2023) showed that the existence of such an Ey is 
essential to overcome the plasma suppression ef- 
fect for bunched coherent emission (Gil et al., 2004; 
Lyubarsky, 2021). 


Spectrum: The radiation spectrum of coherent cur- 
vature radiation by three-dimensional bunches in 
a realistic pulsar magnetosphere was calculated in 
detail by Yang and Zhang (2018) and Wang et al. 
(2022c). The spectrum is found to be in the form 
of a broken power law separated by a few charac- 
teristic frequencies defined by the length and the 
opening angle of the bunch. The spectral indices of 
different segments depend on the relative ordering 
among the characteristic frequencies, and the high- 
frequency spectral index depends on the power law 
index of the emitting electrons p. The possible self- 
absorption effect by other bunches was studied by 
Ghisellini and Locatelli (2018). If charges are spa- 
tially separated, the shape of the coherent spectrum 
would have a narrower peak than the regular case 
(Yang et al., 2020b). In general, to achieve narrow 
spectra for the bunch models, one needs to invoke 
convolution of the intrinsically broad spectrum of 
individual bunched charges and their spatial distri- 
bution (Katz, 2018a). 


Polarization: Both O-mode and X-mode polarized 
waves can be generated with curvature radiation 
(Kumar et al., 2017; Wang et al., 2012). Nearly 
10096 linear polarization degree is expected if the 
observer within the 1/7. cone of the electron emis- 
sion beam, but circular polarization can develop 
outside the emission cone (Wang et al., 2022b,c) 
(see also Tong and Wang (2022)). Depending on 
the location of the emission region, the polarization 
angle can either display a swing (for an inner mag- 
netospheric location and/or a rapid rotation of the 
magnetosphere), as seen in radio pulsars within the 
framework of the rotating vector model (Radhakr- 
ishnan and Cooke, 1969), or stay nearly flat (for an 
outer magnetospheric location and/or a slow rota- 
tion of the magnetosphere). 


Radius to frequency mapping and frequency down- 
drifting: Wang et al. (2019) showed that there 
is a simple interpretation to the frequency down- 
drifting feature observed in some FRBs. Since 
charged bunches need to be radiation-reaction lim- 
ited within this model (balance between E, acc- 
celeration and curvature cooling), y. may main- 
tain a roughly constant value along field lines. 
Since the curvature radius continuously increases 
as the bunches move away from the magnetosphere, 
the curvature radiation frequency continuously de- 
creases with the increasing height. Suppose several 


bunches in adjacent field lines were launched si- 
multaneously from the base, as the magnetosphere 
rotates, the line of sight always catches emission 
from lower altitudes (and hence, with a higher fre- 
quency) first and the emission from higher altitudes 
(and hence, with a lower frequency) later, so that 
frequency downdrifting should be commonly ex- 
pected!®. Allowing that the bunches can be ejected 
at somewhat different times, Wang et al. (2020b) 
showed that occasionally a frequency updrifting 
FRB may be observed, but the downdrifting pat- 
tern should prevail. This is consistent with the ob- 
servations (Zhou et al., 2022a). 


Despite the success of this simple model to interpret 
a broad range of pulsar and FRB phenomenology, the 
mechanism has been criticized by several authors: 


e Melrose (1978) pointed out that coherent curvature 
radiation by bunches suffers from the difficulties of 
bunch formation and maintenance. The bunch for- 
mation mechanisms have been explored extensively 
in the pulsar field. The common ingredient of the 
models is a two-stream instability (e.g. Melikidze 
et al., 2000), which is likely realized in the vio- 
lent event that powers an FRB. The maintenance 
of bunches is more difficult to realize. Strong re- 
pulsion force within the bunches tend to disperse 
the bunch spatially and radiation reaction tends to 
disperse the bunch in the momentum space (Katz, 
2018a, 2020b). However, since the FRB duration 
is short, the maintenance mechanism only needs to 
apply for a millisecond duration. 


e Lyubarsky (2021) emphasized that the plasma ef- 
fect, which tends to limit brightness temperature 
and is moderately severe for pulsar radio emission 
(Gil et al., 2004), becomes substantial in suppress- 
ing coherent emission from FRBs. If the bunch 
moving with ye is surrounded by a plasma mov- 
ing with y,, Lyubarsky (2021) suggests that the 
emission power of the bunch N?,P. is suppressed 
by a huge factor of the order of 107!°. Qu et al. 
(2023) revisited the arguments of (Gil et al., 2004; 
Lyubarsky, 2021) and found that the plasma sup- 
pression effect is not important in FRB problems. 
If a strong E exists in the emission region, as ex- 
pected in the realistic FRB models (Kumar et al., 
2017), there is essentially no suppression if the 
bunch is in the radiation-reaction-limited regime 
for coherent curvature radiation. 


16 Lyutikov (2020b) later proposed a similar idea to interpret fre- 
quency downdrifting using radius-to-frequency mapping, even 
though the radiation mechanism was not specified. 
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e The curvature radiation spectrum, similar to the 
synchrotron spectrum, might be too broad to in- 
terpret the narrow-band spectrum as observed in 
some FRBs, especially the repeaters. Charge sep- 
aration can alleviate this criticism (Katz, 2018a; 
Yang et al., 2020b). 


3. Coherent ICS emission by bunches, free electron laser and 
linear acceleration emission 


Besides curvature radiation, there are another family 
of models that invoke vacuum-like coherent mechanisms 
that do not depend intrinsically on the dispersive prop- 
erties of the plasma. Within these models, bunched par- 
ticles resonate coherently in some low-frequency waves, 
either of electromagnetic or electrostatic types, and in- 
verse Compton scatter the waves to higher frequencies 
to make FRBs. A relatively simple model is to invoke 
low frequency electromagnetic waves, which might be 
excited near the neutron star surface by crustal oscil- 
lations (Zhang, 2022). Usually it is believed that crustal 
oscillations would excite Alfvén waves, and indeed the 
bulk of the energy that eventually powers an FRB is 
likely carried by Alfvén waves. On the other hand, if a 
small amount of oscillation energy would be converted to 
electromagnetic waves by coherently oscillating charges 
in the near-surface magnetosphere, then the waves (all 
modes for a quasi-parallel configuration and X-mode only 
for a quasi-perpendicular configuration) would penetrate 
through the magnetosphere unimpeded. Suppose that 
there are relativistic bunched charges moving with a 
bulk Lorentz factor y, the low-frequency electromagnetic 
waves with angular frequency wo and frequency vo ~ 104 
Hz would be upscattered to a frequency 


(135) 
(136) 


w = ywo(1 — B cos6,), 
v = (1 GHz) 42 ;vo4(1 — B cos 0;), 


where 0; is the incident angle. Such an inverse Comp- 
ton scattering (ICS) model has been considered to in- 
terpret pulsar radio emission (Qiao and Lin, 1998; Qiao 
et al., 2001; Xu and Qiao, 2001) and its promise to inter- 
pret FRB emission has been discussed in detail in Zhang 
(2022). 

The advantage of the ICS mechanism is that the 
emission power of a single particle, PICS ~ (1.6 x 
10-7 ergs-!)(0Bos)?rg? (where àBo is the oscilla- 
tion amplitude of magnetic field of the electromagnetic 
waves and r is radius of the emission region) is much 
greater than that of curvature radiation, PYR ~ (4.6 x 
10715 erg s) spg ^. As a result, the required degree of 
coherence to interpret the FRB high brightness temper- 
ature is greatly reduced. Indeed, even if one adopts the 
most conservative bunch cross section Ay = 7(4A)?, 
the required Ne» is so low that a charge number den- 
sity of the order of ngaj is already enough to account 


for the FRB luminosity. As a result, the bunches do 
not need to have a large plasma density and the criti- 
cism of bunch emission suppression due to the plasma 
effect (Lyubarsky, 2021) is greatly alleviated (Qu et al., 
2023). Even a small fluctuation in charge number density 
with respect to the background Goldreich-Julian density 
(Yang and Zhang, 2018) would be adequate to produce 
bunched ICS radiation, so that the criticisms of bunch 
formation and maintenance (Melrose, 1978) are also alle- 
viated. The frequency of the scattered waves (Eq.(136)) 
depends on wo, y and 6;. All could be nearly constant 
within an FRB (y is radiation-reaction limited), so the 
bunched coherent ICS mechanism has the advantage of 
generating narrow-band spectra than curvature radia- 
tion, as is observed from most repeating FRBs (Zhou 
et al., 2022a). The frequency down-drifting feature may 
be also produced via a radius-to-frequency mapping fea- 
ture or the intrinsic damping of the low frequency waves 
toward longer wavelengths (Zhang, 2022). This mecha- 
nism can also produce intrinsic circular polarization with 
a proper viewing geometry (Qu and Zhang, 2023). 

In the case that vacuum-like electromagnetic waves are 
not excited from the near surface region, coherent radio 
emission may be still excited by relativistic particles scat- 
tering off an oscillating electric field along the magnetic 
field line via the amplified linear acceleration emission 
(ALAE) or off Alfvén waves via the free-electron laser 
(FEL) mechanism. 

ALAE was introduced in early years by Melrose (1978) 
as a mechanism to replace bunched curvature radiation 
to interpret pulsar radio emission. It was further devel- 
oped by Rowe (1995). It invokes an oscillating Ej along 
the direction of particle motion, with particles radiating 
coherently in such an accelerating field. This mechanism 
has not been investigated in detail within the FRB con- 
text. Since the emission occurs nearly the neutron star 
surface, it is unclear whether the mechanism can pro- 
duce the observed high Tp of FRBs without being ab- 
sorbed or scattered within the inner magnetosphere (e.g. 
Beloborodov, 2021a; Ioka, 2020). 

The FEL mechanism invokes the interaction between 
an Alfvén wave disturbance (also called wiggler) and a 
relativistically moving bunch. This mechanism has been 
studied intensively in the laboratories (e.g. Benford and 
Weatherall, 1992) and was discussed by Fung and Kui- 
jpers (2004) within the pulsar radio emission context. It 
was investigated in detail by Lyutikov (2021) within the 
FRB context. The characteristic angular frequency of 
emission is defined by 


w c Ay? (ck), (137) 


where ky is the wavenumber of the low-frequency wiggler 
waves, similar to the direct ICS case Eq.(135) where ck, 
is replaced by wo. The trajectory of bunched relativistic 
electrons in wiggling Alfvén waves can be solved. The 
resulting emission spectrum has a narrow band, which 
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can interpret the spectral feature of the Crab pulsar and 
narrow spectra of repeating FRBs (Lyutikov, 2021). The 
mechanism is also more powerful than curvature radia- 
tion, and hence, easier to satisfy the brightness temper- 
ature constraint. In general, the FEL mechanism (Lyu- 
tikov, 2021) and coherent ICS mechanism by bunches 
(Zhang, 2022) are intrinsically similar mechanisms and 
share several common features and advantages in inter- 
preting FRB coherent radio emission. 


4. Magnetospheric maser mechanisms 


The bunching mechanisms discussed above do not in- 
voke negative absorption or growth of plasma modes that 
depend on the dispersive properties of the plasma. In this 
section, we discuss several magnetospheric maser mech- 
anisms for coherent radio emission. These mechanisms 
include “vacuum masers" that invoke negative absorp- 
tion of electromagnetic radiation as if in an vacuum, and 
^plasma masers" that invoke growth of plasma modes. 
In order for the latter models to work, several require- 
ments are needed: (1) The plasma should support modes 
whose frequency falls into the observed frequency band 
(e.g. GHz); (2) there should be an unstable particle dis- 
tribution in the relevant frequency band, which should 
resonate the plasma mode; (3) the mode should grow 
rapid enough to reach the desired amplitude to account 
for the high brightness temperature; and (4) the plasma 
mode should eventually escape the region as electromag- 
netic waves. To date, none of such magnetospheric mod- 
els have been found suitable to interpret FRB radio emis- 
sion. As a result, even though the following mechanisms 
have been introduced to interpret coherent radio emis- 
sion from other astronomical objects, so far they have 
not been found successful in interpreting FRB emission 
(see Lu and Kumar (2018) for a critical study of maser 
mechanisms within the context of FRBs). Nonetheless, 
they are listed below for completeness. 


e Relativistic plasma emission: Plasma emission is a 
multi-stage process (Melrose, 2017). The first step 
is to drive a longitudinal Langmuir waves through 
a streaming instability between a fast beam and 
a background plasma. The subsequent stages in- 
clude amplification of the Langmuir turbulence and 
the conversion of the plasma mode to electromag- 
netic waves that could eventually escape the mag- 
netosphere. The characteristic emission frequency 
should be the plasma frequency w, or the boosted 
plasma frequency lw, in the relativistic version. 
Even though plasma emission has been identified 
as the main mechanism producing coherent solar 
radio bursts, it has not been successful in interpret- 
ing pulsar radio emission. The main reason is that 
the instability growth rate is too small due to the 


limitation of the relatively small plasma property 
gradients. 


e Electron cyclotron maser emission: This mecha- 
nism involves plasma maser emission at the non- 
relativistic cyclotron frequency wg and its harmon- 
ics swg, with decreasing amplitudes at higher s 
values. The mechanism was found to be respon- 
sible to the decametric radio emission of Jupiter 
and the auroral kilometric radiation of Earth (Mel- 
rose, 2017). However, within a neutron star mag- 
netosphere (especially for a magnetar), wp is usu- 
ally much higher than the radio frequency so that 
this mechanism is usually not relevant for magne- 
tospheric models!”. 


e Curvature radiation maser: Curvature radiation 
maser is not possible for a rotating dipole because 
there is no solution for negative absorption. How- 
ever, if the pulsar magnetosphere is distorted, un- 
der certain conditions negative absorption would be 
possible for curvature radiation (Luo and Melrose, 
1995). Such a model is not attractive to interpret 
FRBs since it is unclear how the specific magneto- 
spheric configuration might be realized in an FRB 
emitting source. 


e Anomalous cyclotron-Cherenkov and Cherenkov 
drift resonances: Instabilities occur when a dis- 
persion relation has a term whose denominator ap- 
proaches zero, termed as a resonance. In a pulsar 
magnetosphere, maser-type plasma instabilities can 
operate at the anomalous cyclotron-Cherenkov res- 
onance q— kj vj +wp/7 = 0 (7 is the plasma Lorentz 
factor in the pulsar frame) and the Cherenkov drift 
resonance w — kj vy) — kiua = 0 (ua is the drift 
velocity). Even though this mechanism is a plau- 
sible candidate to account for pulsar radio emis- 
sion (Lyutikov et al., 1999), they are not favored in 
interpreting FRBs because the conditions for the 
maser mechanisms to operate either cannot be re- 
alized or demand unreasonable parameters (Lu and 
Kumar, 2018). 


5. Other magnetospheric mechanisms 


Two more magnetospheric coherent mechanisms have 
been proposed to interpret FRB emission, which deserve 
special discussion. 

The first model was proposed by Lyubarsky (2020). 
This model invokes a large scale magnetic perturbation 


17 For slow rotating non-magnetar pulsars, the condition may be 
satisfied. However, the energetics of the neutron star would not 
be large enough to power FRBs. 
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to form a magnetic pulse, which strongly compresses the 
magnetospheric plasma and pushes it away. The pulse 
propagates from the flare site within the magnetosphere 
outwards and eventually reaches the current sheet that 
separates the oppositely oriented magnetic fields beyond 
the light cylinder. The FRB is powered by the enhanced 
magnetic reconnection in the current sheet region. It is 
conjectured that coalescence of magnetic islands in the 
reconnection current sheet produces magnetosonic waves, 
which propagate on top of the magnetic pulse and eventu- 
ally escape as electromagnetic waves. The characteristic 
frequency is defined by the dimension of the magnetic 
islands €a’ where a’ is the width of the current sheet, so 
that 
j c 
w = dw =ô Ea" 
where the primed quantities are measured in the rest 
frame of the magnetic pulse, and ó is the Doppler fac- 
tor of the pulse. In order to match the observed FRB 
frequency, € ~ (10 — 100) is required. The emission is 
polarized along the rotation axis of the magnetar. The 
advantage of this model is that the emission site is slightly 
beyond the light cylinder, which avoids the criticisms re- 
garding FRB propagation within the magnetosphere (Be- 
loborodov, 2021a). PIC numerical simulations of such a 
scenario have been carried out, which show the excite- 
ment of narrow-band GHz emission (Mahlmann et al., 
2022). An alternative radiation mechanism within this 
scenario is bunched ICS of relativistic electrons acceler- 
ated from the current sheet off the low frequency waves 
generated from the inner magnetosphere (Zhang, 2022). 
The second model is the direct electromagnetic wave 
generation from non-uniform pair production across dif- 
ferent field lines. Philippov et al. (2020) showed from 
large-scale two-dimensional kinetic plasma simulations 
that non-steady pair production and screening of electric 
fields along the magnetic field lines by freshly produced 
pairs naturally generate electromagnetic waves that can 
escape the magnetosphere. They propose that such a 
mechanism could be responsible for the coherent radio 
emission of radio pulsars. Yang and Zhang (2021) devel- 
oped an analytical toy model for this process and showed 
that the mechanism can also apply to FRBs given that 
non-steady, non-uniform pair production could be real- 
ized in an FRB environment. They argued that crustal 
cracking of a magnetar could be the engine for such non- 
steady, non-uniform pair production processes. 


(138) 


6. Transparency of FRBs from magnetospheres 


One criticism to the magnetospheric mechanism of 
FRBs is that the FRB waves may undergo strong scat- 
tering by the magnetospheric plasma so that the high 
brightness temperature would not be achievable (Be- 
loborodov, 2021a). Due to their high intensities, FRB 


waves have a large oscillation amplitude for the wave E, 
field (also wave B, field) such that the dimensionless 
amplitude parameter (e.g. Luan and Goldreich, 2014) 


eEw eBy WB.w 
a= ~ = — >1 
MeCW — mec w 


(139) 


in the magnetosphere (see §VIII.C for more discussion on 
the large-amplitude wave effect). This amplitude factor 
a denotes how fast an electron moves in response to the 
waves and when a > 1, the electron speed approaches the 
speed of light. Without an external magnetic field, the 
electron would move under both the oscillating E, field 
and the Lorentz force due to Bw, making a “8” shape tra- 
jectory (Sarachik and Schappert, 1970; Yang and Zhang, 
2020). The electron is accelerated to a Lorentz factor 
of the order of a. The scattering cross section can be 
generally defined as ø = P/S, where P is the emitting 
power and S is the received photon flux. Because of the 
relativistic motion of the electron, the emitted power is 
enhanced by a factor of ~ a? with respect to Thomson 
scattering, so that (Sarachik and Schappert, 1970; Yang 
and Zhang, 2020) 
c c aor. (140) 
In a neutron star magnetosphere with background B, 
the situation is more complicated. When B > Bw, the 
electron motion is dictated by B rather than B,,, so that 
the enhancement of o (Eq.(140)) would not occur. In a 
magnetosphere, since B œ r^? for a dipolar configuration 
and since By, = \/L/cr? x r-! for an EM wave, there 
will be a point where B drops below B,. Recall wg = 
eB/mec, this condition can be also written as 


a > —. 
[i 


(141) 


'The scattering cross section is greatly increased and the 
optical depth is greatly enhanced: 


Tes NOT, ~ UAGIAEB. ous P RSS. (142) 
Here n = £naj, 6 = co/([a(r.)|?ov) is the cross sec- 
tion normalized to a?o4, and re = (B, R3)! (c/ L)/4 ~ 
(4.2x 108 em) L3] ^ BE Rg!”, which is the critical radius 
at which Bw = B. This estimate suggests that the FRB 
waves would indeed become opaque to Thomson scatter- 
ing in a magnetar magnetosphere, if L and £ are large 
(Beloborodov, 2021a). The situation worsens since the 
relativistic motion of electrons in complicated trajecto- 
ries would radiate y-rays, which may produce additional 
pairs to increase the opacity. 

However, the above arguments are based on two as- 
sumptions: (1) the magnetospheric plasma is essentially 
at rest and the angle between wave propagation and the 
local B field, 0p is nearly 90°, i.e. the FRB is trying to 
penetrate through the closed field line region. Qu et al. 
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(2022) argued that both assumptions are likely invalid in 
realistic magnetospheric emission models for FRBs. Var- 
ious mechanisms (the standard pulsar mechanism, Alfvén 
wave propagation, and ponderomtive force acceleration) 
likely drive a relativistically moving plasma in the open 
field line region of a magnetar magnetosphere. The prop- 
agation of the intense FRB waves also tend to align the 
k and B vectors so that Ürp is likely < 1. Both ef- 
fects would reduce the scattering optical depth signifi- 
cantly and it is shown that FRBs are transparent in a 
magnetar magnetosphere even for high-luminosity FRBs 
with a large pair multiplicity, if the plasma Lorentz factor 
Yp > 10? (Qu et al., 2022). 

FRBs are likely associated with X-ray and y-ray pho- 
tons emitted from a magnetar magnetosphere. The 
transparency of FRBs depends on the competition be- 
tween the FRB and X-ray luminosities. Ioka (2020) 
showed that FRB photons can break out of the pair- 
rich magnetopshere with radiation pressure if the FRB 
emission radius is larger than a few tens of NS radii. 
As long as the work done by the FRB waves on the 
e^ is small compared with the initial FRB energy, the 
FRB can successfully break out the magnetosphere. Ioka 
(2020) showed that the breakout condition is satisfied in 
the high Lang low Lx regime. According to this result, 
SGR giant flares may not be associated with successful 
FRBs since the bright X-ray emission would likely choke 
the FRB jet (see Katz (2016) for discussion of alternative 
possibilities of non-detection of an FRB associated with 
SGR 1806-20 giant flare). This is consistent with the ra- 
dio luminosity upper limit of the SGR 1806-20 giant flare 
(Tendulkar et al., 2016). 


C. Relativistic shock models 


The second general type of models invoke relativistic 
shocks to generate coherent radio emission. The term 
“synchrotron maser” has been adopted to describe sev- 
eral very different scenarios. We discuss the three ver- 
sions of the model below, with the decreasing order of 
their relevance to “synchrotron maser”, which inciden- 
tally, is also the reverse order of popularity. 


1. Vacuum synchrotron maser 


The first model is literally “synchrotron maser”. For 
a synchrotron emitting source, the synchrotron absorp- 
tion coefficient can be written in the form of (Ghisellini, 
2017; Lu and Kumar, 2018; Rybicki and Lightman, 1979; 
Waxman, 2017) 


u 1 a dN/dy 


where j,(v, y) is the viewing-angle-dependent emissivity 
for a single electron (in erg s^! Hz^! sr}, different from 


the volume emissivity commonly defined). One may also 
write the net absorption cross section per particle (Ghis- 
ellini, 2017; Ghisellini and Svensson, 1991; Lu and Ku- 
mar, 2018) 


1 1 ð 
as 2mev? «y? Oy 


[^ 3, Cy, p). (144) 


A vacuum maser is possible when either c,,, is negative 
(ju (y, Y) is a decreasing function of y), or a, is positive 
(dN/dy distribution is steeper than y~?, i.e. population 
inversion). 

Ghisellini (2017) found that if the emission region has 
an extremely ordered magnetic field and if the emitting 
electrons have a very narrow distribution for both pitch 
angle and energy, Ca, < 0 is possible in a certain range 
of the viewing angle w > 1/7., where ye is the electron 
Lorentz factor. Even though relativistic shocks are not 
specified in the model discussed by Ghisellini (2017), the 
required magnetic field strength B and electron energy ye 
for the characteristic synchrotron frequency to fall into 
the FRB band are consistent with the typical values for 
shock models. The plasma effect is not important in this 
model, so the emission is the “vaccum” type. 

Even though the mechanism is clean and straightfor- 
ward, the difficulties of the model include how to main- 
tain extremely ordered B field (within 1/7. angle), how 
to accelerate particles to maintain a narrow pitch angle 
distribution (again within 1/ye angle) and how to accel- 
erate particles to maintain a narrow energy distribution. 
Known astrophysical particle acceleration mechanisms, 
e.g. relativistic shocks and magnetic reconnection, usu- 
ally accelerate particles to a power law energy distribu- 
tion and the accelerated relativistic electrons typically 
have a wide angular distribution with respect to the lo- 
cal B field. Perturbations usually introduce wiggles of 
magnetic field lines. As a result, this mechanism may 
not be realized in nature due to the contrived physical 
conditions required. 


2. Plasma synchrotron maser in non-magnetized relativistic 
shocks 


Accelerated particles in relativistic shocks usually have 
an energy distribution dN/dy x ^? with p ~ —2 above 
the minimum Lorentz factor ym. Maser emission is 
therefore impossible for the frequency range defined by 
^ > Ym; Since a, is positive. Nonetheless, population 
inversion (p > 2) may be possible at y < Ym (Sagiv 
and Waxman, 2002; Waxman, 2017). In the extreme 
case, a sharp cutoff of 4-distribution below ym mimics 
a 6-function, which is much steeper than 77. In the fre- 
quency space, maser emission occurs at (Sagiv and Wax- 
man, 2002) 


v < vg. = mins, (vp / vg) ?]vp, (145) 
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where vp» is the modified Razin frequency below which 
the plasma effect becomes dominant, v, = w,/27, and 
vg = wp/27. The relativistic beaming effect for syn- 
chrotron radiation is suppressed because of the role 
played by the refractive index n, (Rybicki and Light- 
man, 1979). The traditional synchrotron radiation is 
suppressed, but the possibility for maser emission is 
opened. This mechanism is a plasma version of syn- 
chrotron maser, and it applies to a weakly magnetized 
plasma with wg « wp. In a hydrodynamical shock, one 
usually defines microscopic parameters e, and eg as the 
fraction of shock internal energy that are distributed in 
electrons and magnetic fields. Observations of GRBs 
show that typically eg < ee < 1 (Kumar and Zhang 
(2015) and references therein). Since v,/vg = wp/wp ~ 
(c«/ep)!?, the condition for plasma synchrotron maser 
emission is satisfied. 

For this model to work, a weakly magnetized central 
engine is preferred. Demanding model parameters to 
satisfy FRB observational constraints, Long and Pe’er 
(2018) showed that neutron stars with surface magnetic 
fields B, < 10!! G is preferred. This is at odds with the 
observational constraint that magnetars are responsible 
to at least some FRBs. Also, since the emission region 
is weakly magnetized, such a model does not predict an 
extremely high polarization degree as is observed in the 
majority of FRBs. As a result, this mechanism, if rele- 
vant, would not be responsible for the majority of FRBs. 


3. Bunched coherent cyclotron/synchrotron radiation in highly 
magnetized relativistic shocks 


Another version of the relativistic shock models invokes 
a highly magnetized upstream. The upstream magnetic 
field lines are highly ordered. As the shock propagates 
into the magnetized medium, magnetic fields are ampli- 
fied and particles coherently gyrate around these field 
lines, forming a “ring” in the momentum space, even 
though they can spread in a wide position space. They 
then radiate coherently as a global bunch at the gyra- 
tion frequency ~ w'4 = eB’/mec, where B’ is the down- 
stream magnetic field strength in the comoving frame of 
the fluid. The observed frequency is Doppler boosted by 
a factor of the bulk Lorentz factor T if the shock moves 
towards the observer relativistically. Such a mechanism, 
even still termed as “synchrotron maser”, is in fact more 
analogous to bunched coherent cyclotron/synchrotron ra- 
diation mechanism (the electron Lorentz factor is typ- 
ically a few), even though bunching occurs in the mo- 
mentum space. The mechanism was introduced to the 
FRB field by Lyubarsky (2014) and studied by various 
teams to interpret FRB observations (Beloborodov, 2017, 
2020; Lu et al., 2020; Margalit et al., 2020a,b; Metzger 
et al., 2019; Yu et al., 2021). The physical process of this 
mechanism has been verified via particle-in-cell numer- 


ical simulations (Babul and Sironi, 2020; Plotnikov and 
Sironi, 2019; Sironi et al., 2021). This mechanism is phys- 
ically robust (with some requirements such as ordered 
B field and cold plasma) and user friendly in interpret- 
ing observations. As a result, it is the most competitive 
mechanism within the relativistic shock model category. 

The features, strengths and weaknesses of this model 
can be summarized as follows. 


e The most important condition for such a mecha- 
nism to operate is the existence of ordered magnetic 
fields in the upstream. Such a feature allows elec- 
trons to gyrate coherently in momentum space so 
that their cyclotron/synchrotron radiation power 
could be coherently added. A commonly suggested 
scenario is that an FRB magnetized pulse collides 
with a magnetized magnetar wind that carries a 
global ordered B field. The FRB is emitted in the 
forward shock region. Tere are two versions of this 
model: the external shock type in which the up- 
stream is an electron-ion wind produced from a pre- 
vious magnetar flare (Metzger et al., 2019), and the 
internal shock type in which the upper stream is a 
relativistic rotationally-powered electron-positron 
pair wind (Beloborodov, 2020). In any case, be- 
cause of the highly ordered magnetic field, very 
high linear polarization degree is expected. The lin- 
ear polarization angle is expected to stay constant 
during each burst, as has been observed in some re- 
peaters (e.g. Jiang et al., 2022; Michilli et al., 2018). 
Because of the same reason, the condition for this 
maser mechanism to operate is also demanding. Ir- 
regularities in the field configuration would greatly 
suppress coherent emission. Also a rapid swing of 
linear polarization angle across individual bursts 
(e.g. Luo et al., 2020a) poses a great challenge to 
such a model. 


e Another condition for such a mechanism to operate 
is that the upstream media should remain “cold”. 
Random motion of electrons in a hot plasma would 
smear up or even destroy the “ring” in the mo- 
mentum space, leading to suppression of coherent 
emission (Babul and Sironi, 2020). As a result, this 
feature poses a constraint on the waiting time of 
successive FRBs within the external shock model. 
Shortly after a collision, both the shocked wind 
and the shocked FRB ejecta would be hot. The 
magnetic field configurations may be also distorted 
due to the irregularities introduced during the col- 
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However, active repeaters also have another peak 
in the waiting time distribution, which is of the 
order of milliseconds (Li et al., 2021b; Xu et al., 
2022). These closely connected bursts, also known 
as “burst storms” (Hewitt et al., 2022) or “burst 
clusters” (Zhou et al., 2022a), pose a challenge to 
the external shock version of this model. This is not 
an issue for magnetospheric models, since different 
pulses are related to different emission regions in 
a rotating magnetosphere as they sweep across the 
line of sight. 


The magnetization parameter 


B? p" 
7 4nloc? Arpe ee) 


is defined as the Poynting-flux-to-kinetic-flux ra- 
tio in the lab frame or the magnetic internal en- 
ergy density (magnetic energy density plus mag- 
netic pressure) over mass density in the comoving 
frame. For an electron-positron plasma, one also 
hasc = wh fut? = w/w where wg = eB/Tm,c 
and wp = (Amne?/Tm,)'?. For this mechanism 
to operate efficiently, the upstream c value should 
be in the Goldilocks zone with a value o ~ 1. At 
smaller c values, since magnetic energy is not dom- 
inant, global magnetic fields are likely subject to 
turbulent perturbation so that the field lines tend 
to be more tangled. The coherent mechanism can- 
not operate efficiently. At higher o values, the frac- 
tion of energy carried by particles reduces (most 
energy is still in magnetic fields) so the efficiency 
of making coherent emission also drops. PIC sim- 
ulations suggest that the maser efficiency scales as 
n ~ 10-807! (Sironi et al., 2021)!8. One inter- 
esting question regarding these models is to ad- 
dress why at the FRB emission radius, o ~ 1 is 
by chance achieved. Pulsar wind theories suggest 
that a pulsar wind with an initial magnetization 
oo > 1 tends to reach o ~ o” at the sonic point 
where the wind speed is as high as the fast sonic 
wave speed so that the magnetic “piston” losses 
pressure to accelerate the outflow (e.g. Li et al., 
1992). Beyond this radius, magnetic acceleration 
is rather slow (unless there exists an external pres- 
sure confinement to maintain a significant magnetic 


lision. If another FRB pulse collides into this rem- 
nant of previous collision, strong coherent emission 
would be likely suppressed. A long waiting time 
of the order of ~ 100 s would be reasonable (Met- 
zger et al., 2019), which is consistent with the sec- 
ond peak of the waiting time distribution of ac- 
tive repeaters (Li et al., 2021b; Xu et al., 2022). 


18 Plotnikov and Sironi (2019) suggested 7 = 7 x 107407? from 
an earlier 1D simulation. The results of 3D simulations by 
Sironi et al. (2021) are generally consistent with the 1D re- 
sults. The difference in -dependence is different frames used. 
The o~?-dependence applies to the shock frame, while the o~1!- 
dependence applies to the downtream frame, which is more rel- 
evant in estimating the maser efficiency (L. Sironi, 2021, private 
communication). 


pressure gradient) so that it is difficult to reduce c 
further down to unity. In general, this mechanism 
predicts a relatively low radio emission efficiency 
n & 1073, suggesting that FRBs should be accom- 
panied by bright high-energy emission in X-rays 
(Margalit et al., 2020a) or optical (Beloborodov, 
2020). It also suggests that the total energy bud- 
get required in the shock models is generally higher 
than that required in the magnetospheric models. 
It turned out that the Galactic FRB 200428 has 
an X-ray-to-radio luminosity ratio of the order of 
~ 10*, which can be accounted for from both mod- 
els (Lu et al., 2020; Margalit et al., 2020a). How- 
ever, active repeaters rFRB 20121102A and rFRB 
20201124A already have a very high total energy 
budget in the radio band during their active burst- 
ing periods. This demands that the radio efficiency 
cannot be much smaller than 10^? in order to sat- 
isfy the total energy budget of magnetars (Li et al., 
2021b; Xu et al., 2022; Zhang et al., 2022). 


In the downstream comoving frame the character- 
istic frequency for maser emission is w% (Sironi 
et al., 2021), which is defined by the strength of 
upstream magnetic field B at the emission radius 
Rypnp, bulk Lorentz factor I, and the central en- 
gine parameters (e.g. surface magnetic field Bs, 
spin period P of the magnetar). Demanding the 
observed frequency Tw to be in the ~ GHz regime 
and combining with other constraints (e.g. RFRB 
as the deceleration radius defined by FRB energy, 
ambient density, and I', and the duration of the 
FRB defined by w = Rprp/cl?), one can place 
some interesting constraints on model parameters. 
This has been done for the Galactic FRB 20200428 
(Lu et al., 2020; Margalit et al., 2020a; Yu et al., 
2021). The general conclusion is that the obser- 
vations can be reproduced, even though some spe- 
cial physical conditions have to be satisfied. To 
overcome such a fine-tuning issue, Metzger et al. 
(2019) argued that the peak of the FRB spectrum 
sweeps across a wide frequency range as it decel- 
erates and the observer only sees them when the 
peak is in the radio band. This idea is also used to 
interpret the spectral down-drifting observed in re- 
peating FRBs (Metzger et al., 2022). On the other 
hand, observationally there is no systematic peak- 
frequency time evolution among adjacent bursts or 
a duration-spectral width correlation to suppose 
this speculation (Zhou et al., 2022a). 


D. Summary 


The discussion in this section can be summarized as 
follows: 
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e There are many coherent radio emission models 


proposed in the literature to interpret FRB emis- 
sion, which generally fall into two categories: mag- 
netospheric (closer-in, pulsar-like) models and rel- 
ativistic shock (far-out, GRB-like) models. Some 
models (e.g. the bunched curvature radiation 
model and magnetized synchrotron maser model) 
have been extensively studied and demonstrated 
ability of interpreting certain FRB data. Some 
other models (e.g. magnetospheric maser mod- 
els and two other versions of shock maser models) 
suffer from some significant criticisms so may not 
be strong candidates to power FRBs. Some other 
models (e.g. bunched ICS and FEL mechanisms, 
reconnection in current sheet, non-steady pair pro- 
duction induced radiation) deserve closer investi- 
gations and confrontation with the data. Current 
observations cannot pin down exactly which mech- 
anism is at play to power FRBs. 


Purely from the theoretical perspective, none of the 
proposed models are free of issues or difficulties. 
Within the magnetospheric models, the bunching 
coherent curvature radiation or ICS models de- 
mand an Ej to continuously inject energy into the 
bunches to satisfy the energy budget constraint. 
The origin of the Ej is not well identified. Vari- 
ous particle - low frequency wave interaction mod- 
els beg the existence of these low frequency waves, 
whose existence can be justified for a neutron star 
model (e.g. through star quakes or glitches) but 
may not be justified all types of central engine 
models (e.g. black hole engines). Magnetospheric 
models in general need to address the opacity of 
high-luminosity bursts, which demands a relativis- 
tically moving plasma in the magnetopshere (Qu 
et al., 2022). The magnetized synchrotron maser 
model in relativistic shocks need to address the ori- 
gin of the demanding requirements including very 
ordered, cold upstream plasma, the Goldilocks c 
value, as well as special model parameters required 
from the data. 


From the observational perspective, data can be 
used to differentiate among some models. In par- 
ticular, the following four criteria (Zhang, 2020c) 
would be helpful: 


1. Polarization angle swings: even though a flat 
PA curve can be accounted for by both shock 
and magnetosphericc models, a significant PA 
swing is consistent with magnetospheric mod- 
els but poses a great difficulty to the shock 
models; 


2. Radio efficiency: a high/low radio emission 
efficiency may offer support to the magneto- 
spheric/shock models. The constraints on ef- 


ficiency may be based on the energy of the 
high-energy counterpart of the FRB (e.g. the 
X-ray burst associated with FRB 200428) or 
theoretically derived total energy budget; 


3. Beaming angle: Magnetospheric models pre- 
dict a narrower emission beam than the shock 
model. Therefore, the identification of nar- 
row beaming for certain FRBs may offer a 
support to the magnetospheric models. Ev- 
idence in support of narrow beaming may in- 
clude the lack of FRBs associated with most 
X-ray bursts from SGR 1935+2154 (Lin et al., 
2020), possible detection of off-beam FRBs, or 
“slow radio bursts” (Chen and Zhang, 2022; 
Zhang, 2021), and frequency-dependent pe- 
riodic window of rFRB 20180916B (Li and 
Zanazzi, 2021). 


4. Rapid variability: Since dt ~ R/T? 
(Eq.(98)), a very small ót would point to- 
ward a small R (if I is constrained) (Beni- 
amini and Kumar, 2020). 'The 60-ns variabil- 
ity (Nimmo ef al., 2021) observed in rFRB 
20200120E from a globular cluster in M81 dis- 
favors a shock origin of the FRB (Lu et al., 
2022). 


Looking ahead, upcoming abundant FRB data may 
shed light on the radiation mechanism of FRBs. It is 
optimistic that data may provide clues on the location of 
the FRB emission (magnetospheres vs. shocks), but the 
identification of the very coherent mechanism(s) may not 
be easy, as the experience in understanding pulsar radio 
emission mechanism speaks itself. The current available 
data seem to support the magnetospheric origin of at 
least some FRBs. It is possible that both magnetospheric 
and shock models operate (e.g. the latter works for the 
most energetic bursts while the former works for less en- 
ergetic ones), but the current data of burst properties 
have not demanded a dichotomy explanation yet. 


VI. SOURCE MODELS 


In this section, we discuss various source models for 
FRBs. Since repeaters seem to be common and since 
there is no proof that intrinsically one-off FRBs (those 
associated with cataclysmic events) exist (but see Mo- 
roianu et al. (2022)), all the subsections except the last 
one in this section discuss sources for repeating FRBs. 
Different from the previous theory review (Platts et al., 
2019) that lists models in a stamp-collecting manner (see 
also the FRB theory Wiki page!?), we attempt to pro- 


19 https: //frbtheorycat.org 
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vide critical comments on these models. §VI.A-VI.C dis- 
cuss the neutron star models, which are the most likely 
models. This is followed by other non-neutron-star astro- 
physical models (§VI.D) and more exotic models (§VI.E) 
for repeating FRBs. Finally, cataclysmic models are dis- 
cussed in 8VI.F. 


A. Magnetars 


The leading source model for FRBs is the magne- 
tar model. Magnetars (Duncan and Thompson (1992); 
Thompson and Duncan (1995, 1996), see also Katz 
(1982)) may be generally defined as neutron stars with 
dipolar surface magnetic fields exceeding ~ 1014 G, but 
there is no clear separation line between magnetars and 
high-B pulsars. Observationally they appear as soft 
y-ray repeaters (SGRs) and anomalous X-ray pulsars 
(AXPs), both having quiescent X-ray luminosities ex- 
ceeding their spin-down luminosities, with the former 
displaying repeated soft-y-ray /hard-X-ray bursts (Kaspi 
and Beloborodov, 2017). Later observations suggest that 
some neutron stars emit SGR-like bursts but have sur- 
face dipolar magnetic fields below 10!^ G (Rea et al., 
2010). These sources may have strong multi-polar mag- 
netic fields near the surface and are also included in the 
magnetar population. There are 30 magnetars currently 
known”, including 16 SGRs and 14 AXPs. In another 
research front in transient astrophysics, a type of mil- 
lisecond magnetar has been hypothesized (Usov, 1992), 
which has been widely discussed as the central engine 
of gamma-ray bursts (GRBs) and superluminous super- 
novae (SLSNe) (e.g. Kasen and Bildsten, 2010; Metzger 
et al., 2011; Woosley, 2010; Zhang and Mészáros, 2001). 

The connection between FRBs and magnetars has been 
discussed by many authors within different contexts. The 
earliest suggestion was by Popov and Postnov (2010) 
who interpreted the Lorimer burst (Lorimer et al., 2007) 
as SGR hyperflares. In the paper by Thornton et al. 
(2013) who reported four additional FRBs, the authors 
discussed several possibilities and pointed out that the 
inferred FRB rate is consistent with the rate of SGR 
flares. The SGR-like model was later further discussed by 
Kulkarni et al. (2014) and Katz (2016). Interaction be- 
tween magnetar flares and ambient wind was introduced 
by Lyubarsky (2014) as a mechanism to generate FRBs. 
Prompted by the discovery of the active repeater rFRB 
20121102A that resides in a dwarf star forming galaxy 
similar to the hosts of long GRBs and SLSNe, Metzger 
et al. (2017) suggested that millisecond magnetars could 
be the engine of active repeating FRBs. This model was 
further developed by Beloborodov (2017, 2020); Margalit 


20 http: //www.physics.mcgill.ca/~pulsar/magnetar/main.html 


and Metzger (2018); Metzger et al. (2019) within the 
framework of the synchrotron maser model. Kumar et al. 
(2017); Yang and Zhang (2018), on the other hand, con- 
sider the requirement of producing FRB emission from 
neutron star magnetospheres, and drew the conclusion 
that the isolated neutron stars that can power FRBs 
are likely magnetars. Wadiasingh and Timokhin (2019) 
proposed that magnetars with a low-twist of magnetic 
fields would initially not have enough pairs to screen Ej 
so that a pair cascade may be triggered to eventually 
power an FRB. Wadiasingh et al. (2020) discussed the 
line of death of FRB emission from magnetars and sug- 
gested that FRB emission is favored in magnetars with 
long periods. Lyubarsky (2020) proposed that enhanced 
magnetic reconnection in the current sheet region of a 
magnetar could power FRBs. Recent developments in 
magnetar FRB models include coherent inverse Comp- 
ton scattering model (Zhang, 2022), the free electron 
laser model (Lyutikov, 2021), and a direct emission model 
from a magnetized shock (Thompson, 2022). Prompted 
by the discovery of FRB 200428 associated with SGR 
J1935+2154 (Bochenek et al., 2020; CHIME/FRB Col- 
laboration et al., 2020), many studies have been carried 
out to investigate how the magnetar model can produce 
FRBs within the magnetophere (Lu et al., 2020; Yang 
and Zhang, 2021; Yang et al., 2020b) or in relativistic 
shocks (Margalit et al, 2020b; Yu et al., 2021). The 
0.286-s period of FRB 20192112A offers a strong support 
to the magnetospheric magnetar models (Beniamini and 
Kumar, 2022; Chime/Frb Collaboration ef al., 2022) at 
least for this special source. 

Various versions of the magnetar models have the fol- 
lowing common ingredients: 


e Energy budget: These models make use of two en- 
ergy reservoirs: either the rotation energy of the 


magnetar 

E, = PI ~ (2.0 x 1076 erg) I44 P7? (147) 
or the magnetic energy of the magnetar?! 
En € TBR ~ (1.7 x 10%" erg) B2, Fig, (148) 


where J is the moment of inertia, P is the spin 
period, B, is the surface dipolar magnetic field at 
the pole, and R is the radius of the neutron star. 
One can immediately see that the rotation energy 
reservoir becomes smaller than the magnetic en- 


ergy reservoir when P > (0.34 BI EIE 


21 This estimate includes the dipolar magnetic field only. Magne- 
tars may store a toroidal magnetic field component, which may 
be much stronger than the poloidal component. So this estimate 
is a conservative lower limit. 


44 


is satisfied. The total energy of bursts for re- 
peating FRBs should be bound by these limits. 
For example, rFRB 20121102A emitted a total 
amount of energy ~ 3.4 x 104! erg in the ra- 
dio band assuming isotropic emission from 1652 
bursts detected in 59.5 hours in a 47-day time 
span (Li et aL, 2021b). Correcting the observa- 
tional duty cycle, the total energy emitted would 
exceed (6.4 x 109 erg) Fy,-17,,+4 assuming a radio 
efficiency 5 ~ 1074 and a global beaming factor 
F, = 0.1. This is a substantial fraction of the mag- 
netic energy available from a magnetar. The mag- 
netar models involving a wide beaming angle and a 
low radiative efficiency (e.g. the synchrotron maser 
model) are greatly constrained by the data. Mag- 
netospheric models invoke a smaller global beaming 
factor and a higher 7,, which are favored (Zhang, 
2020c). 


Energy loss/dissipation rate: The average FRB 
emission luminosity should be bound by the av- 
erage energy loss/dissipation rate of the magnetar. 
The energy loss rate due to magnetic dipole spin- 


down is 
B?R®O4 
È, = g ~ (10% erg 8") Bi asP~*RG, (149) 


and the average energy dissipation rate of the mag- 
netic energy may be estimated as 


Eg = — ~ (3.2 x 10% erg 8 D) Eg arty) 
Td , 


(150) 


where r4 = (10* yr)ta,4 is the characteristic de- 
cay time scale of magnetic fields (e.g. Colpi et al., 
2000). Note that since FRB emission has a very 
low duty cycle (even for very active repeaters), the 
luminosities of individual bursts are not subject to 
these average energy loss/dissipation rate bounds 
as long as the average FRB energy emission rate 
is below this. For example, the active episode of 
rFRB 20121102A occurred in 2019 emitted of the 
order of 10** erg energy in radio band in 47 days. 
Consider that the source has a ~ 160-d period (Ra- 
jwade et al., 2020) and that the source is not active 
in some of the projected cycles, one may roughly es- 
timate the average radio-band energy emission rate 
as ~ 10%! erg/2 yr ~ (1.6 x 10?? erg s-1). This is 
smaller than both È, and Ep. However, the re- 
quirement on 7; is tight, i.e. 


(5 x 1073) Eg 4774,4, Eg budget 
(1.6 x 1074) B, 1, P^R; 9, E, budget. 


3 


Nr > (151) 


Again models with a low 7, are disfavored, unless 
Eg is much larger or 74 much shorter??. 


22 Beloborodov (2017) and Margalit et al. (2019) argued that mag- 


e Triggering mechanism: All magnetar FRB models, 
regardless of how the radio waves are emitted, rely 
on some common trigger mechanisms. One com- 
monly discussed trigger mechanism is crust crack- 
ing at the neutron star surface (e.g. Beloborodov 
and Thompson, 2007; Dehman et al., 2020; Thomp- 
son and Duncan, 2001; Wadiasingh and Timokhin, 
2019; Wang et al., 2018; Yang and Zhang, 2021), 
even though some authors (Levin and Lyutikov, 
2012) suggested that crust cracking may not pro- 
ceed in an abrupt way. An alternative trigger mech- 
anism may be fast ambipolar diffusion in the core 
(Beloborodov, 2017). In any case, oscillations of 
the crust would send Alfvén waves to the magneto- 
sphere, triggering various processes that might be 
related to FRB production (e.g. bunched curvature 
radiation (Cooper and Wijers, 2021; Kumar et al., 
2017; Yang and Zhang, 2018) or inverse Compton 
scattering (Zhang, 2022), direct electromagnetic 
wave generation due to non-uniform pair produc- 
tion (Philippov et al., 2020; Yang and Zhang, 2021), 
enhanced reconnection in the current sheet region 
outside the magnetosphere (Lyubarsky, 2020), as 
well as ejecting magnetic pulses outside the mag- 
netosphere to produce FRBs via magnetized rela- 
tivistic shocks (Yuan et al., 2020). Alternative trig- 
ger mechanisms include sudden magnetic reconnec- 
tion events in the magnetosphere (Popov and Post- 
nov, 2010), sudden discharge of vaccum gaps (Katz, 
2017a), or sudden triggers from an external event 
(Dai, 2020; Zhang, 2017). 


The proposed magnetar models also differ in several 
aspects: 


e Emission site: From small to large distance from 
the neutron star surface, there are four versions of 
magnetar models: 1. models invoking FRB emis- 
sion region inside the magnetosphere, typically 10s 
to 100s of the neutron star radii (e.g. Kumar and 
Bošnjak, 2020; Kumar et al., 2017; Lu et al., 2020; 
Lyutikov, 2021; Wadiasingh and Timokhin, 2019; 
Yang and Zhang, 2018, 2021; Zhang, 2022); 2. mod- 
els invoking the current sheet region outside the 
light cylinder as the FRB emission site (Lyubarsky, 
2020; Mahlmann et al., 2022); 3. models invoking 
internal shocks due to collisions between magnetic 
blobs (Beloborodov, 2017, 2020); and 4. models 
invoking external shocks?? formed when magnetic 
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shells are decelerated by the magnetar wind (Mar- 
galit et al., 2020b; Metzger et al., 2019; Thompson, 
2022). 


e Radiation mechanism: Many mechanisms dis- 
cussed in 8V have been proposed for various ver- 
sions of the magnetar models. The magneto- 
spheric models invoke bunched curvature radiation 
(Cooper and Wijers, 2021; Kumar et al., 2017; Lu 
et al., 2020; Yang and Zhang, 2018), bunched in- 
verse Compton scattering (Zhang, 2022), free elec- 
tron laser (Lyutikov, 2021), or direct EM genera- 
tion due to non-uniform pair production (Philip- 
pov et al, 2020; Yang and Zhang, 2021) as ra- 
diation mechanisms. In the current sheet region, 
magnetosonic waves excited by coalescence of mag- 
netic islands (Lyubarsky, 2020) or coherent inverse 
Compton scattering (Zhang, 2022) are invoked to 
produce FRB emission. In both the magnetic in- 
ternal and external shock regions, the specific ver- 
sion of the synchrotron maser (bunched coherent 
cyclotron/synchrotron radiation) mechanism is in- 
voked to produce FRB emission (Plotnikov and 
Sironi, 2019; Sironi et al., 2021). 


There are many open questions regarding the magnetar 
models for FRBs. Besides the question regarding trigger 
mechanism, emission site, and radiation mechanism dis- 
cussed above, the following are some other examples of 
open questions; 


e Can magnetars produce all FRBs in the universe? 
Shortly after the discovery of FRB 200428 in associ- 
ation with the Galactic magnetar SGR 1935+2154, 
the enthusiasm and confidence of interpreting all 
FRBs in the universe as being generated by magne- 
tars have grown tremendously (e.g. Lu et al., 2020; 
Margalit et al., 2020a). The assumption is that all 
FRBs are intrinsic repeaters. Regular magnetars 
such as those observed in Milky Way may be re- 
sponsible for the apparently non-repeating FRBs 
and those repeaters with a low repetition rate, 
while young magnetars may be responsible for ac- 
tive repeaters observed in cosmological distances. 
The fact that there is no active repeating FRB 
sources from the Milky Way is interpreted as the 
lack of very young magnetars in the Galaxy (or 
if there is any, the FRB emission beam does not 
point towards Earth). The repeaters in association 
with the old population such as globular clusters 
(Kirsten et al., 2022; Nimmo et al., 2022) were in- 
terpreted as a new population of young magnetars 


netars different from the Milky Way known population with a born from accretion induced collapse or mergers of 


le tic field and a short tic d ti l : ^ 
EEE NEN. MENU UENRAQEFIR QN MAC NIA MEE Ld binary neutron stars (Margalit et al., 2019; Wang 
may exit in other galaxies to power active repeaters. 


?3 Under certain conditions, the emission radius of the Metzger et al., 2020a), binary white dwarfs (Kremer et al., 
et al. (2019) model can be smaller than that of the internal shock 2021), or NS-WD binaries (Zhong and Dai, 2020). 


model (Beloborodov, 2020). However, growing evidence suggests that this most 


conservative, “magnetars make them all” sugges- 
tion for the FRB origin may not be adequate to 
account for all the FRB observational data. For 
example, the rFRB 20200120E-like sources may be 
very common. However, none of the known Galac- 
tic magnetars are associated with globular clusters. 
The ‘magnetars make them all” scenario likely runs 
into the event rate issue. The general delay with 
respect to star formation rate required for the in- 
ferred FRB redshift distribution (Hashimoto et al., 
2022; Qiang et al., 2022; Zhang and Zhang, 2022) 
also raises a flag to this simple scenario. 


Does FRB emission favor young or old magne- 
tars? Active repeaters are widely interpreted as 
being produced by new-born magnetars. The argu- 
ments in support of this idea include the association 
of a dwarf star-forming galaxy (Tendulkar et al., 
2017) with rFRB 20121102A, the associations of 
a persistent radio source with rFRB 20121102A 
(Chatterjee et al., 2017; Marcote et al., 2017) and 
rFRB 20190520B (Niu et al., 2021), as well as a 
larger energy reservoir (both magnetic and spin en- 
ergies) and probably a faster decaying rate (Be- 
loborodov, 2017; Metzger et al., 2017) in young 
magnetars. The issues of having very young mag- 
netars as prolific FRB emitters include significant 
free-free absorption and induced Compton scatter- 
ing in a dense environment associated with super- 
nova remnants or pulsar wind nebulae around the 
new-born magnetars. On the other hand, charge- 
starvation seems to be favorable for magnetars to 
make FRBs within their magnetopsheres. Older 
magnetars tend to more easily reach charge starva- 
tion because of the reduced pair production due to 
slow spin and low twist (e.g. Beniamini et al., 2020; 
Wadiasingh et al., 2020; Wadiasingh and Timokhin, 
2019). 


What is the mechanism for Ej, in magnetar mag- 
netospheres? A charge starved region in a magne- 
tar magnetosphere is where Ej is developed and 
particles accelerated. Magnetospheric FRB mod- 
els require the existence of an E to continuously 
supply energy to otherwise rapidly cooling particle 
bunches (Kumar et al., 2017; Zhang, 2022). The 
exact mechanism to generate Ej in the FRB emis- 
sion region is not identified. One possibility is 
that Ej can be developed as Alfvén waves prop- 
agate to the outer magnetosphere where e+ den- 
sity is not sufficient to supply the current required 
to sustain the Alfvén waves (Kumar and Bošnjak, 
2020). Another mechanism is the traditional pul- 
sar mechanism that opens various types of gaps 
in the pulsar magnetosphere (Arons and Scharle- 
mann, 1979; Cheng et al., 1986; Harding and Mus- 
limov, 1998; Muslimov and Harding, 2004; Ruder- 
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man and Sutherland, 1975; Zhang et al., 1997). The 
energetics of these gaps, on the other hand, are 
limited by the spindown power of the magnetars, 
which is not large enough to power FRBs for slow 
rotators. Sudden crust cracking may excite global 
readjustment of the magnetospheric configuration, 
leading to temporarily enhanced gaps with large 
Ej, which could be another mechanism to power 
FRBs. 


What is the role of Alfvén waves? Various FRB 
models invoke Alfvén waves as an important ingre- 
dient. The role of Alfvén waves varies in differ- 
ent models. Kumar and Bošnjak (2020); Lu et al. 
(2020) invoked Alfvén waves as the agent to pro- 
duce Ej at a large enough radius to accelerate 
bunched particles to power FRB emission. Chen 
et al. (2022) questioned this possibility by a nu- 
merical simulation that shows that particles are ad- 
vected without forming a significant charge starved 
region. Kumar et al. (2022a) performed simulations 
in a longer duration and found that an Ey of the 
order of a few percent of the Alfvén wave ampli- 
tude can be indeed generated. Yuan et al. (2020) 
showed that low-amplitude Alfvén waves from a 
magnetar quake propagate to the outer magneto- 
sphere and convert to “plasmoids” (closed magnetic 
loops). The plasmoids are accelerated from the 
star, driving blast waves into the magnetar wind. 
Lyubarsky (2020) invokes Alfvén waves to signif- 
icantly compress the current sheet region outside 
the light cylinder to enhance relativistic magnetic 
reconnection, which may facilitate the generation of 
FRBs in the reconnection region (Lyubarsky, 2020; 
Mahlmann et al., 2022; Zhang, 2022). 


B. Other isolated neutron star models 


Besides magnetars, other types of isolated neutron 
stars have been discussed as the source of FRBs. 


e Giant pulses from young pulsars. Giant radio 


pulses have been observed from some young pul- 
sars, such as the Crab pulsar. The brightest giant 
pulse (GP) observed so far has a peak amplitude 
Symax = 2.2 MJy at 1 GHz and a pulse width 
< 0.4 ns, corresponding to a brightness tempera- 
ture Tj > 10*43 K (thanks to its very short du- 
ration) (Hankins and Eilek, 2007). An immediate 
inference is that similar GPs from nearby galax- 
ies would be detected as FRB-like events by Earth 
observers (Connor et al., 2016; Cordes and Wasser- 
man, 2016). Unlike magnetar-powered FRBs that 
possibly consume magnetic energy of the parent 
star, these GP-like FRBs likely consume spin en- 
ergy of the parent star. Placing GP-emitting pul- 


sars to larger distances suggests that they could 
be detected up to ~ 100 Mpc, but not to larger 
cosmological distances as suggested by the DM ex- 
cess of most FRBs. So FRBs in these models 
are also called *ERBs" — extragalactic radio bursts 
(Cordes and Wasserman, 2016). Now it has been 
confirmed that most FRBs originate from cosmo- 
logical distances greater than 100 Mpc (e.g. Bhan- 
dari and Flynn, 2021; Macquart e£ al., 2020; Ten- 
dulkar et al, 2017). The simplest version of this 
model is incapable of interpreting the data unless 
much brighter GPs are invoked. 


Pulsar lightening. Katz (2017a) argued that the 
FRB phenomonology is similar to atmosphere light- 
ening and conjectured that FRBs are produced 
when vacuum gaps in pulsar magnetospheres break 
down to suddenly drive currents in the magneto- 
sphere. The FRB energetics in this model is also 
limited by spindown power of the underlying pul- 
sar. 


C. Interacting neutron star models 


A number of FRB models invoke neutron stars inter- 
acting with an external agent. These interacting neutron 
star models come in different flavors depending on the 
energy budget that is invoked to explain FRB emission. 
In the extreme versions of the interaction models, the 
ultimate energy comes externally from the gravitational 
energy of a falling object, or the kinetic energy of an ex- 
ternal moving fluid. In milder versions of the models, 
the ultimate energy still comes from the neutron star it- 
self (e.g. the spin or magnetic energy), but the external 
agent may play a role of triggering FRB emission or shap- 
ing the detectability of FRBs. We now discuss several of 
such models in the literature. 


e Comet/asteroid interaction models: One suggested 
way of making FRBs is through interactions be- 
tween comets or asteroids and a neutron star. The 
direct impact model (Bagchi, 2017; Dai, 2020; Dai 
et al., 2016; Dai and Zhong, 2020; Geng and Huang, 
2015; Smallwood et al., 2019) invokes a gravita- 
tional energy budget 


GM 
E, = i ~ (1.9 x 1019 erg) 
M m EE! 
152 
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to power FRBs, where M and R are the mass and 
radius of the neutron star, and m is the mass of 
the small body. One can see that the Galactic FRB 
20200428 from SGR 19354-2154, which has an radio 
luminosity/energy smaller by orders of magnitude 
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than cosmological FRBs, already requires a small 
body mass? m ~ 10% g. Scaling the required mass 
up based on the luminosities to cosmological FRBs, 
one finds that the demanded comet /asteroid mass 
is immense. Take the 2019 active episode of rFRB 
20121102A as an example (Li et al., 2021b). The 
total energy emitted in radio in 1652 bursts de- 
tected in 59.5 hours during 47 days is ~ 3.4 x 10+! 
erg. Counting on the missed FRBs outside of the 
FAST observing window, the total radio energy 
would be ~ 6.4 x 104? erg. Recall that the Galactic 
FRB 20200428 with a radio energy of a few 10°° 
erg (Bochenek et al., 2020; CHIME/FRB Collab- 
oration et al., 2020). The total small body mass 
to power the rFRB 20121102A for that emission 
episode is already a few times of 10?" g, which is of 
the order of the Earth mass. So the comet/asteroid 
collision model is a “very expensive" mechanism 
which consumes a lot of mass. The total mass in 
the Kuiper belt of our own solar system is about 
296 of Earth mass (Pitjeva and Pitjev, 2018). Fur- 
thermore, a significant fraction of comets/asteroids 
are dynamically ejected when a neutron star enters 
the comet/asteroid belt (Smallwood et al., 2019). 
The huge mass budget and the very short wait- 
ing times (as short as several miliseconds) between 
some bursts (Li et al., 2021b; Xu et al., 2022; Zhang 
et al., 2022) essentially rule out the direct impact 
model at least for rFRB 20121102A. 


Another version of the comet/asteroid interaction 
model does not invoke direct impact. Mottez and 
Zarka (2014) suggested that small bodies orbiting a 
pulsar at low orbits could periodically interact with 
the pulsar winds to drive two stationary Alfvénic 
structures called Alfvén wings. The destablisation 
of the plasma by the Alfvén wing's current may 
excite coherent radiation and make FRBs. Mottez 
and Zarka (2014) estimated that a multi-Jy level 
radio burst may be generated if the source is at a 
distance of D — 1 Mpc and if the small body is 
r = 1 AU from the pulsar (flux depends on r^ ?). 
Interpreting cosmological FRBs within this model 
require narrow beaming and low orbits. In general, 
this model is energetically much more efficient than 
the direct impact models, since it does not require 
the small body being destroyed. Nonetheless, since 
the ultimate emission power comes from the spin 
energy of the pulsar, the same energy requirements 
for single neutron stars also applies to this model. 
Observationally, it may be difficult to distinguish 
this model from some isolated neutron star models. 


24 Beaming correction is usually not considered in these models 
because the emission solid angle is expected to be large. 


e Cosmic comb model: Zhang (2017) suggested that a 
sudden interaction between a fluid flow (also called 
an astronomical stream) from a nearby source of 
an otherwise isolated neutron star can make co- 
ehrent radio emission. An FRB is observed by an 
Earth observer when the “combed” magnetosphere 
sweeps across the line of sight. The sources of the 
stream could be energetic events such as super- 
novae, gamma-ray bursts, tidal disruption events, 
or more moderate events such as AGN flares or 
even erratic outflows from a companion. As a re- 
sult, an FRB may or may not be associated with 
bright counterparts depending on the source of as- 
tronomical stream. Note that the specific version 
of this model invoking interaction between a su- 
pernova and a neutron star was proposed earlier 
(Egorov and Postnov, 2009) which was overlooked 
by Zhang (2017). rFRB 20121102A was interpreted 
as a repeatedly combed regular neutron star near 
a massive black hole (Zhang, 2018b). The ultimate 
energy power of this model comes from the kinetic 
energy of the astronomical stream. The kinetic lu- 
minosity received by the neutron star may be esti- 


mated as 
L PN” 
T (5 ) (153) 
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where L is the luminosity of the source of the astro- 
nomical stream, r the distance of the neutron star 
from that source, and P is the spin period of the 
neutron star (cP/27 is the light cylinder radius). 
'The condition for a cosmic comb event to happen 
is that the ram pressure exceeds the magnetic pres- 
sure at the light cylinder, i.e. 


B2 (2nR\° 
2 s 
v > — | — | , 154 

E 8c ( cP ) a 
where p and v are the density and velocity of the 
stream at the interaction radius, B, and R are the 
surface magnetic field and radius of the neutron 
star, respectively. 


The 16-day period of rFRB 20180916B (The 
CHIME/FRB Collaboration et al., 2020) may be 
interpreted as the orbital period of a binary sys- 
tem containing an FRB pulsar (or magnetar) and 
a massive star or neutron star companion (Ioka 
and Zhang, 2020; Lyutikov et aL, 2020). For a 
total mass Mio. = 10M, in the binary, the sep- 
aration between the two stars is 4 x 10/7 cm. 
The kinetic luminosity received by the FRB pul- 
sar is Ly = (3.6 x 108? erg s!) L39P?, where 
L = (109? erg s 1)Lag is the companion's kinetic 
luminosity normalized to its Eddington luminosity. 
Such a luminosity is too small to interpret the re- 
peated FRBs from rFRB 20180916B unless an ex- 
tremely narrow beam or a much greater luminosity 
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than the Eddington value are assumed. As a re- 
sult, the original version of the cosmic comb model 
is not adequate to interpret the observations of at 
least rFRB 20180916B. 


Binary comb models: Ioka and Zhang (2020) pro- 
posed the binary comb model for periodically re- 
peating FRBs. The role of the companion is 
no longer to directly provide the power of FRBs. 
Rather, the interaction between the companion 
wind and the FRB pulsar magnetosphere defines 
a funnel from which FRBs, intrinsically produced 
by the FRB pulsar itself, can escape and be de- 
tected from Earth. The similar scenario was inde- 
pendently proposed by Lyutikov et al. (2020), who 
also displayed the companion-wind-defined funnels 
through numerical simulations. Wada et al. (2021) 
expanded on the funnel mode of Ioka and Zhang 
(2020) and identified two more modes (7-crossing 
mode and inverse funnel mode) to define the FRB 
escaping window for periodic FRBs. Even though 
within the binary comb model the companion wind 
only plays a passive role of defining the detectabil- 
ity of the bursts, it was nonetheless speculated 
(Ioka and Zhang, 2020) that the so-called aurora 
particles entering the magnetopshere of the FRB 
pulsar may play an active role in driving coherent 
radio emission and powering FRBs. 


Magnetospheric interaction models: It is possible 
that direct interaction between the magnetospheres 
of two neutron stars may make FRBs. Possible 
FRB-like electromagnetic field signals have been 
discussed within the context of binary neutron star 
mergers shortly before the merger (e.g. Hansen and 
Lyutikov, 2001; Lai, 2012; Piro, 2012; Wada et al., 
2020; Wang et al., 2016). The commonly discussed 
energy release mechanism is the unipolar effect as 
a neutron star with a weak magnetic field trav- 
els in the magnetosphere of another neutron star 
with a stronger magnetic field. Gourgouliatos and 
Lynden-Bell (2019) studied several configurations 
between the two magnetospheres of the two neu- 
tron stars in the pre-merger phase and discussed 
possible energy dissipation. They mentioned the 
possibility of connecting these interactions with re- 
peating FRBs. Zhang (2020b) showed that for typ- 
ical parameters similar to the double pulsar sys- 
tem (Kramer and Stairs, 2008), strong magneto- 
sphere interactions between the two inspiring neu- 
tron stars occur decades to centuries before the 
merger. He argued that such systems could be ideal 
candidates for producing repeating FRBs through 
magnetic reconnection with the expense of the mag- 
netic energy (and ultimately the spin energy) of 
the two neutron stars. Invoking the beaming effect 
(which is expected from magnetospheric interaction 


induced events), he argued that the energy budget 
in the system is more than enough to power rFRB 
20121102A-like active repeaters. The model has 
several predictions (Zhang, 2020b): 1. The activ- 
ity level elevates with time as the two neutron stars 
get closer with time; 2. Active repeaters could be 
mHz gravitational wave sources detectable by fu- 
ture space-born GW observatories such as LISA 
(Amaro-Seoane et al., 2017), Taiji (Ruan et al., 
2018) and TianQin (Luo et aL, 2016); 3. There 
could be quasi-periodic signals at the orbital pe- 
riod, which is typically 100s of seconds. The en- 
vironment of the globular cluster rFRB 20200120E 
(Kirsten et al., 2022; Nimmo et al., 2022) is con- 
sistent with that of a binary neutron star merger, 
even though models invoking BNSs require that the 
source lasts for a much longer duration, e.g. 10° yr 
(Kremer et al., 2021; Lu et al., 2022). 


White-dwarf-fed neutron star model: Gu et al. 
(2016, 2020) delineated a scenario that invokes a 
compact NS-WD binary in which the WD already 
fills its Roche lobe so that matter from the WD can 
be channelled towards the NS. The authors specu- 
lated that magnetic reconnection may be triggered 
by episodic accretion of WD materials approaching 
the NS surface. Curvature radiation is then envis- 
aged to happen as relativistic particles stream out 
along the magnetic field lines. Such a scenario is 
speculative since known neutron star accreting sys- 
tems (e.g. X-ray pulsars) tend to produce thermal 
emission in the accretion column. In general, it is 
difficult to produce delicate magnetospheric coher- 
ent radio emission in an accreting system. 


D. Non-neutron-star astrophysical models 


Stellar-mass black hole sources. Besides neutron 
stars, the only other kind of objects whose sizes 
are small enough to accommodate millisecond du- 
rations of FRBs are stellar mass black holes. The 
difference between a black hole engine and a neu- 
tron star (e.g. magnetar) engine is that the for- 
mer may or may not have a clean magnetosphere 
as the latter does because of the dirty accretion 
environment, so that the magnetospheric radiation 
mechanisms associated with neutron star models 
may not be straightforwardly applied. Nonetheless, 
Katz (2017b, 2022a) speculated that the accretion 
disk of a black hole may collimate a “funnel” from 
which jet-like emission may be released. In order to 
account for the short duration of FRBs, he further 
speculated that the jet may be rapidly wandering 
and that the duration of the FRB corresponds to 
the duration when the very narrow FRB jet sweeps 


49 


across the line of sight. Katz (2020a) further ar- 
gued that the lack of periodicity in repeating FRBs 
at the typical neutron star spin period favors the 
black hole origin of FRBs (but see Section IV.C for 
counter arguments). These papers did not specify 
the FRB emission site and the coherent radiation 
mechanism. Li et al. (2018a) suggested that the ac- 
cretion system involving a black hole and a white 
dwarf with Roche lobe overflow may launch mag- 
netic blobs and produce FRBs via the synchrotron 
maser mechanism. Sridhar et al. (2021a) proposed 
a detailed model for periodic FRBs invoking an ac- 
creting black hole binary similar to the BH ultralu- 
minous X-ray (ULX) sources. The FRB mechanism 
is hypothesized as the synchrotron maser mecha- 
nism in relativistic, magnetized shocks, similar to 
GRB-like models for magnetars. Such a model 
makes some specific predictions (e.g. some known 
ULX sources will produce FRBs someday). Since 
these models rely on the synchrotron maser mod- 
els as the radiation mechanism, the general theo- 
retical and observational caveats discussed in §V.C 
also apply to these models. 


Supermassive black hole (SMBH) sources: The 
supermassive black holes in the center of galax- 
ies or AGNs have a characteristic timescale, i.e. 
r/c ~ 10°s(M/108 Mo) (where rs = 2GM/c? is 
the Schwarzschild radius of a black hole with mass 
M), much longer than milliseconds. So it is not 
straightforward to invoke a SMBH to power FRBs 
unless emission is confined in a region much smaller 
than event horizon. Nonetheless, some suggestions 
have indirectly made use of SMBHs to power FRBs. 
Romero et al. (2016); Vieyro et al. (2017) proposed 
that FRBs may be produced through interactions 
between a relativistic electron beam from an AGN 
jet and a turbulent plasma. The emitters (called 
cavitons) have a much smaller scale than the SMBH 
so that they can make millisecond-duration bursts. 
In this model, the coherent radio emission is pro- 
duced through Langmuir-wave-driven intense elec- 
trostatic soliton emission, which may be broadly 
defined as one kind of “bunching” mechanisms dis- 
cussed in Section V. Within the framework of the 
cosmic comb model (Zhang, 2017), Zhang (2018b) 
invoked the episodic wind from an SMBH inter- 
acting with a neutron star to interpret the large 
RM and persistent radio emission associated with 
rFRB 20121102A. Gupta and Saini (2018) applied 
the similar scenario to make episodic AGN winds 
to interact with a Kerr stellar-mass black hole to 
launch episodic jets that power FRBs. Wada et al. 
(2021) studied rFRB 20121102A within the frame- 
work of the binary comb model and constrained the 
allowed parameters of the companion of the FRB 


pulsar. They found that a SMBH could be a plausi- 
ble companion of this FRB source. The host galaxy 
data of localized FRBs already rule out AGNs or 
galactic centers as the sources of the majority of 
FRBs. As a result, models attempting to inter- 
pret the bulk FRB population invoking AGNs or 
galactic centers are ruled out. Nonetheless, invok- 
ing AGNs or galactic centers for individual FRB 
sources within the scope of broader models (e.g. 
binary combs) remains possible. 


e Stellar flares: For a short period of time, flares 
from Galactic stars were considered as the sources 
of at least some FRBs (Loeb et al., 2014). The 
suggestion faced the issue of the free-free absorp- 
tion constraint (Luan and Goldreich, 2014) and the 
duration limit (R/c >> 1 ms for stars). The local- 
ization of rFRB 20121102A in a distance galaxy 
(Chatterjee et al., 2017; Marcote et al., 2017; Ten- 
dulkar et al., 2017) quickly put away this model 
and any model invoking origins inside the Galaxy. 


E. Exotic repeater models 


Many FRB models have invoked hypothetical objects 
or phenomena to interpret FRBs, which we summarize in 
this subsection. The confirmation of the existence of any 
of these objects/phenomena would have profound impli- 
cations for astrophysics and physics in general. However, 
since some aforementioned models (e.g. magnetars, iso- 
lated or interacting neutron stars) have provided reason- 
able interpretations to most of the FRB phenomenology, 
we regard these models exotic. Instead of critically com- 
menting on the validity of each model, we simply list 
them below. The only comment on all these models is 
the famous quote by Carl Sagan: “extraordinary claims 
require extraordinary evidence” . 


e Strange quark stars: Strange quark stars are hypo- 
thetical compact stars made up of three flavor (u, d, 
s) quarks (Alcock et al., 1986). Ouyed et al. (2020, 
2021) suggested that conversion from neutron stars 
to quark stars would make quark novae that can ac- 
count for an array of astrophysical transients such 
as GRBs and FRBs. In particular, FRBs are pro- 
duced when the quark nova ejecta chunks collide 
with the ambient medium. Strange stars may have 
a thin normal-matter crust (Alcock et al., 1986). 
Episodic accretion induced collapses of the crust 
have been also suggested to power repeating FRBs 
(Geng et al., 2021; Zhang et al., 2018a). 


e Primordial black holes: Primordial black holes 
(PBHs) are hypothetical black holes formed shortly 
after the Big Bang, which can carry a wide mass 
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distribution not subject to stellar evolution includ- 
ing masses much smaller than Mc. Abramowicz 
et al. (2018) proposed a PBH-NS interaction model 
for repeating FRBs. After a PBH enters the cen- 
ter of a NS, the NS will be accreted and eventually 
swallowed by the PBH. During the process, the NS 
magnetosphere is continuously reconfigured, mak- 
ing repeating FRBs. 


Superconducting cosmic strings: Cosmic strings are 
hypothetical string-like topological structures in 
the universe which are the macroscopic manifesta- 
tion of string solutions in field theories. Like elas- 
tic, current-carrying wires, cosmic strings are en- 
visaged to carry energy, to be dynamically evolv- 
ing and super-conducting. Vachaspati (2008) sug- 
gested that oscillations at the “cusps” (points on an 
idealized string that reach speed of light for a brief 
instant) would radiate FRB-like emission. Other 
suggestions include collisions of string structures 
(cusps and kinks, Cai et al. (2012a,b); Ye et al. 
(2017)), interaction of a current carrying loop in 
the local magnetic field (Yu et al., 2014), and de- 
cay of string cusps (Brandenberger et al., 2017). 


Axion stars, azion clumps, and asion quark 
nuggets: The axion is a hypothetical elementary 
particle and a promising candidate for cold dark 
matter in the universe. If axions exist, it is hypoth- 
esized that they can form gravitationally bound ax- 
ion clumps or axion stars (typically with a mass of 
~ 10-!?M,). It has been suggested that FRBs 
could be generated when axion stars collide with 
neutron stars or black hole accretion disks (Iwazaki, 
2015, 2021). Other ideas include induced collapse 
of axion clumps (“miniclusters”) by the strong 
magnetic field of a compact star (Tkachev, 2015), 
a black hole laser powered by axion superradi- 
ent instabilities or “BLAST” (Rosa and Kephart, 
2018)), and even magnetic reconnection in a neu- 
tron star magnetosphere triggered by the falling of 
“axion quark nuggets (Van Waerbeke and Zhitnit- 
sky, 2019). 


Macroscopic dipole collisions: Thompson (2017a,b) 
conjectured macroscopic, superconducting mag- 
netic dipoles might have formed around the time 
of cosmic electroweak symmetry-breaking. The 
collisions of these “large superconducting dipoles” 
(LSDs) may make tiny explosions to power FRBs. 
The collisions more preferably happen near massive 
black holes where LSDs have higher densities. Both 
repeaters and non-repeaters may be produced with 
this mechanism. 


Dicke’s superradiance: In quantum optics, Dicke’s 
superradiance (DSR, Dicke (1954)) is a phe 


nomenon that occurs when a group of excited (pop- 
ulation inverted) atoms or molecules interact with a 
triggering event (e.g. a light) to radiate coherently. 
The phenomenon was well tested in the laborotory 
(Skribanowitz et al., 1973). Houde et al. (2018) 
hypothesized that DSR can occur in the Galactic 
scale involving ~ 10?? — 10°? entangled molecules 
over distances spanning 100-1000 AU, which can 
power FRBs. Houde et al. (2019) further suggested 
that a pulsar located from ^ 100 pc away from the 
entangled molecules could serve as the trigger for 
DSR. 


e Alien technology: Lingam and Loeb (2017) spec- 
ulated that FRBs may be artificial beam-powered 
light sails of extragalactic aliens. Zhang (20202) 
suggested that FRBs we observe are of astronom- 
ical origins, but communicative extraterrestrial in- 
telligences (CETIs) in the Milky Way galaxy may 
choose to emit FRB-like signals if they want to 
broadcast their existence. The non-detection of 
any artificial FRB-like signals from the galaxy in 
a decade with all-sky radio monitors may place a 
meaningful upper limit on the average emission rate 
of such signals by CETIs in the Galaxy. To produce 
a 1 ms-Jy signal on Earth, the required emission 
power for aliens at a typical distance of 10 kpc is 
~ (1077 W)fe,-3(d/10 kpc)?, where fy ~ 10^? is 
the beaming angle. 


F. Cataclysmic progenitor models 


Even though the majority of detected FRBs are not 
observed to repeat, the cataclysmic progenitor models 
for FRBs are not taken as seriously as repeater mod- 
els. The main arguments against these ideas to become 
the main stream FRB models include: (1) Since the en- 
ergy budget of FRBs is much smaller than the energy 
available in cataclysmic events, and since repeaters have 
been detected, it is essentially impossible to prove that 
the apparently one-off FRBs will never repeat; (2) The 
FRB event rate density is much greater than the rate 
densities of all known cataclysmic events (Ravi, 2019). 
'The most common cataclysmic events in the universe is 


core-collapse supernovae with Roc ~ 10° Gpc^? yr-!, 


whereas the FRB event rate density above 10?" erg s~! is 
a few Rrrp(L > 1037 erg s7!) ~ (107 — 105) Gpe~3 yr7! 
Lu et al. (2020). If any cataclysmic channels are relevant, 
they must only account for a small fraction of FRBs, 
maybe above a particular luminosity (where the rate den- 
sity becomes smaller, e.g. Rrrp(L > 109? erg s^!) ~ 
(3.5 x 104) Gpc^? yr71, (Luo et al., 2020b)), or spread 
out in a wider luminosity range but with negligible con- 
tribution to the observed event rate density. 
Nonetheless, some cataclysmic models are quite attrac- 
tive, since these events are destined to produce brief elec- 
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tromagnetic radiation signals, whether or not they are 
FRBs. Two leading models include the “blitzar” sce- 
nario invoking implosion of supramassive neutron stars 
and various compact binary coalescence (CBC) models 
invoking mergers of neutron stars (NSs) and black holes 
(BHs). We highlight these models below. 


e Blitzars: Supramassive neutron stars (SMNSs) are 
spin-supported massive neutron stars whose non- 
spinning mass already exceeds the maximum NS 
mass allowed by the NS equation of state. The ex- 
istence of an SMNS is therefore temporary. The 
NS will inevitably collapse to a BH as it is spun 
down via magnetic dipolar radiation or even grav- 
itational wave radiation. As the bulk of the NS 
enters the horizon during its collapse to a BH, 
the closed magnetic field lines would detach from 
the star and get ejected (the open field lines pen- 
etrating the hole may stay longer). Falcke and 
Rezzolla (2014) suggested that such a magneto- 
sphere ejection process would power an FRB and 
termed the phenomena “blitzars”. The process was 
numerically simulated (Most et al., 2018), and a 
millisecond-duration episode of significant Poynt- 
ing flux injection was indeed observed, suggesting 
the robustness of the mechanism. Falcke and Rez- 
zolla (2014) envisaged that a significant amount of 
SMNSs may be produced from a few percent of core 
collapse supernovae. Assuming that these SMNSs 
do not carry a strong magnetic field, they suggested 
that collapse happens thousands to million years 
after the birth of the SMSNs, so that no bright 
counterpart is expected in association with FRBs. 
Zhang (2014) pointed out that the so-called “inter- 
nal X-ray plateaus" observed in both long and short 
GRBs are best interpreted as collapse of SMNSs 
born during the GRB events. He therefore sug- 
gested that FRBs should be produced hundreds to 
thousands of seconds after some GRBs if the blitzar 
mechanism is valid. One concern is whether the 
produced FRB can escape the messy environment 
near a GRB. Zhang (2014) suggested that this is 
not a concern since the relativistic GRB jet has 
cleared a funnel to facilitate the propagation of the 
FRB. Since internal plateaus are more commonly 
observed following short GRBs (e.g. Lü et al., 2015; 
Rowlinson et al., 2013, 2010), Zhang (2014) sug- 
gested that there could be intriguing tripple as- 
sociations among FRBs, short GRBs and gravi- 
tational waves. Follow-up radio observations to 
search FRB-like events have been carried out for 
some FRBs (e.g. Bannister et al., 2012; Bouwhuis 
et al., 2020; Palaniswamy et al., 2014; Rowlinson 
and Anderson, 2019), even though no confirmed 
association has been reported (but see Bannister 
et al. (2012) for two untriggered events whose oc- 


curring epochs are consistent with the suggested 
epoch of Zhang (2014)). An intriguing association 
between a non-repeating FRB 20190425A and a 
BNS merger gravitational wave event GW190425 
was claimed by Moroianu et al. (2022) with a 
chance coincidence of ~ 1.9 x 1074. The FRB is 
delayed by 2.5 h from the GW event. Such an as- 
sociation, if indeed physical, is consistent with the 
suggested scenario by (Zhang, 2014). The potential 
host galaxy of FRB 20190425A is also found con- 
sistent with that of a BNS merger (Panther et al., 
2022). 


NS-NS mergers: There are many suggested asso- 
ciations between one-off FRBs and NS-NS merg- 
ers. Most of these suggested processes occur right 
before the merger. Hansen and Lyutikov (2001) 
considered possible electromagnetic precursor emis- 
sion before the merger caused by magnetospheric 
interactions between the two NSs and estimated 
the X-ray and radio luminosities. Piro (2012) and 
Lai (2012) studied the pre-merger magnetospheric 
interaction processes using the unipolar inductor 
model and estimated a brief EM signal with lu- 
minosity up to 10*6 erg s^! within 1 s before the 
merger. Prompted by the discovery of four more 
FRBs (Thornton et al., 2013) that established the 
possible astronomical origin of FRBs, Totani (2013) 
suggested that an FRB may be made right be- 
fore the merger as the magnetospheres of the two 
merged NSs synchronize to orbital motion. The 
unipolar inductor model was specifically applied 
to interpret FRBs by Wang et al. (2016), who 
showed that many of the observed FRB properties 
could be reproduced within this model. Sridhar 
et al. (2021b) proposed a pre-merger FRB model 
invoking specifically the synchrotron maser model. 
The general charged compact binary coalescence 
(cCBC) model (Zhang, 2016, 2019) (see below) also 
applies to NS-NS mergers since NSs are generally 
charged, even though its signal may be outshone 
by other signals discussed above. One common is- 
sue of all the pre-merger NS-NS models for FRBs is 
that the neutron-rich ejecta launched due to tidal 
effect would make the environment “dirty” so that 
the FRB emission may escape only in a small solid 
angle. This further reduces the detection rate of 
these events, making NS-NS mergers incapable to 
interpret the majority of FRBs. 


It is worth restating that the post-merger blitzar 
scenario (Zhang, 2014) gives another possibility 
of NS-NS merger association with one-off FRBs. 
Short GRB observations and theoretical modeling 
suggest that the collapse of the post-merger SMSN 
happen 100 to 104 s after the merger (Gao et al., 
2016; Lasky et al., 2014; Lü et al., 2015; Ravi and 
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Lasky, 2014). This provides a time window of inter- 
est for the search of FRBs associated with NS-NS 
mergers. The GW190425/FRB 20190425A associa- 
tion with a 2.5 hr time difference is consistent with 
this scenario (Moroianu et al., 2022). 


Finally, if a NS-NS merger leaves behind a stable 
massive magnetar. The standard magnetar mech- 
anism may operate and powers repeating FRBs 
(Margalit et al., 2019; Wang et al., 2020a; Yamasaki 
et al., 2018). 


WD-WD mergers: Even though the size of a white 
dwarf is too large to accommodate the millisecond- 
duration of FRBs, Kashiyama et al. (2013) pro- 
posed that mergers of two white dwarfs would lead 
to synchronization of the magnetic fields and pro- 
duce millisecond radio bursts from the polar region 
of a post-merger magnetized white dwarf. This 
model predicts Type Ia supernova - FRB associ- 
ations, which has never been observed. 


NS-BH mergers: For CBCs, if at least one of the 
members is charged, one naturally gets a Poynting 
flux with luminosity rising sharply towards merger 
(Zhang, 2019). Such cCBCs would naturally give 
rise to an FRB-like signal in association with the 
merger (with the FRB observationally delayed due 
to the plasma dispersion). Since NSs (and all 
spinning magnetized objects) are globally charged 
(Michel, 1982), the cCBC signal must exist for neu- 
tron star mergers. Since the NS-NS merger sys- 
tems are messier (see discussion above), NS-BH 
mergers, especially those “plunging events" with- 
out tidal disruption of the NS, are ideal systems to 
observe these cCBCs. Zhang (2019) estimated that 
the total cCBC electromagnetic luminosity of these 
systems can reach 5 x 104? erg s~! for a dimen- 
sionless charge (charge normalized to the critical 
charge defined by the mass of the merging member) 
4 ^ 1077. Another channel to power an EM coun- 
terpart in the plunging NS-BH merger systems is 
to invoke a charged BH due to its interaction with 
the magnetic field of the companion NS, making 
the system a black hole battery (Levin et al., 2018; 
Mingarelli et al., 2015). Dai (2019) showed that 
such a mechanism can produce a detectable EM 
transient (probably in the X-ray band), especially 
if the BH carries a rapid spin. The post-merger 
system of such NS-BH mergers may also release 
the BH spin energy to power a brief EM transient 
(Pan and Yang, 2019; Zhong et al., 2019). 


For non-plunging events, the standard magneto- 
spheric interaction effect may not be important 
(unless the BH is charged, see below). However, if 
jet-like materials can be released before the merger, 
the synchrotron maser mechanism may still operate 


to produce FRB-like events (Sridhar et al., 2021b). 


e BH-BH mergers: BH-BH mergers are not sup- 
posed to produce any EM counterparts unless they 
are either surrounded by matter or electromagnetic 
fields. For the former case, the EM signals should 
typically have long durations unless the matter den- 
sity is close to the nuclear density??. For the latter 
case, brief, FRB-like events may be emitted if at 
least one of the BHs is charged through the cCBC 
process (Zhang, 2016). Such a process has been 
robustly supported from numerical relativity sim- 
ulations (Liebling and Palenzuela, 2016). Other 
related ideas include merger-induced discharge of 
Kerr-Newmann BHs (Liu et al., 2016) and direct 
electric dipole radiation from merging charged pri- 
mordial BHs (Deng et al., 2018). One commonly 
asked question is how BHs attain and retain signifi- 
cant charges. One interesting fact is that collapse of 
a spinning neutron star would leave behind a spin- 
ning BH with charge, i.e. a Kerr-Newmann BH, 
and that the charge does not appear to rapidly de- 
plete (Nathanail et al., 2017). The charge may be 
retained, if a force-free magnetosphere is formed 
around the KN BH. Another possibility is that two 
BHs are merging in a magnetized environment (e.g. 
an AGN disk). The BHs will gain charges via the 
Wald mechanism (Wald, 1974) and launch a Poynt- 
ing flux whose luminosity rapidly increases towards 
the merger (Kelly et al., 2017). 


There are other one-off FRB models. They are either 
exotic or have been significantly constrained by the ob- 
servational data. We list some examples in the following. 
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merger. Such a model would then falls into the 
broad category of NS-NS merger models discussed 
above. 


Primordial black hole evaporation: Primordial BHs 
with mass M, ~ 5 x 10" g (Rice and Zhang, 2017) 
are supposed to evaporate now. Besides making 7- 
rays (Hawking, 1974), these events were suggested 
to emit radio waves as well (Rees, 1977). Keane 
et al. (2012) suggested that this mechanism could 
be one possibility to interpret FRBs. Since the 
total energetics of such an event is ~ 10?! M, ~ 
5 x 10?? erg (Rees, 1977), this model is relevant 
only if FRBs are nearby (e.g. within the Galaxy). 
The cosmological distance of FRBs and their much 
greater isotropic energies rule out this mechanism 
to interpret FRBs. 


White hole explosions: White holes (WHs) are hy- 
pothetical objects in general relativity that have 
opposite properties as black holes. Some quantum 
gravity theories predict black-to-white transition 
as a vast amount of energy falls into a black hole 
reaching the Planck density. The quantum gravity 
pressure would push the matter backwards making 
a white hole. Primordial BHs with mass ~ 1078 g 
are expected to explode today as WHs, which may 
generate non-repeating FRBs (Barrau et al., 2018, 
2014). 


G. Summary 


e Schwinger pairs at the birth of magnetars: Lieu 
(2017) suggested that at the birth of a rapidly spin- 
ning magnetar, abundant pairs would be produced 
from the polar cap region with the Schwinger mech- 
anism, i.e. pairs are drawn from vacuum by strong 
electric fields. The pairs produce FRBs by bunched 
coherent curvature radiation. The star would be 
quickly spun down in milliseconds, and the source 
is not expected to repeat. The main difficulty of 
such a model is that a new magnetar born from 
massive star core collapse is buried inside the ex- 
ploding star and the source is highly opaque when 
the suggested process happens. Even though it is 
not discussed in the original paper, one way to pro- 
duce a naked magnetar might be through a NS-NS 


Even if there have been more than 50 FRB source mod- 
els discussed in the literature, current observational con- 
straints and “Occam’s razor" principle have actually nar- 
rowed down the model options quite significantly. One 
may summarize the state-of-the-art of the source models 
as follows: 


e Repeaters are very likely powered by neutron stars 
that can provide a large enough energy budget and 
frequent enough triggers, either from isolated sys- 
tems or interacting systems. 


e Among isolated neutron star sources, the leading 
candidate is magnetars. However, it is unclear 
whether younger (rapid rotators) or older (slow 
rotators) objects are more favorable to produce 
FRBs. Arguments in favor of both cases have been 
discussed in the literature. More data are needed 
to draw a conclusion. 


25 The free fall timescale, which is the shortest timescale in an ac- 
cretion system is proportional to p-1/? where p is the mass den- 
sity. For a typical stellar density, this timescale is or the order 
of 10? — 10? s (Zhang, 2018c), much longer than milliseconds 
relevant to FRBs. 


Certain interaction processes may play a role in 
defining the observed properties of FRBs and prob- 
ably even in triggering the bursts. 


e Non-neutron-star repeating sources are not needed, 
but not excluded. If these sources exist, likely they 
involve stellar-mass BHs. 


e The existence of a small population of cataclysmic 
FRBs is not robustly established. If they exist, 
blitzars and CBCs are the best guesses. 


VII. ENVIRONMENTAL MODELS 


Since FRBs are extragalactic phenomena, their lo- 
cal environments are not well observed. Nonethe- 
less, the association with a persistent radio source for 
rFRB 20121102A (Chatterjee et al, 2017) and rFRB 
20190520B (Niu et al., 2021) and the large DMyss in 
these two and several other sources suggest that there 
could be compact nebulae near some FRB sources. This 
led to the speculation of the association of a supernova 
remnant or a pulsar/magnetar wind nebula with at least 
some FRB sources. Also the apparent periodicity ob- 
served in rFRB 20180916B (The CHIME/FRB Collabo- 
ration et al., 2020) and probably rFRB 20121102A (Ra- 
jwade et al., 2020) raised the speculation of a binary en- 
vironment at least for some FRB sources. Rapid RM 
variations in some active repeaters rFRB 20201124A (Xu 
et al., 2022) and rFRB 20180520B (Anna-Thomas et al., 
2022; Dai et al., 2022) suggested a dynamically evolving 
magnetized environment of these FRBs. In this section, 
we discuss several environmental models for FRBs. 


A. Persistent radio sources 


Persistent radio sources (PRSs) are associated with at 
least two active repeaters, rFRB 20121102A (Chatter- 
jee et al., 2017) and rFRB 20190520B (Niu et al., 2021). 
These two sources also possess the highest absolute val- 
ues of RM among FRBs: ~ 10° rad m? for the former 
(Michilli et al., 2018) and ~ 10* rad m? for the latter 
(Anna-Thomas et al., 2022; Dai et al., 2022). Leading 
scenarios to interpret the PRSs include supernova rem- 
nants (SNRs), FRB-heated sources, or pulsar (magnetar) 
wind nebulae (PWNe or MWNe), all even mini-AGNs, 
which will be discussed in next subsections. Regardless 
of the detailed models, it is possible to present a generic 
discussion of the emission properties of PRSs, which gives 
a relation between the specific luminosity of the PRS 
LPRS and the RM associated with the FRBs (Yang et al., 
2020a, 2022a). 

The radiation mechanism of PRSs is very likely syn- 
chrotron radiation of relativistic particles from a nebula 
in the vicinity of the FRB source. Since the PRS emis- 
sion is not rapidly varying, it is likely that the PRS does 
not possess a relativistic bulk motion (unlike GRB after- 
glows). One interesting property of synchrotron emission 
is that the specific emission power of each particle only 
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depends on the magnetic field strength B (this is because 
the total emission power P, ~ (4/3)52ec8(B? /81) and 
the characteristic frequency vsr ~ (3/47)y2(eB/mec) 
(Eq.(92)) are both proportional to ^? so the dependence 
on electron Lorentz factor ye is canceled out), e.g. (Ry- 
bicki and Lightman, 1979) 


; (155) 


e 
P, ~ 36 —— B, 
mec 
where ¢ is a factor of the order unity. The peak specific 
luminosity of the PRS can be then estimated as 


3 
e ; B) 
mec 
(156) 
where ne is the total ionized electron number density, fb 
is the fraction that are accelerated to relativistic speeds, 
and the nebula is assumed as a filled sphere with radius 
Rprs and a uniform magnetic field B. From Equation 
(74) and assuming that the observed RM of the source is 
mostly contributed from the nebula (which is reasonable 
since other RM components are typically much smaller 
than the RM of the PRS), one can approximately write 
RM as 


4 
Ly max ind NEP, PM (nias) (v30 
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€ 
[RM| ~ zai (eC) (bj B) Rens, 


(157) 
where CN® is the fraction of ionized electrons that mainly 
contribute to RM and bj = Bj/B < 1 is a fractional 
number to denote the parallel component of the mag- 
netic field. One immediately sees that both Ly max and 
[RM] linearly depend on ne and B, so that their ratio is 
independent of two key parameters of the PRS, i.e. 


Ly max 877 p 2| p2 
RM ^ v3 (zig, ) (me prs: 


Assuming that Rprs does not differ significantly among 
sources, Yang et al. (2020a, 2022a) suggested that the 
reason of a detectable PRS for rFRB 20121102A was be- 
cause of its relatively large |RM|. The non-detection of 
PRSs for the majority of repeating FRBs is simply due 
to their relatively small |RM| values. This suggestion is 
supported by the recent detection of a PRS from rFRB 
20190520B, with a relatively large |RM| (Niu et al., 2021). 


(158) 


B. Supernova remnants 


The association of FRBs with SNRs was suggested 
in the early FRB literature (e.g. Connor et al., 2016; 
Kashiyama and Murase, 2017; Murase et al., 2016; Piro, 
2016). Prompted by the discovery of the PRS of the 
first repeater rFRB 20121102A, Metzger et al. (2017) 
suggested that repeating FRBs are powered by new-born 
magnetars from extreme explosions such as long GRBs or 


superluminous supernovae (SLSNe). Within such a pic- 
ture, an FRB source should be surrounded by an SNR, 
which itself makes radio emission and whose evolution 
dictates the secular DM and RM evolution of the FRB 
source (Margalit and Metzger, 2018; Metzger et al., 2017, 
2019; Piro and Gaensler, 2018; Yang and Zhang, 2017). 
Since such expected coordinated DM/RM evolution is 
not observed, it is now clear that most FRB sources are 
not associated with dwarf star forming galaxies or active 
star-forming regions within the host galaxies, which are 
typical for long GRBs and SLSNe (Bhandari et al., 2020; 
Heintz et al., 2020; Li and Zhang, 2020). The global FRB 
redshift distribution also seems not follow the star forma- 
tion history of the universe (Zhang and Zhang (2022), see 
also Hashimoto et al. (2022); Qiang et al. (2022), but see 
Shin et al. (2022)). So, probably most FRBs are not 
associated with SNRs. In any case, a small fraction of 
FRBs, especially the active repeaters (Chatterjee et al., 
2017; Niu et al., 2021), may be associated with SNRs. 

A dense SNR initially blocks FRBs due to various ab- 
sorption/attenuation processes. The detailed optically- 
thinning conditions depend on the explosion parameters 
(ejecta energy, ejecta mass, ejecta speed) and the ambi- 
ent medium density profile (a constant density Piro and 
Gaensler 2018; Yang and Zhang 2017 or a pre-explosion 
wind profile with n o r^? Metzger et al. 2017; Piro and 
Gaensler 2018). In any case, the general condition is that 
the SNR’s age needs to be of the order of year or decade in 
order to allow FRBs escape freely without suffering from 
various attenuation processes, as discussed in §IV.F. 

The interaction between an SNR blastwave and an am- 
bient medium could be one source of synchrotron emis- 
sion that powers the observed PRS emission as observed 
from rFRB 20121102A and rFRB 20190520B. Metzger 
et al. (2017) applied a parameterized self-absorbed syn- 
chrotron spectrum in the form of F, = Fo(v/v,)9?(1 — 
exp[- (v/v,) -(**9)/?]) to fit the observed spectrum of the 
PRS of rFRB 20121102A and showed that it can roughly 
interpret the data. 

An SNR around an FRB source provides a testable 
prediction about the secular evolution of DM and RM, 
as well as their temporal evolution rates (Metzger et al., 
2017; Piro and Gaensler, 2018; Yang and Zhang, 2017). 
The detailed scaling relations, on the other hand, de- 
pend on several factors, including the density profile of 
the ambient medium, whether the ejecta is fully ionized, 
the density profile of the ejecta itself, as well as the ion- 
ization status of the pre-shocked medium. In general, the 
evolution of an SNR includes four stages: 1. the free ex- 
pansion stage when the ejecta velocity remains constant, 
ie. v x tł; 2. the Sedov-Taylor stage when the ejecta 
accumulates enough mass from the medium and adiabat- 
ically decelerates with the total energy in the blastwave 
conserved; 3. the snow-plow phase when the ejecta de- 
celerates with significant radiative cooling, which is char- 
acterized by momentum conservation; and 4. the disap- 
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pearance stage when the SNR is mixed with ISM. The 
transition radius Rae; between the free-expansion phase 
and the Sedov- Taylor phase occurs when the swept mass 
from the medium becomes comparable to the original 
mass in the ejecta, with the transition time defined by 
taec = Rae /v where v = (2E/M)!/? is the velocity of 
the blastwave with kinetic energy E and mass M. For 
a medium number density n and the mean molecular 
weight um ~ 1.2, one has 


3M I 
c 
AT m 
Arx 
~ (0.43 pc) (Gr) n3”, (159) 
© 
fine = i 
- MNNS 
~ (42 yr) Eg? (x) nj? (160) 


for a constant density medium (Yang and Zhang, 2017) 


and 
M VS 
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tdec = 


E M NP 
~ (1.0 x 104 yr) Ez? (xz) A;! (162) 


for a wind medium (see also the expressions in Met- 
zger et al. (2017) in terms of v rather than E), where 
A = My/(4mv,,) is the wind parameter, A, = A/(5 x 
1011) gcm™! is the typical value of A (Chevalier and 
Li, 1999). Note that in reality the wind profile would 
not extend to infinite distances. It is very likely that 
the medium density profile already returns to the con- 
stant case way before reaching Rae. of the wind model. 
The transition from the Sedov- Taylor phase to the snow- 
plow phase occurs in thousands years after the explosion 
(Draine, 2011; Yang and Zhang, 2017). If FRBs can be 
only made when the neutron star engine is young, only 
the transition from the free expansion phase to the Sedov 
Taylor phase is relevant. 

In general, an SNR may be separated in four regions. 
From outer to inner, they are: 1. unshocked medium 
(ISM or wind); 2. shocked medium; 3. shocked ejecta; 
4. unshocked ejecta or the inner boundary of the ejecta 
if the reverse shock already crosses the shell. Denote 
regions with their respective numbers and the separa- 
tion radii using the two adjacent numbers (ie. Ri» as 
the forward shock radius, R53 as the contact discontinu- 
ity radius, R34 as the reverse shock radius or the inner 
boundary of the ejecta). The total DM from an SNR 


system can be in general calculated as 


R23 Rı2 Ri 
DMsnr = nadr+ | nodr + f nıdr, (163) 
R34 R23 Rio 


where R; is the ionization front in the unshocked 
medium, n; is the total electron number density in re- 
gion i, and f is the ionization fraction in region 1. After 
delineating how R12, R23, R34, ng and ma evolve with 
time, one can derive the t-dependence of DMsnr. 

The strengths of the magnetic field in regions 2 and 3 
can be also estimated by assuming that a fraction ep of 
the internal energy in the respective region is converted 
to (ordered) magnetic fields. Making one additional as- 
sumption that (By) is of the same order as B in the re- 
spective region, one can then calculate the total absolute 
value of RM in the SNR system via 


R23 Rı2 
[RMsun| =} (By,3) nadr + f (Bj 3) nadr 


34 R23 
Ri 
(164) 


and delineate its temporal evolution. 
The predicted scaling laws by various authors and their 
assumptions can be summarized as follows: 


e For a constant density medium, Yang and Zhang 
(2017) assumed that the entire region 3 is ionized 
and obtained 


DMEkR et", dDMSNg/dt x t? (165) 
for the free-expansion phase, and 
DMST p œx 17/5,  qDMST. /dt oc t 9/5 (166) 


for the Sedov- Taylor phase. Note that the DM evo- 
lution scaling does not depend on the medium pro- 
file during the free expansion phase, so that the 
same scaling Eq.(165) also applies to the case of a 
wind medium profile with n « r^? (Metzger et al., 
2017). The assumption of fully ionized region 3 
may be reasonable in view of the existence of a 
repeating FRB source at the center so that any re- 
maining neutral materials between the engine and 
R4 should have been ionized by X-ray emission as- 
sociated with the repeated bursts (see §VII.C for 
more discussion). One interesting finding is that 
DMS$i, increases with time. This is because the 
DM increase rate in shocked medium (Region 2) is 
larger than the DM decrease rate in the unshocked 
medium (Region 1) during the self-similar deceler- 
ation phase. 


e Piro and Gaensler (2018) argued that not the whole 
ejecta is fully ionization. Rather, only the region 
between the reverse shock and the forward shock is 
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ionized?6. Properly following the evolution of the 
reverse shock and assuming an ordered magnetic 
field in the ejecta, they considered the DM and 
RM evolution relations for both a constant den- 
sity medium and a wind medium. For the constant 
density (ISM) case, they obtained 


DM$NR oct, dDMSNR  /dt oc t 9/?,(167) 
|RMEBISM| x 4-1/2.— dg RMEEASM /dt ox 173/1168) 


in the free expansion phase, and 


DMEM œ t?/5, dDMann /dt c 79/5, (169) 
|RMETISM| gc 4-1/5.— dg RMSTISM| /at oc 176/3170) 


in the Sedov-Taylor phase. Note that the scaling in 
the ST phase is the same as Yang and Zhang (2017) 
who assumed full ionization, since in the ST phase, 
the shocked medium (Region 2) is the dominant 
region to contribute to the observed DM. 


e For a wind medium, Piro and Gaensler (2018) ob- 
tained 


DME weit, dDMERR™ /dt x t7?, (171) 
Rie ocd, aR dE œ t7? (172) 


in the free expansion phase, and 


DM 472/3, aM /dt œ t75/3 (173) 
RMSE | o 474/3, dIRMEDTPS | /dt x t77/3174) 


in the Sedov-Taylor phase. 


C. Pulsar wind nebulae & FRB-heated nebulae 


The FRB source, likely a young magnetar, would eject 
a wind through spindown and may eject even stronger 
winds during flaring activities. The wind would interact 
with the surrounding supernova remnant to form a pulsar 
wind nebula (PWN). Such a PWN may play an impor- 
tant role in powering the FRB emission itself through 
synchrotron maser emission (Lyubarsky, 2014; Metzger 
et al., 2019), may contribute to the observed DM or RM 
(Margalit et al., 2019; Metzger et al., 2019), and may 
contribute to the emission of PRS as well. 

Dai et al. (2017) argued that a repeating FRB source 
does not necessarily need to have a surrounding super- 
nova remnant to generate a PRS. The wind from the 


26 This assumption needs scrutiny because a new-born SNR. is 
likely very hot that regions outside the shocked region are also 
likely ionized. FRB-associated X-rays will also ionize any neutral 
atoms in the region. 


FRB pulsar may interact with the surrounding medium 
to form a pulsar wind nebula and power persistent radio 
emission. However, in order to power a detectable PRS as 
observed from FRB 20121102A, the central pulsar needs 
to be rapidly spinning (e.g. P S 10 ms, to allow a large 
energy budget) and does not possess a strong magnetic 
field (to allow a long spindown timescale). 

Yang et al. (2016) noticed that the interaction between 
the FRB ejecta and a surrounding synchrotron nebula 
could play an important role in both nebular emission 
and FRB emission. In particular, for certain parameters, 
the FRB frequency could be below the synchrotron self- 
absorption frequency of the nebula. These FRBs would 
be absorbed and could not reach the observer. Rather, 
they would heat up the synchrotron nebula and make a 
bump in the synchrotron spectrum near the absorption 
frequency. This prediction was found suitable to inter- 
pret the spectrum of PRS of FRB 20121102A after the 
latter was discovered (Li et al., 2020a). 


D. Binary systems 


A widely discussed FRB source environment is binary 
systems, in which the companion (a massive star, an- 
other neutron star, or even a massive black hole) of the 
FRB source (likely a pulsar or magnetar) plays a no- 
ticeable role in shaping the properties of the detected 
bursts. Binary interaction was invoked as one of the 
mechanisms to trigger FRBs within the cosmic comb 
model (Zhang, 2017). The discussion of binary systems 
becomes popular after the discovery of the ~ 16d period 
of rFRB 20180916B (The CHIME/FRB Collaboration 
et al., 2020) as the observed period may be interpreted 
as the orbital period of the binary system. It was quickly 
realized that the massive companion of the FRB pul- 
sar could provide a strong, opaque wind to block FRBs 
in certain directions, so that repeated bursts could be 
only observed in certain orbital phases (Ioka and Zhang, 
2020; Lyutikov et al., 2020). More generally, Wada et al. 
(2021) discussed three possible modes for companion - 
FRB source interactions: (1) The funnel mode is the 
mode in which companion wind is stronger than the FRB 
pulsar wind, so that the latter can only open a funnel as 
the pressures of the two winds balance. The funnel is vis- 
ible by the observer at certain orbital phases (Ioka and 
Zhang, 2020; Lyutikov et al., 2020); (2) The r-crossing 
mode is the mode in which the active window is defined 
by the orbital phases where the optical depth of FRB 
against Thomson scattering, free-free absorption and in- 
duced Compton scattering becomes less than a few (the 
photosphere radius due to induced Compton scattering 
is usually defined by 7 ~ 10 rather than 7 ~ 1). The 
FRB source pulsar crosses the photosphere twice dur- 
ing the orbibal motion and only when the orbit is above 
the photopshere could the FRB emission be observed; 
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(3) The inverse funnel mode is the opposite case of the 
funnel mode, in which the FRB pulsar wind is stronger 
than the companion wind and the active phase is greater 
than half of the period. Zhang and Gao (2020) studied 
various binary systems including one NS companion us- 
ing population synthesis models and found that a 16-d 
period is common and the companion is likely a B-type 
star. The frequency-dependent periofic window of rFRB 
20180916B has been raised as evidence against the sim- 
ple binary comb scenario (Pastor-Marazuela et al., 2021). 
However, several scenarios have been proposed to account 
for the observations within various binary scenarios (Li 
and Zanazzi, 2021; Li et al., 2021c; Wada et al., 2021). 

The complicated RM evolution as well as apparent 
Faraday conversion observed in rFRB 20201124A (Xu 
et al., 2022) does not directly point toward a binary sys- 
tem (due to the lack of periodicity). However, a detailed 
modeling of the polarization properties of the system 
seems to require multiple layers of plasma to contribute 
to RM and radio wave absorption and a binary system 
is a likely possibility to account for the data (Li et al., 
2022a; Wang et al., 2022a; Yang et al., 2022b). 

An extreme version of binary systems is to have the 
FRB pulsar orbiting a massive or even a supermas- 
sive black hole. Zhang (2018b) suggested that rFRB 
20121102 may reside near a supermassive black hole 
whose AGN-like-activities may be powering the persis- 
tent radio emission of the source. It is interesting that 
the parameter space allowed for the binary comb model 
to interpret its ~ 157-d period also prefers a supermas- 
sive black hole as the companion (Wada et al., 2021). 
'The large absolute RM value and sign change observed 
in rFRB 20190520B may be also interpreted by invoking 
a massive black hole in the vicinity of the source (Dai 
et al., 2022). 


VIII. PROPAGATION EFFECTS 


Besides the standard dispersion and Faraday rotation, 
FRB radio waves undergo additional interesting propaga- 
tion effects before being detected on Earth. The propa- 
gation effects may leave imprints on the observed signals 
and observed information may in turn be used to diag- 
nose the physical properties of the medium where FRB 
waves propagate through. 


A. Multi-path effects: scattering, scintillation, and RM 
scatter 


One important feature of radio wave propagation is 
that the observed radio waves at a particular time is 
likely the superposition of rays from multiple paths. This 
is because the frequency-dependent propagation speed of 
radio waves depend on plasma density the waves prop- 


agate through and because the densities along the mul- 
tiple lines of sight likely have fluctuations, mostly likely 
because of turbulence that is ubiquitous in astrophysical 
environments. These fluctuations would spread the rays, 
blur the image, broaden the radio pulse, and smear the 
bandwidth. All these effects are characterized as scatter- 
ing (describing pulse broadening) and scintillation (de- 
scribing intensity fluctuation and bandwidth smearing) 
(Rickett, 1977, 1990). 

Scattering is often manifested as a temporal scatter- 
ing tail in FRB pulses. Let the FRB and a thin plasma 
screen (lens) be located at the angular diameter distances 
D, and Di from Earth, respectively. Let the angular di- 
ameter distance between the source and the screen be 
Dı, which is close to but not equal to D, — D; for cos- 
mological sources. The scattering half angle 0, and the 
scattering timescale 7, can be calculated as (e.g. Cordes 
et al., 2016; Macquart and Koay, 2013; Rickett, 1977, 
1990; Xu and Zhang, 2016; Yang et al., 2022a), i.e. 
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where A is the observed wavelength (longer by a factor 
(1+ z;) than that at the scattering screen), and z; is the 
redshift of the screen (lens). Because in the FRB case 
the screen is usually in the host galaxy, when relevant 
we also write the simpler expression in the last equation 
for the case of D, = Dı. Here there are two important 
length scales. One is the Fresnel scale 
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which is the geometric mean of the effective distance 
Deg = DıDıs/D,s and the rest-frame reduced wavelength 
As = A/[2a(1 + z;j) For a spherical wave, this is the 
transverse scale of the wave front where the light path 
difference is A, at a distance of Deg. 

A more important distance scale is the so-called diffrac- 
tive lengthscale, raim, which is the transverse scale of 
the wave front where the root-mean-square difference be- 
tween the two rays is às. Let us assume that the scatter- 
ing effect is introduced by electron density fluctuations 
that arise from a turbulent cascade and the relevant spec- 
trum takes the power-law form in wave number k (e.g. 
Cordes and Lazio, 2002; Cordes et al, 1985; Macquart 
and Koay, 2013; Rickett, 1977; Xu and Zhang, 2016) 

Pon. (k) E Ok $, 


2n/L < k < 2r/lo, (178) 


where lo and L are the inner (dissipation) and outer (in- 
jection) scales of the turbulent energy, C? is the spectral 
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coefficient (the amplitude of turbulence) that describes 
the significance of the density fluctuations, and f is the 
spectral index, which equals 11/3 for the *Kolmogorov" 
turbulent spectrum but can take a more general value. 
The turbulence is short-wave-dominated when f > 3 and 
long-wave-dominated when £ < 3. From the density vari- 
ance ((dne)”) = f Psn. (k)d?k and L > lo, one can write 
(Xu and Zhang, 2016) 


C2 ~ xg One LP, ES, (179) 
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It is convenient to define a scattering measure as the line 
integration of C? along the line of sight (Cordes and 
Lazio, 2002; Cordes et al., 1985) 


D 
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where in the second equation we have assumed that scat- 
tering only happens in a thin screen with thickness A. 
One can finally write the expression of rgjg in the two 
regimes (Xu and Zhang, 2016) 
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(183) 
where re = e?/m,c? is the classical radius of the electron. 
With all these preparations, one can finally derive the 


observed scattering timescale that has dependence as (Xu 
and Zhang, 2016; Yang et al., 2022a) 
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regardless of the regime of 8. For Kolmogorov turbu- 
lence with 6 = 11/3, the numerical value of the index 
is 20/(8 — 2) = 22/5 = 4.4. The value of 72° depends 
on the SM and for FRB parameters, the contribution of 
7255 from the host galaxy or the immediate environment 
of the FRB source is much greater than those from the 
IGM and from the Milky Way (Cordes et al., 2016; Xu 
and Zhang, 2016). 

With the scattering timescale, one can immediately 
define a scintillation bandwidth 

Av; ~ 1/729 ~ (1 kHz)r-1, (185) 

which is too small to be identified in the observing band 
of the telescopes. On the other hand, scintillation band 
smearing is detected in the radio band, which should 
have a very different origin. For FRBs, the detected 
scintillation bandwidth smearing is likely dominated by 
the multi-path propagation effect within the Milky Way 
galaxy. 


The multi-path effect can also affect the observed po- 
larization properties. For linearly polarized FRB emis- 
sion, the multi-path effect can introduce an RM scatter 
(Feng et al., 2022), i.e. different lines of path undergo dif- 
ferent Faraday rotations so that the final observed emis- 
sion is depolarized (Beniamini et al., 2022; Yang et al., 
2022a). The RM scatter may be estimated as (Yang 
et al., 2022a) 
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where ó(n,Bj)i, is ó0(n&Bj) on the scale of ls, and 
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is the maximum transverse scale of the multi-paths. The 
effect of orm is to introduce a frequency-dependent po- 
larization degree, with the fractional reduction of the lin- 
ear polarization amplitude defined by fmw,aepoi = 1 — 
exp(—2A*o2,,) (Feng et al., 2022; O'Sullivan et al., 2012). 
'This effect presents an interpretion to the frequency- 
dependent linear polarization degree of a sample of re- 
peating FRBs (Feng et al., 2022). 

One prediction of the RM scatter theory is that it is 
positively correlated to the observed scattering timescale, 
ie. opm X TS, with a ~ (0.5 — 0.8) (Yang et al., 2022a). 
'This is qualitatively consistent with the observational 
data (Feng et al., 2022). 

Another mechanism to scatter FRB emission is 
through fragmentation of the FRB waves in a magne- 
tar wind. This may induce additional modulation in the 
emission with a Tse x v~? scattering dependence, which 
is not widely observed (Sobacchi et al., 2022). 


B. Plasma lensing and gravitational lensing 


An extreme version of the plasma multi-path effect is 
plasma lensing (Cordes et al., 2017). In general, a denser 
lens with a positive electron column density would serve 
as a diverging lens, but rays passing through different 
parts of the lens, especially from voids, may converge 
to generate caustics that amplify burst signals. Since 
plasma lenses may be dynamically evolving, the lensed 
bursts can allow different spectral behaviors, in contrast 
to gravitational lensing that retains the spectral shape. 

A simplest model is a 1D Gaussian plasma lens model 
(Clegg et al., 1998) that can be described as DM(x) — 
DM; exp(—2?/22), where zo is the characteristic trans- 
verse scale of the lens and zx is the transverse coordinate. 
Let the transverse coordinates in the source, lens, and 
observer’s planes are £s, £x, and zops, respectively, and 
define dimensionless coordinates us = £s/£0, U = x/xo, 


59 


and uobs = Vobs/o0, the lens equation in geometric optics 
could be expressed as 


u(14 ae") =u (188) 


through the Kirchhoff diffraction integral of the Gaussian 
lens (Cordes et al., 2017). Here 


u = (Dı/Ds)us + (Dis /Ds)Uobs; (189) 
and 
Xr. DM, / Di, Di 
2 190 
acm ( e (190) 


is a dimensionless parameter. The amplification factor 
can be written as 


G — [1 -- o(1— 2u2)e- |“, (191) 
which has a maximum 
Ginax OY zo/Tp (192) 


at the caustics where @ = amin. Cordes et al. (2017) con- 
strained the lens parameters required to have caustics, 
which reads DM;Dj,/z2 = 0.65 pc? AU? cm^?. They 
argued that the apparently more active repetition behav- 
ior of rFRB 20121102A compared with other sources may 
be a consequence of significant plasma lensing. The dis- 
coveries of several more active repeaters cast the doubt 
to interpret all of them with the plasma lensing effect and 
tend to suggest that different FRBs may have different 
active levels and some of them (maybe young magnetars) 
are intrinsically more active than others. Nonetheless, 
plasma lensing may leave certain imprints in FRB ob- 
servations. For example, Er et al. (2020) argued that 
frequency-dependent delay due to the geometric effect 
could be comparable to the dispersion delay, so that the 
measured DM could be overestimated if sigmals propa- 
gate through a high-density gradient clump of plasma. 

Similar to other astronomical objects, FRBs can un- 
dergo gravitational lensing. The high event rate of FRBs 
makes it plausible that lensed FRBs can be detected as 
the detected sample increases quickly with time (Li and 
Li, 2014). Since gravitational lenses are not dynami- 
cally evolving, multi-images of the lensed bursts would 
be more analogous with each other with the a strict de- 
lay timescale for all the bursts from the same repeater 
source. The combination of observing multiple images 
with VLBI and the time delay of the images would allow 
a direct probe of the proper motion of a repeating FRB, 
which will directly constrain the physical conditions at 
the source (Dai and Lu, 2017). 


C. Large-amplitude wave effects 


One unique property of FRB waves, thanks to their 
very high luminosites in radio frequencies, is that at a 


small enough radius from the engine, the amplitude of 
the electromagnetic waves is so large that electrons in- 
teracting with the waves would move with a relativistic 
speed. For an FRB with luminosity L, the Poynting flux 
at a distance r from the source is F = L/(4zr?), which 
can be also written in terms of the EM wave amplitude 
F = cE? /(8m) ~ cB?,/(87). As a result, the wave am- 
plitude can be written as 


2L 
Ey ~ By = \{ > = (820 esu or G)L? rig. (193) 
cr 


One can define a dimensionless parameter 


a= = (194) 
MeCw w 
of a wave for its amplitude (where wz, = eBwu/Mec = 
eEw/mec), which is essentially the dimensionless oscilla- 
tion velocity vosc/c of an electron in response to the wave 
when a « 1l. Plugging in the typical FRB parameters, 
one has 


a = 231? rust. (195) 


One can see that for an L = 104? erg s^! FRB, the ampli- 
tude factor is a >> 1 when r « 10!? cm. In such a large- 
amplitude wave regime, a series of propagation effects 
not shared by low-amplitude radio waves are introduced. 
Similar effects apply to laboratory lasers which can have 
very large intensities to reach the relativistic regime. The 
importance of the large-amplitude effects within the con- 
text of FRBs was first pointed out by (Luan and Goldre- 
ich, 2014) and later discussed by various authors within 
various contents (e.g. Beloborodov, 2020; Gruzinov, 2019; 
Kumar and Lu, 2020; Lu and Phinney, 2020; Yang and 
Zhang, 2020). In analogy to the large-amplitude wave 
effects for laboratory lasers, Yang and Zhang (2020) sys- 
tematically studied the large-amplitude effects for FRBs, 
which can be summarized as follows: 


e Enhancement of emission cross section. In the 
a > 1 regime, an electron moves in a “figure-of- 
eight” trajectory because besides the traditional 
harmonic motion due to the oscillating Ew, it is 
also affected by the Lorentz force from the oscil- 
lating B, (Sarachik and Schappert, 1970). In the 
oscillation-center rest frame, the electron moves 
with a Lorentz factor 5' = a/V2. Similar to syn- 
chrotron radiation, the emission power of the elec- 
tron is P ~ a?Prp, where Pr = e*E2/3m26 is 
the received power given by the Thomson formula. 
Considering the Poynting energy flux in the waves 
is S = cE? /81 and that the cross section is defined 
as c = P/S, one gets (e.g. Yang and Zhang, 2020) 


P 
a aor. (196) 
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With the existence of a background magnetic field 
B, as is the case of FRBs emitted from a magnetar 
magnetosphere, the problem becomes more compli- 
cated. In the inner magnetosphere where B > B,, 
is satisfied, the large amplitude effect is suppressed, 
since the electron is confined by the much stronger 
background B. In a dipolar field, one has B œ r~3, 
which decays faster than B, œ r-l. As a re- 
sult, the large-amplitude effect would become im- 
portant when B becomes smaller than B, (Be- 
loborodov, 2021a,b). Detailed numerical results 
suggest that c/c is typically greater than a”, with 
a dependence on the angle between the wave vec- 
tor k and the B vector and the relationship between 
wp/w and a (Qu et al., 2022). When the plasma 
is streaming outwards relativistically, bright FRBs 
can propagate through it and escape the plasma 
successfully (Qu et al., 2022). 


Transparency of strong waves. In the a > 1 and 
weak magnetic field (far away from magnetosphere) 
regime, the dispersion relation for a circularly po- 
larized wave is modified as (Yang and Zhang, 2020) 


w2 


w = ke +, (197) 
y 


where y = (1 + a?/2)!/?. This effectively reduces 
the near-source plasma frequency by a factor of 
y7, or reduce the plasma density by a factor of y. 
This would reduce the DM contribution from the 
vicinity of the FRB source (e.g. within 1 AU for a 
L = 10? erg s^! burst) by a factor of ~ y, making 
the FRB more transparent (Lu and Phinney, 2020; 
Yang and Zhang, 2020). 


Relativistic self-focusing. In the a >> 1 regime, 
the non-linear refractive index is n, = c/vp = 

1— w? /y(a)w?, which is intensity-dependent. 
Consider a beamed FRB with a decreasing inten- 
sity from the center. The propagation effect nat- 
urally “squeezes” the light, making the FRB more 
beamed (Yang and Zhang, 2020). Such an effect is 
especially important for a high-density emitter, e.g. 
in the synchrotron maser scenario. The squeezing 
effect becomes negligible in a magnetosphere envi- 
ronment (Lyutikov, 2020a). 


Ponderomotive force electron acceleration in wake- 
field waves. An electromagnetic pulse with a non- 
uniform energy density (which is the case of an 
FRB) would exert a ponderomotive force (F, = 
—mec?V(1 + (d2))!/? in the relativistic regime, 
where d = eA/m,c? (A is the vector potential, i.e. 
B = V x A) is a dimensionless vector whose am- 
plitude is comparable to a) to the ambient plasma. 


Electrons would be more easily expelled away from 
equilibrium due to the radiation pressure, form- 
ing an oscillating electrostatic field in the plasma. 
This is the so-called wakefield wave. Such a field 
would accelerate electrons. However, such an effect 
is too small to be observational interesting (Yang 
and Zhang, 2020). 


IX. FRBS AS ASTROPHYSICAL AND COSMOLOGICAL 
PROBES 


Regardless of their physical origin(s), FRBs are 
tremendous cosmic probes that can be used to study 
various problems in astrophysics, cosmology, and even 
fundamental physics. In this section we summarize some 
proposed applications of FRBs as cosmological probes. 
Reviews on these subjects can be also found in Bhandari 
and Flynn (2021) and Xiao et al. (2021). 


A. Missing baryons: Q, and fiaw 


Most of the following probes make use of the salient 
feature of the (DMigm) — z relation (Eq.(10)), which 
makes a connection between two observables, DM and 
z. The complication is that there are multiple compo- 
nents that contribute to DM (Eq.(9)). However, in most 
cases, DMiqm is the dominant term. If one can prop- 
erly deduct other components, one can directly measure 
OQ» fiam from the data (Eq.(10)). This has been done with 
a small sample of FRBs (Macquart et al., 2020). The re- 
sults are consistent with indirectly inferred Qa from cos- 
mic microwave background and Big Bang nucleosynthesis 
measurements (Boesgaard and Steigman, 1985; Planck 
Collaboration et al., 2020). This solves the long-standing 
^missing baryon problem" and suggests that the majority 
of the missing baryons are in the intergalactic medium. If 
one adopts the best-fit Q, from the CMB measurements, 
one can directly constrain friem. The results inferred 
from FRBs (Li et al. (2020b) and Figure 4 right panel of 
this review) are generally in agreement with the previous 
results using other methods (Fukugita et al., 1998). 


B. IGM inhomogeneity 


Equation (10) is an average relationship. For indi- 
vidual lines of sight, the measured DM at the same z 
could be very different because the IGM is inhomoge- 
neous. Numerical simulations (McQuinn, 2014) showed 
that the standard deviation o[DM] of the DM distribu- 
tion ranges from 180 to 400 cm^? pc at z = 1 pending 
on whether the “missing” baryons lie around the virial 
radius of 1014 — 10? M; halos or further out. Jaroszyn- 
ski (2019) showed ~ 13% scatter of DM at z = 1 and 
~ 7% scatter at z = 3 using Illustris simulation, see also 
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Takahashi et al. (2021). Macquart et al. (2020) presented 
a sample of 8 FRBs with z measurements, which indeed 
showed a large scatter and the authors expected that the 
range of scatter should increase with redshift. The cur- 
rent data with 21 z-known FRBs (Fig.4 right panel of this 
review) do not show such a trend. Li et al. (2019) recon- 
structed the DM-z relation for nearby FRBs using the 
observed optical galaxy data and the halo baryon distri- 
bution models and found that the inferred DMjay values 
for individual FRBs indeed deviate significantly from the 
predicted values based on the average relation Eq. (10). 
A more detailed study of FRB 20190608 making use of 
both optical and radio data led to a reconstruction of the 
cosmic web along the line of sight (Simha et al., 2020). 
With a much larger sample of localized FRBs with z mea- 
surements, the scatter of the DMraw — z relation will be 
mapped directly from the data. This scatter is also very 
important to decide how good FRBs are to serve as other 
types of probes as discussed below. 


C. Circum-galactic medium 


Individual galaxies are surrounded by a circum- 
galactic medium (CGM), which is the gas surrounding 
the galaxies outside their disks or ISM but inside the 
virial radii. The properties of the CGM are poorly stud- 
ied. The amount of mass in the CGM would affect the 
scatter of the DMraw — z relation. FRBs can probe the 
CGM directly, either for the halo of our own Milky Way 
Galaxy or the halo of foreground galaxies along the line of 
sight of some FRBs. Low DM FRBs from nearby galax- 
ies can be used to directly constrain DMhajo of the Milky 
Way (Prochaska and Zheng, 2019). Analyses of the radio 
data of FRB 20181112 posed strong constraints on the 
properties of the halo of a foreground galaxy, which has 
low net magnetization and turbulence (Prochaska et al., 
2019). The studies in this direction will flourish as more 
data are accumulated. 


D. FRB host galaxy and the surrounding medium 


Another uncertainty that hinders the application of the 
DM iam — z relation to probe the universe is the DM con- 
tribution from the FRB host galaxy as well as the imme- 
diate medium around the FRB source. Both are poorly 
known and difficult to measure because they are degen- 
erate with DMygm, which itself has a large uncertainty. 
Nonetheless, DMyost and DMs;- have been studied from 
different aspects. Theoretically, Xu and Han (2015) sim- 
ulated the DM distributions for three types of FRB hosts 
and different viewing angles. The DM contribution from 
a dense medium (e.g. supernova remnant) around FRBs 
has been extensively modeled (Metzger et al., 2017; Piro 
and Gaensler, 2018; Yang and Zhang, 2017). Observa- 


tionally, some FRBs with an apparent excess DM (e.g. 
rFRB 20121102A, Tendulkar et al. (2017); and rFRB 
20190520B, Niu et al. (2021)) have shown evidence of 
a large DMpost + DMsrc. Information of the host galaxy 
type and relative position of the FRB in its host galaxy 
(e.g. Bannister et al., 2019; Bhandari et al., 2020; Ten- 
dulkar et al., 2017; Xu et al., 2022) can also help to esti- 
mate the DM contribution to the host galaxy. If one as- 
sumes that the DMpost J- DM. of a large sample of FRBs 
follow a normal distribution (which may be the case if 
the outliers such as rFRB 20190520B are removed), the 
average DM contribution from the host/source may be 
inferred statistically using the observed DM-fluence re- 
lation (Yang et al., 2017) or DM-z relation (Li et al., 
2020b). With 5 FRBs with z measurements, Li et al. 
(2020b) estimated the local value of DMpost + DMsgre as 
~ 107722 pe cm-? (the measured value is smaller by 
a factor of (1+ z)). The larger sample of the current 
21 FRBs leads to the similar constraint (see Fig.4 right 
panel and related discussion). With a large enough sam- 
ple, DMhost +DMsrc can be also directly inferred through 
differential increase of the observed DMg with z (Yang 
and Zhang, 2016). From cosmological simulations, it was 
found that DMhost is redshift-dependent and the me- 
dian value ranges from ~ 35 pc cm ? at z = 0.1 to 
~ 106 pc cm? at z = 1.5 (Zhang et al., 2020b). 


E. Dark energy 


Suppose a large sample of FRBs are localized and 
z measured, the IGM inhomogeneity and host/source 
DM contribution can be better quantified. This would 
open an opportunity to compare the data with different 
(DMicgm) — z models and constrain relevant model pa- 
rameters. The first exciting prospect is to use FRBs to 
constrain the evolution of the universe, in particular, the 
nature of dark energy as delineated by the E(z) function 
in Equation (7). Simulations (Gao et al., 2014; Walters 
et al., 2018; Zhou et al., 2014) suggest that depending on 
the degree of IGM inhomogeneity, meaningful constraints 
on dark energy may be achieved with a large enough 
sample, especially in combination with other cosmologi- 
cal probes such as Type Ia supernovae, cosmic microwave 
background, and baryon acoustic oscillations. The chal- 
lenges for robustly extracting distance and the quantita- 
tive estimates of the systematics control needed for FRBs 
to be competitive distance probes have been discussed by 
Kumar and Linder (2019). 


F. Reionization history 
Another prospect of using the (DMjay) — z relation as 


cosmological probes is to probe the reionization history 
of the universe. This is because the observed DM is only 
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contributed by free electrons. The relation (Eqs.(10) and 
(11)) carries the ionization fraction for both H and He 
(Deng and Zhang, 2014; Zheng et al., 2014). Theoreti- 
cal modeling and observational constraints suggest that 
He might be fully ionized at z ~ 3 (Zheng et al., 2014), 
whereas H is ionized at z > 6 (Fan et al., 2006). The de- 
tailed ionization history, especially that of H ionization 
in the so-called *dark ages", is not well constrained, and 
FRBs can potentially probe it directly. Detailed simu- 
lations (Bhattacharya et al., 2021; Caleb et al., 2019b) 
showed that He ionization from z — 3 to z — 6 can be 
differentiated with 1.6 x 10? — 10* FRBs. For H reion- 
ization, the epoch of reionization may be constrained via 
an observed DMmax or 40 FRBs detected at redshifts 
z € (6,10) (Beniamini et al., 2021). 


G. Large-scale structure and turbulence 


With the DM and spatial distribution of a large sam- 
ple, one can perform a study of the angular correlation of 
DMs for FRBs, extracting their structure function and 
correlation function to probe the large-scale structure 
(Shirasaki et al., 2022) or even turbulence at very large 
scales. The pre-CHIME sample showed a preliminary ev- 
idence of possible large-scale turbulence (Xu and Zhang, 
2020), which is not confirmed with the larger CHIME 
sample (Xu et al., 2021). Nonetheless, the results are 
broadly consistent with the statistical modeling of the 
cosmological DM from numerical simulations (Takahashi 
et al., 2021). Rafiei-Ravandi et al. (2021) found a statisti- 
cally significant cross-correlation between CHIME FRBs 
and galaxies in the redshift range z € (0.3, 0.5). 


H. Host/source and intergalactic magnetic fields 


Besides using DM to perform various constraints, a 
combination of DM and RM may place a constraint on 
magnetic fields under ideal situations. Similar to Eq.(9), 
one may decompose the observed RM to several terms 


RMhost + RMsre 

(1+2)? 
(198) 
where RM;,, is the contribution from the Earth iono- 
sphere that gives a measurable small contribution, and 
the (1+ z)? correction factor in the last two terms comes 
from the 0 = \?RM relation, where 0 is the polarization 
angle. In general, the observed RM is likely dominated 
by the near-source medium, which is likely a dynamically 
evolving magnetized environment (Feng et al., 2022; Luo 
et al., 2020a; Michilli et al., 2018; Xu et al., 2022). The 
observed DM, on the other hand, is dominated by the 
IGM term. As a result, RM/DM is not a good probe of 
the average DB, along the line of sight (unlike pulsars). 
One may remove the Milky Way and IGM contributions 


RM = RMion + RMmw + RMiem + 


to DM (with the caveat of a large uncertainty in DMjcm), 
and estimate the average line-of-sight magnetic field in 
the host and source (most likely in the source region) as 


RMhost,obs T RMsrc,obs 
DMhost,obs T DMsrc,obs 


| ~ (1.284G)(1+ 2) , (199) 
with the observed RM as a proxy of the numerator. The 
derived (Bj) values for FRBs are on average consis- 
tent with those of pulsars and magnetars observed in 
the Milky Way (Wang et al., 2020c), with the excep- 
tion of rFRB 20121102A, which has a much higher value 
(Hilmarsson et al., 2021). 

Another way of estimating B near the FRB source is 
to make use of the observed variations of DM and RM, 
i.e. 


ARM 


zn (200) 


Bw (1.234G)(14- 2) 


This already assumed that the variation of Bj is not 
the dominant factor for RM variations. For rFRB 
20201124A, the detection of significant ARM and the 
non-detection of ADM led to a constraint of By > 
0.2mG (Xu et al., 2022). Some repeating FRBs (e.g. 
rFRB 20190520B) show significant RM reversals, sug- 
gesting the reversal of the magnetic field directions. In 
such cases, one has (Yang et al., 2022b) 


óRM |. 5DM | óBj 


RM DM B, eu 


Since both óRM/RM and óBj/B, are of the order of 
unity, the value of Bj cannot be constrained. 

If the RM contribution from the host and source is 
small, or its behavior can be well quantified for a large 
FRB sample, one may combine the observed DM and RM 
information to make a constraint on the poorly known 
IGM magnetic field. Hackstein et al. (2019) showed that 
less than 100 FRBs from magnetars in a stellar-wind en- 
vironment hosted by starburst dwarf galaxies at z 2 0.5 
would be able to differentiate different IGM magnetic 
field models. Recent observations of more complicated 
FRB surrounding medium in terms of RM variations 
(Luo et al., 2020a; Michilli et al., 2018; Xu et al., 2022) 
and RM scatter (Feng et al., 2022) make it difficult to 
correct for the dominant RM contribution from the near- 
source region, rendering constraining the IGM magnetic 
fields much more challenging. 


|. Additional probes with gravitationally lensed FRBs: Ho, 
Qk, and dark matter 


The high event rate of FRBs makes it likely to detect 
gravitationally lensed FRBs in the future. These lensed 
sources, especially the lensed repeating sources, offer new 
opportunities to probe cosmology using FRBs. Thanks 


63 


to their very short durations, the time delays between 
the images can be measured with an unprecedented pre- 
cision. Since the gravitational lensing geometry involves 
the measurements of the angular diameter distances of 
the source and lens, which depend on the Hubble con- 
stant Ho through z and the curvature of the universe Qk, 
lensed FRBs can be used to directly measure Ho and Qk 
(Li et al., 2018b). Simulations showed that with about 
10 lensed repeating FRB systems, Ho can be measured 
to a sub-percent precision level and Qg can be measured 
to a precision of ~ 0.076 in a model-independent manner 
(Li et al., 2018b). 

FRBs can be micro-lensed by massive compact halo 
objects (MACHOs) which have been proposed as one 
type of contributor to the dark matter. For MyacHo Z 
20Mo, the delay time would be longer than one millisec- 
ond. If such lensed events are observed, one FRB with a 
single pulse would be observed as a double-pulse (lensed 
by one MACHO object) or triple-pulse (lensed by a MA- 
CHO binary) bursts. The non-detection of these events 
would place an upper limit on the abundance of these 
MACHOs (Muñoz et al., 2016; Wang and Wang, 2018). 
As a type of MACHO, the abundance of primordial black 
holes is already loosely constrained using the CHIME 
catalog database (Zhou et al., 2022b). The constraints 
will be further improved as the FRB sample continues to 
grow. 

A search for lensed FRBs has been carried out with the 
1st CHIME FRB catalog with no detection (Kader et al., 
2022). This posed a novel constraint on the abundance of 
primordial black holes (Leung et al., 2022). Connor and 
Ravi (2022) forecast the detection rates of gravitational 
lensing of FRBs with delay timescales from microseconds 
to years, corresponding to a wide range of the lens mass 
spanning fifteen orders of magnitude. 


J. Fundamental physics: Weak equivalence principle, 
photon mass, and Lorentz invariance violation 


Thanks to their very short durations, FRBs have been 
also suggested as probes for fundamental physics because 
of the lack of spreading in time in contrast to the predic- 
tions of some theories. 

The first test is Einstein's weak equivalence principle 
(WEP), which states that all point-like structureless par- 
ticles fall along the same path within a gravitational field. 
This is the foundation of the general theory of relativ- 
ity, a geometric description of gravitation. According to 
this principle, photons with different energies from the 
same source should travel with the same trajectory with 
the same speed to reach the observer. In the parame- 
terized post-Newtonian (PPN) description, the deviation 
from the WEP is the PPN parameter y deviating from 1. 
FRBs cannot be used to directly constrain y, but can be 
used to test the difference of y values between two fre- 


quencies vı and v2, which are usually the boundaries of 
the detection frequency band (Wei et al., 2015). Thanks 
to their large distances and short durations, one can con- 
strain Ay with FRBs to be as small as 10715 — 10720 
(Hashimoto et al., 2021; Tingay and Kaplan, 2016; Wei 
et al., 2015; Xing et al., 2019). 

Another interesting constraint FRBs can offer is the 
photon mass (Bonetti et al., 2016; Wu et al., 2016). If 
photons indeed have a non-zero rest mass, the lower- 
frequency photons (with a lower *Lorentz factor") should 
travel slightly slower than higher-frequency photons. The 
duration of an FRB therefore presents an absolute max- 
imum delay due to such an effect. With a more sophis- 
ticated method by combining non-zero photon mass de- 
lay and the plasma dispersion delay (it turns out that 
the two dispersion relations have the similar forms with 
different normalization factors and slightly different z- 
dependences), a more stringent constraint can be reached 
with FRBs of known redshifts, especially with a sample 
of z-known FRBs using a Bayesian approach (Shao and 
Zhang, 2017). The most stringent upper limit of the pho- 
ton mass posed by FRBs already reached m., X 5x 10-45 
g (Bonetti et al., 2016; Shao and Zhang, 2017; Wu et al., 
2016; Xing et al., 2019). 

Another widely discussed fundamental physics con- 
straint is Lorentz invariance violation due to the delay of 
high-energy photons as they travel through the foam-like 
space in very small scales. The effect is most significant 
at high-energies, so that short-duration GRBs are much 
more suitable to pose meaningful constraints than FRBs, 
which have very low photon energies. 


X. PROBLEMS AND PROSPECTS 


The rapid progress in the FRB field is accompanied by 
many open questions, which continue to drive the field 
forward. We discuss three most pressing questions at the 
time when this review is written. 


A. Do all FRBs repeat? 


This question is interesting from both observational 
and theoretical aspects. Observationally, it is much more 
difficult to prove that an FRB does NOT repeat than it 
does. If you have not detected a repeated burst from the 
source yet, it could well be that (1) it repeated but the 
telescope has missed it; (2) it repeated, but the burst is 
below the telescope sensitivity; or (3) it simply has not re- 
peated and the waiting time is longer than the observing 
time. If one adopts the two sub-bursts of FRB 20200428 
from the Galactic magnetar SGR 19354-2154 (Bochenek 
et al., 2020; CHIME/FRB Collaboration et al., 2020) as 
one burst, then the source may not be regarded as a 
repeating FRB source yet (many repeated radio bursts 
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from the source are not bright enough to be detected 
as FRBs at cosmological distances) even though we are 
certain that it should be an FRB repeater, because the 
magnetar source itself did not show significant difference 
before and after the FRB and there is no reason why the 
physical conditions to make FRB 20200428 would not 
be satisfied again to make another one. Theoretically, 
this question is very interesting because it is related to 
whether any of the cataclysmic FRB models are relevant. 

There have been great efforts in addressing this ques- 
tion. 1. From the observational side, although repeater 
bursts are found to display some interesting character- 
istics (e.g. longer duration, down-drifting subpulses, 
narrower spectra (CHIME/FRB Collaboration et al., 
2019b)), there is still no definite clue to suggest that ap- 
parent non-repeaters are indeed different. 2. Machine- 
learning methods have been proposed to differentiate re- 
peaters and non-repeaters (Chen et al., 2022; Luo et al., 
2023; Zhu-Ge et al., 2022), and the results seem to sug- 
gest that most apparently non-repeating FRBs are in- 
deed different from the repeating bursts. 3. A statistical 
study of the observational properties of repeaters and ap- 
parent non-repeaters suggested that there might be two 
populations (Zhong et al., 2022). 4. With limited data 
in the pre-CHIME era, arguments have been made that 
rFRB 20121102A is much more active than any other 
non-repeaters (Caleb et al, 2019a; Palaniswamy et al., 
2018) and so there might be two distinct classes (or at 
least two classes of repeaters with distinct activity lev- 
els). These arguments need to be revisited with the uni- 
form, much larger database from CHIME. 5. One inter- 
esting test is to study the observed fraction of repeating 
sources from all FRBs, Fyops. Ai et al. (2021) showed 
that if there indeed exist non-repeaters and if repeaters 
repeat forever, F} obs should approach a maximum af- 
ter a certain observing time (when most repeaters are 
discovered) and then decline with time afterwards. How- 
ever, uncertainties in the repetition rate and its distribu- 
tion in repeaters make this criterion not clean. In some 
cases, the required time to reach the maximum is longer 
than astronomers' timescale (e.g. longer than 1000 yr). 
When the lifetime of the repeaters is considered, there is 
essentially no achievable maximum within astronomers' 
timescale. In any case, continuously monitoring Fy obs 
may provide important clues to address this open ques- 
tion. The long-term CHIME observations seem to sug- 
gest a constant F; ops over time (Z. Pleunis, 2022, Cornell 
FRB workshop), which may suggest the existence of non- 
repeating FRBs. 

In long terms, besides refining the above analyses with 
a much larger data set, a detection of an FRB robustly 
associated with a cataclysmic event (e.g. a gravita- 
tional wave event) would offer a strong support to the 
existence of these special types of FRBs. The plausi- 
ble GW190425/FRB 20190425A association (Moroianu 
et al., 2022) might be the first such case. Based on the 


event rate density arguments, these FRBs must only be 
a small fraction of all FRBs and may have some special 
properties. Another caveat is for individual cases, the 
robustness of the association must be addressed through 
various (e.g. temporal, spatial and distance) chance co- 
incidence probabilities as well as theoretical arguments 
(Moroianu et al., 2022). The new population may be es- 
tablished only after a sample of such association events 
are detected. 


B. Are there more than one class of repeating FRBs? 


This question actually has two aspects: First, observa- 
tionally, do we see different clustering properties among 
the observed repeaters? Second, physically, are there 
more than one type of engine sources that power dif- 
ferent repeaters? From the observational side, I’d argue 
that there have been already three types: 1. regular ac- 
tive repeaters in the cosmological distances (e.g. rFRB 
20121102A, rFRB 20180916B, rFRB 20190520B, rFRB 
20180301A, rFRB 20201124A), which have not been 
found in the Milky Way galaxy yet; 2. less energetic and 
less active magnetar repeaters such as SGR 1935+2154 
that produced two sub-pulses in FRB 20200428; and 
3. the globular cluster FRB 20200120E in M81, which 
has a high activity level but produces bursts with much 
lower luminosities than other cosmological active re- 
peaters (Nimmo et al., 2022). Since the central source of 
the second type is already identified as a magnetar, the 
general trend in the community is to attribute all three 
observationally identified types to magnetars, with differ- 
ent evolutionary stages and probably different formation 
channels as well. For example, the first type (active cos- 
mological repeaters) might be younger magnetars formed 
from recent supernova explosions and the third type may 
be magnetars produced from older formation channels 
such as WD-WD and NS-NS mergers or AIC of WDs. 
Even though this “magnetars make them all” hypothe- 
sis is theoretically attractive and passes some observa- 
tional constraints, it nonetheless suffers from some draw- 
backs. For example, the detection of FRB 20200120E 
from the M81 globular cluster suggests that these sys- 
tems are quite common (Kremer et al., 2021; Lu et al., 
2022). This seems to be inconsistent with the fact that 
none of the 30 discovered magnetars from the Milky Way 
or LMC/SMC are associated with globular clusters. The 
fact that the CHIME DM distribution demands a dom- 
inant delayed population of FRBs with respect to star 
formation Zhang and Zhang (2022), see also Hashimoto 
et al. (2022); Qiang et al. (2022), but see Shin et al. 
(2022)) also suggests that if magnetars do it all, the old- 
population magnetar channel should be the dominant 
one, in the contrary to the known magnetar population 
data. So, the current data may have already suggested 
the existence of other non-magnetar FRB engines. 
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C. FRB radiation mechanisms: where and how? 


Within the magnetar model of FRBs, there exist un- 
certainties regarding the location of the emission region 
(e.g. magnetospheres vs. relativistic shocks) and the 
radiation mechanism (bunched emission, plasma insta- 
bilities vs. vacuum maser mechanisms). As discussed 
in Sects. VLA, V.B and V.C, active studies and intense 
debates exist in the field, and growing evidence suggests 
that the magnetospheric origin is relevant for at least 
some FRBs. It remains unclear whether more than one 
emission site and more than one coherent mechanism are 
operating in FRBs. It is foreseen that the investigations 
in this direction will continue for years to come and the 
debates may not be settled in the near future, as the 
history of the study of the radiation mechanism of radio 
pulsars suggested. 


D. Prospects 


In a young and rapidly growing field, it is fun to make 
predictions. The three authors of Petroff et al. (2019), 
Emily Petroff, Jason Hessels, and Duncan Lorimer, made 
their respective predictions about the field in five years 
in the original review paper and also in their later up- 
dated review paper (Petroff et al., 2022). It is amusing 
to see that even though some of the predictions were real- 
ized, some unpredicted surprising discoveries were made 
within less than a three-year period since the first pre- 
dictions were made. These include a periodically modu- 
lated FRB (rFRB 20180916B) with a 16-d period, a MJy 
low-luminosity FRB (FRB 20200428) from a Galactic 
magnetar, and a repeating low-luminosity FRB (rFRB 
20200120E) from a globular cluster in M81. The FRB 
field seems to discourage conservative predictions. Just 
for fun, I close this review with ten predictions for the 
next 5-10 years. 


1. The detected FRB number will continue to grow 
rapidly, reaching ~ 10* different sources (includ- 
ing both non-repeating and repeating FRB pop- 
ulations) from the survey programs such as the 
CHIME FRB Project and reaching ~ 10^ bursts 
from a few active individual sources from dedicated 
observational campaigns such as the FAST FRB 
Key Project. 


2. The FRB community will continue to grow and the 
numbers of papers and citations per year will keep 
rising for another 5-10 years. 


3. Surprises will continue to arrive, which will shake 
the FRB theoretical framework a few times before 
a standard paradigm is established. 


4. X-ray counterparts of FRBs from nearby galaxies 
will be discovered, which are consistent with an 


SGR origin of FRBs. 


5. Despite active searches, prompt optical flashes in 
coincidence with FRBs will NOT be discovered, be- 
cause of the intrinsic faintness of the prompt optical 
emission. 


6. Claims about the associations between a progeni- 
tor of a magnetar (e.g. a long GRB, a superlumi- 
nous supernova, a short GRB, or a regular Type 
II supernova) and a repeating FRB source will be 
made, but a firm association cannot be established 
because of the uncertainties in chance coincidences. 


7. More claims about the associations between non- 
repeating FRBs and gravitational wave sources will 
be made, but the sample is not large and consistent 
enough to draw a definite conclusion. 


8. More Galactic FRBs will be detected, most likely 
from SGR. 1935+2154 or other magnetars, but also 
possible from sources other than magnetars, such 
as the Galactic center, young or old neutron stars, 
or even black hole binary systems. 


9. Multiple channels of repeating FRBs will be widely 
accepted. The ansatz that “all FRBs repeat” still 
cannot be completely ruled out. 


10. Debates on the physical mechanisms of FRBs will 
continue among theorists, not only because “a com- 
petent theorist can make any model to match any 
observational data”, but also because there might 
be indeed several physically plausible mechanisms 
that operate together. 
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